Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,745,627 C→T 5.6% G167G (GGG→GGA)  ydhQ ← autotransporter adhesin‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,745,6270CT 5.6% 228.4 / 2.4 71G167G (GGG→GGA) ydhQautotransporter adhesin‑related protein
Reads supporting (aligned to +/- strand):  ref base C (33/34);  new base T (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

ACGATAAATGGATCCATTCTCAACAATAGTGTCAGAGTCTTCACCACCGAGATTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGATTTTTTCCGCGCGATGTCCTTCCAGTACACGCAGGTTACCGCCATTTTCCAG  >  NC_000913/1745483‑1745773
                                                                                                                                                |                                                                                                                                                  
aCGATAAATGGATCCATTCTCAACAATAGTGTCAGAGTCTTCACCACCGAGATTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCa                                                                                                                                                 <  8:406176/148‑1 (MQ=255)
      aaTGGATCCATTCTCAACAATAGTGTCAGAGTCTTCACCACCGAGATTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGa                                                                                                                                           <  2:43933/148‑1 (MQ=255)
       aTGGATCCATTCTCAACAATAGTGTCAGAGTCTTCACCACCGAGATTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGAt                                                                                                                                          <  4:300742/148‑1 (MQ=255)
       aTGGATCCATTCTCAACAATAGTGTCAGAGTCTTCACCACCGAGATTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGAt                                                                                                                                          <  3:252831/148‑1 (MQ=255)
         ggATCCATTCTCAACAATAGTGTCAGAGTCTTCACCACCGAGATTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATa                                                                                                                                         >  8:354486/1‑147 (MQ=255)
              cATTCTCAACAATAGTGTCAGAGTCTTCACCACCGAGATTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACca                                                                                                                                    >  6:263059/1‑147 (MQ=255)
                    caacaaTAGTGTCAGAGTCTTCACCACCGAGATTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCCACAAtt                                                                                                                              >  5:288476/1‑147 (MQ=255)
                    caacaaTAGTGTCAGAGTCTTCACCACCGAGATTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAAtt                                                                                                                              >  3:385017/1‑147 (MQ=255)
                     aacaaTAGTGTCAGAGTCTTCACCACCGAGATTATTGATGGTGCCACCAGCAAGCAACGTATAACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAAt                                                                                                                               <  7:338146/145‑1 (MQ=255)
                       aaaTAGTGTCAGAGTAGTCACCACCGAGATTATTCATGGTGCCACCAGCAAGCACCGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTcacc                                                                                                                          <  6:288476/147‑1 (MQ=255)
                        aaTAGTGTCAGAGTCTTCACCACCGAGATTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTcacca                                                                                                                         <  4:385017/148‑1 (MQ=255)
                          tAGTGTCAGAGTCTTCACCACCGAGATTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACATTTGCTGGCTTCCCCAcc                                                                                                                                               >  6:289284/1‑124 (MQ=255)
                          tAGTGTCAGAGTCTTCACCACCGAGATTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACATTTGCTGGCTTCCCCAcc                                                                                                                                               <  5:289284/124‑1 (MQ=255)
                          tAGTGTCAGAGTCTTCACCACCGAGATTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCaccacc                                                                                                                       <  8:222626/148‑1 (MQ=255)
                                     tCTTCACCACCGAGATTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCa                                                                                                            <  3:110714/148‑1 (MQ=255)
                                      cTTCACCACCGAGATTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACTCCGCCCTGATTAATCTCACAATTGCTTGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCAc                                                                                                           >  7:438593/1‑148 (MQ=255)
                                           ccaccGAGATTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCg                                                                                                      >  4:279862/1‑148 (MQ=255)
                                             accGAGATTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCt                                                                                                    >  3:255734/1‑148 (MQ=255)
                                                gagaTTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAg                                                                                                 <  8:307664/148‑1 (MQ=255)
                                                 agaTTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAgg                                                                                                <  8:106248/148‑1 (MQ=255)
                                                  gaTTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGt                                                                                               >  4:135279/1‑148 (MQ=255)
                                                       ttCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGCTACACCAACAATTCACCACCTTCATCTACAACGCCCGCTGAGGTTGTcc                                                                                          >  5:65940/1‑148 (MQ=255)
                                                          aTGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCAt                                                                                       <  7:354486/148‑1 (MQ=255)
                                                           tGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCAt                                                                                       >  2:411471/1‑147 (MQ=255)
                                                           tGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCAt                                                                                       >  4:128760/1‑147 (MQ=255)
                                                                ccaccaGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCa                                                                                 <  4:255734/148‑1 (MQ=255)
                                                                     agcaagcaACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCt                                                                                                    >  4:410193/1‑124 (MQ=255)
                                                                     agcaagcaACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCt                                                                                                    <  3:410193/124‑1 (MQ=255)
                                                                      gcaagcaaCGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGgc                                                                           >  3:70640/1‑148 (MQ=255)
                                                                        aagcaaCGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGAt                                                                                                                                          <  2:296469/83‑1 (MQ=255)
                                                                        aagcaaCGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGAt                                                                                                                                          >  1:296469/1‑83 (MQ=255)
                                                                          gcaaCGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGgccgcc                                                                       >  4:222573/1‑148 (MQ=255)
                                                                          gcaaCGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGgccgcc                                                                       <  3:279862/148‑1 (MQ=255)
                                                                                 aTCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGa                                                                >  4:99093/1‑148 (MQ=255)
                                                                                  tCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGa                                                                <  1:411471/147‑1 (MQ=255)
                                                                                  tCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTAACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGa                                                                <  4:70640/147‑1 (MQ=255)
                                                                                   cACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCAtt                                                                                      >  3:365190/1‑124 (MQ=255)
                                                                                   cACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCAtt                                                                                      <  4:365190/124‑1 (MQ=255)
                                                                                    aCTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCg                                                             <  6:98573/148‑1 (MQ=255)
                                                                                      tGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACTATTCACCACCTTCATCTACCAGGACCGCTGAg                                                                                                 <  1:100044/110‑1 (MQ=255)
                                                                                       ggCTTTCCCGCCCAGCATAAAAACGCCGCCCTGATAAAACTCACAATTGCTGGCTTCCCCACCTCGATACACCAACAATTCACCACCATCATCTACCACGACCGCTGAgg                                                                                                >  2:100044/1‑110 (MQ=255)
                                                                                          tttCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAAt                                                        >  1:274607/1‑147 (MQ=255)
                                                                                           ttCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGa                                                      >  2:159762/1‑148 (MQ=255)
                                                                                            tCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGACCGAGAATg                                                       <  5:263059/146‑1 (MQ=255)
                                                                                              ccGGACAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGATGGCGTACCCAGCTGGATACACCAACAATTAACCACCTTCATGTTCCACGACCGCTGAGGTTGTCCCATTAAACAACAGGCCGCCctct                                                                   <  6:242878/132‑1 (MQ=255)
                                                                                               cGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCTCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGAtttt                                                  >  7:254585/1‑148 (MQ=255)
                                                                                               cGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCTCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGAtttt                                                  >  7:248116/1‑148 (MQ=255)
                                                                                               cGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTACCCACCTGGATACACCAACAATACACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCtt                                                                  >  5:242878/1‑132 (MQ=255)
                                                                                                        tAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGATTTTTTCcgcgcg                                         <  3:222573/148‑1 (MQ=255)
                                                                                                        tAAAAACGCCGCCCAGATTAATCTCACAAATGCTTGCTTCCCGACCTGTATACCACAACAATTCACAACCTTCATCTACCACGACAGCTGATGTTGTCCCATTAACCAACAGGCCGCCCTCTTGa                                                                <  5:109352/125‑1 (MQ=255)
                                                                                                        tAAAAACGACGCCCTGATTAATCTCACAATTGCTGTCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGt                                                                                               <  8:418001/94‑1 (MQ=255)
                                                                                                         aaaaaCGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAAc                                                                               >  8:401739/1‑109 (MQ=255)
                                                                                                         aaaaaCGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAAc                                                                               <  7:401739/109‑1 (MQ=255)
                                                                                                         aaaaaCGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGATTTTTTCCGCGCGa                                        >  5:291197/1‑148 (MQ=255)
                                                                                                         aaaaaCGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGATTTTTTCCGCGCGa                                        >  3:54025/1‑148 (MQ=255)
                                                                                                         aaaaaCGCCGCCCTGATTAATCTCACAATTCCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGtt                                                                                              >  7:418001/1‑94 (MQ=255)
                                                                                                          aaaaCGCCGCCCTGATTAATCTCACAATTGCTGGCTTCTCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGATTTTTTCCGCGCGa                                        <  8:254585/147‑1 (MQ=255)
                                                                                                          aaaaCGCCGCCCTGATTAATCTCACAATTGCTGGCTTCTCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGATTTTTTCCGCGCGa                                        <  8:248116/147‑1 (MQ=255)
                                                                                                             aCGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCCCCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGATTTTTTCAGCGCGATGTc                                    >  3:331985/1‑148 (MQ=255)
                                                                                                                      tGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGATTTTTTCcg                                             <  5:431152/130‑1 (MQ=255)
                                                                                                                      tGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGATTTTTTCcg                                             >  6:431152/1‑130 (MQ=255)
                                                                                                                      tGATTAATCTCACAATTGCGGGCTACCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCAGCCCTCTTGATCGAGAATGATTTTTTCcgcgcg                                         <  4:222260/134‑1 (MQ=255)
                                                                                                                       gATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGATTTTTTCCGCGCGa                                        >  2:422244/1‑134 (MQ=255)
                                                                                                                       gATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGATTTTTTCCGCGCGa                                        <  1:422244/134‑1 (MQ=255)
                                                                                                                         ttAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGATTTTTTCCGCGCGATGTCCTTCCAGTacac                        >  4:79578/1‑148 (MQ=255)
                                                                                                                            aTCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGATTTTTTCCGCGCGATGTCCTTCCAGTACACGCa                     <  3:128760/148‑1 (MQ=255)
                                                                                                                              ctcACAATTGCAGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGATTTTTTCCGCGCGATGTACTTCCAGTACACGCAgg                   >  1:327637/1‑148 (MQ=255)
                                                                                                                                 acaATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGATTTTTTCCGCGCGATGTCCTTCCAGTACACGCAGGTTa                >  7:181657/1‑148 (MQ=255)
                                                                                                                                    aTTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCa                                                                                 <  2:94966/80‑1 (MQ=255)
                                                                                                                                    aTTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCa                                                                                 >  1:94966/1‑80 (MQ=255)
                                                                                                                                           gCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCTAGAATGATTTTTTCCGCGCGATGTCCTTCCAGTACACGCAGGTTACCGCCAtttt      <  3:135279/148‑1 (MQ=255)
                                                                                                                                                cccACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGATTTTTTCCGCGCGATGTCCTTCCAGTACACGCAGGTTACCGCCATTTTCCAg  <  3:79578/147‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
ACGATAAATGGATCCATTCTCAACAATAGTGTCAGAGTCTTCACCACCGAGATTATTCATGGTGCCACCAGCAAGCAACGTATCACTGGCTTTCCCGGCCAGCATAAAAACGCCGCCCTGATTAATCTCACAATTGCTGGCTTCCCCACCTGGATACACCAACAATTCACCACCTTCATCTACCACGACCGCTGAGGTTGTCCCATTAACCAACAGGCCGCCCTCTTGATCGAGAATGATTTTTTCCGCGCGATGTCCTTCCAGTACACGCAGGTTACCGCCATTTTCCAG  >  NC_000913/1745483‑1745773

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑