Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,5660TC38.5% 119.5 / 148.7 149E139E (GAA→GAG) ompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base T (41/50);  new base C (29/28);  total (70/79)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 153 total aligned reads displayed.

CCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCG  >  NC_000913/986468‑986662
                                                                                                  |                                                                                                
ccATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAAtt                                                                                                 >  1:2263174/1‑100 (MQ=255)
ccATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAAtt                                                                                                 >  2:2562050/1‑100 (MQ=255)
ccATCAACCAGACCAAAGAAGTCGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAAtt                                                                                                 >  1:2211088/1‑100 (MQ=255)
 cATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTc                                                                                                <  2:2542813/100‑1 (MQ=255)
  aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCt                                                                                               <  2:1546912/100‑1 (MQ=255)
   tCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTg                                                                                              <  2:666862/100‑1 (MQ=255)
   tCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTg                                                                                              <  1:628115/100‑1 (MQ=255)
   tCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTg                                                                                              <  2:1640860/100‑1 (MQ=255)
   tCAACCAGACCAAAGAAGTTGGAGTTACGAGAGGTAGCAACGCCGCCAACGCGACCAACGAAGAAGTCAGCGCTGTATGCAGTATCACCACCAAATTCTg                                                                                              <  1:23240/100‑1 (MQ=255)
     aaCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgc                                                                                            <  2:78577/100‑1 (MQ=255)
       ccAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcag                                                                                          >  1:2182259/1‑100 (MQ=255)
       ccAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcag                                                                                          >  2:2047758/1‑100 (MQ=255)
        cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagc                                                                                         <  1:947140/100‑1 (MQ=255)
        cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagc                                                                                         <  1:2219079/100‑1 (MQ=255)
         aGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGgcagca                                                                                        <  2:184754/100‑1 (MQ=255)
         aGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagca                                                                                        <  2:2310366/100‑1 (MQ=255)
            ccAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCatat                                                                                     >  1:2254619/1‑100 (MQ=255)
            ccAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCatat                                                                                     >  1:1025916/1‑100 (MQ=255)
              aaaGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCg                                                                                   >  1:2376687/1‑100 (MQ=255)
              aaaGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCg                                                                                   >  1:617900/1‑100 (MQ=255)
               aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                                  <  1:2376512/100‑1 (MQ=255)
               aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                                  <  1:871177/100‑1 (MQ=255)
               aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                                  <  2:1167459/100‑1 (MQ=255)
                 gaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGgtg                                                                                <  2:375956/100‑1 (MQ=255)
                   agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTa                                                                              <  2:919445/100‑1 (MQ=38)
                    gTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTaa                                                                             <  2:1644077/100‑1 (MQ=255)
                      tGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAAcc                                                                           >  1:2074744/1‑100 (MQ=38)
                      cggAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAAcc                                                                           >  1:1710895/2‑100 (MQ=38)
                          gTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGt                                                                       >  1:512209/1‑100 (MQ=38)
                           ttACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTg                                                                      <  2:1663794/100‑1 (MQ=255)
                            tACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGc                                                                     >  1:415469/1‑100 (MQ=255)
                             aCGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGca                                                                    <  2:1465256/100‑1 (MQ=38)
                              cGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcat                                                                   <  1:1881806/100‑1 (MQ=255)
                               gATAGGTAGCAACGCCGCCAACACGACCAACGAAGTAGTCGTCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatc                                                                  >  2:2532933/1‑100 (MQ=255)
                                aTAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGcatca                                                                 >  1:1167522/1‑100 (MQ=38)
                                aTAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatca                                                                 <  1:2123115/100‑1 (MQ=255)
                                   ggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaa                                                              <  2:544030/100‑1 (MQ=255)
                                   ggTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaa                                                              >  1:845596/1‑100 (MQ=255)
                                    gTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATaaa                                                             >  1:75202/1‑100 (MQ=38)
                                     tAGCAACGCCGCCAACACGACCAATGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaaan                                                            <  1:711453/100‑2 (MQ=255)
                                      ggCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                           >  1:370123/2‑100 (MQ=38)
                                      aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                           <  2:682629/100‑1 (MQ=38)
                                      aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                           <  2:603457/100‑1 (MQ=38)
                                      aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                           <  2:207671/100‑1 (MQ=38)
                                        cAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACcac                                                         >  1:1795486/1‑100 (MQ=255)
                                        cAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACcac                                                         <  2:918968/100‑1 (MQ=255)
                                         naCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAGCACCACCAAATTCTGGCAGCATATCGGGGTAACCCAGTGCATCAGAAACcaca                                                        >  1:2113983/2‑100 (MQ=255)
                                         aaCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACcaca                                                        >  2:1786649/1‑100 (MQ=255)
                                          aCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCacac                                                       <  2:2580217/100‑1 (MQ=255)
                                           cgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACAcc                                                      >  1:2386289/1‑100 (MQ=38)
                                           cgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACAcc                                                      >  1:1484546/1‑100 (MQ=38)
                                            gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                                     >  1:1955816/1‑100 (MQ=38)
                                            gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                                     >  1:1702312/1‑100 (MQ=38)
                                            gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                                     >  1:445487/1‑100 (MQ=38)
                                            gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                                     >  1:449368/1‑100 (MQ=38)
                                            gccgccAACACGACCAACGAAGAAGTCATCGCTGTATACCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                                     <  2:2355914/100‑1 (MQ=37)
                                              cgccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTa                                                   <  2:1687743/100‑1 (MQ=38)
                                               gccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                                  <  2:2311402/100‑1 (MQ=38)
                                               gccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                                  <  2:1487797/100‑1 (MQ=38)
                                               gccAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                                  <  2:1529888/100‑1 (MQ=38)
                                               gccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAg                                                  <  1:2311597/100‑1 (MQ=255)
                                                 cAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGtt                                                <  2:2226196/100‑1 (MQ=38)
                                                  aaCACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTa                                               <  2:2502603/100‑1 (MQ=255)
                                                       gACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGcc                                          >  1:64962/1‑100 (MQ=255)
                                                         ccAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGt                                        >  2:2424157/1‑100 (MQ=255)
                                                         ccAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGt                                        >  2:322852/1‑100 (MQ=255)
                                                         ccAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGt                                        <  2:1320542/100‑1 (MQ=255)
                                                             cGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATc                                    >  2:1404601/1‑100 (MQ=255)
                                                             cGAAGAAGGCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATc                                    <  2:1321000/100‑1 (MQ=255)
                                                              gaagaaGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCg                                   <  2:2341143/100‑1 (MQ=38)
                                                                agaagTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTACCCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGag                                 <  2:1374350/100‑2 (MQ=38)
                                                                agaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaa                                 <  2:1598767/100‑1 (MQ=255)
                                                                 taaatCATCGCTGTAGGCAGTCTCCCCTCCAAATTCTAGCAGCAGGTCGGCGTAACCCAGTGCATCATAATCCCCACCGTAGTTACGGCCGTAATCGaaa                                <  2:2210630/96‑1 (MQ=255)
                                                                 gaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaaa                                <  2:1979841/100‑1 (MQ=255)
                                                                  aagTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTCACCCAGTGCATCATAAACCATACCGTAGTTACTGCCGTAATCGAAAg                               <  2:2428632/100‑1 (MQ=37)
                                                                     tCATCGCGGTATGCCGTGTCACCGCCACACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGTCGTAATCGAAAGAAc                            <  2:1974471/100‑1 (MQ=255)
                                                                       aTCGCTGTATTCAGTATCACCACCAAAGGCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACGCCGTAGTTACGGCCGTAATCGAAAGAACCa                          <  2:1417284/100‑1 (MQ=255)
                                                                         cGCTGTATGCAGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAAc                        >  1:510220/1‑100 (MQ=39)
                                                                          gCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTCCGGCCGTAATCGAAAGAACCAACg                       >  1:1169936/1‑100 (MQ=255)
                                                                           cTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGt                      <  1:1921459/100‑1 (MQ=255)
                                                                             gTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCa                    >  2:1847328/1‑100 (MQ=255)
                                                                              tATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAg                   <  1:2495391/100‑1 (MQ=255)
                                                                               aTGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc                  >  1:998172/1‑100 (MQ=38)
                                                                               aTGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGc                  >  1:1988973/1‑100 (MQ=38)
                                                                               aTGCCGTGTCACCGCCAAACGCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCAAACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAgt                  <  2:1058746/100‑2 (MQ=255)
                                                                                tGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCg                 >  1:1960501/1‑100 (MQ=38)
                                                                                tGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCg                 <  1:1478391/100‑1 (MQ=255)
                                                                                    gTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTAtt             >  1:1009060/1‑100 (MQ=255)
                                                                                    gTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTAtt             <  1:930731/100‑1 (MQ=255)
                                                                                        caccgcCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAg         >  1:691436/1‑100 (MQ=255)
                                                                                         accaccAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGa        >  2:585856/1‑100 (MQ=255)
                                                                                          ccgccaAACTCAGGCAGCATATCGGTGTAACCAAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAc       <  2:2046814/97‑1 (MQ=255)
                                                                                          ccaccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAc       >  1:200388/1‑100 (MQ=255)
                                                                                           ngccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc      >  1:302758/3‑100 (MQ=255)
                                                                                           cgccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc      >  1:1708334/3‑100 (MQ=255)
                                                                                           cgccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc      >  1:78284/3‑100 (MQ=255)
                                                                                            accaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAccc     <  1:2161553/100‑1 (MQ=255)
                                                                                             ccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCcg    <  2:2147981/100‑1 (MQ=255)
                                                                                             ccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCcg    >  1:682141/1‑100 (MQ=255)
                                                                                               aaaTTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCcgcg  <  2:936475/100‑1 (MQ=255)
                                                                                                  |                                                                                                
CCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCG  >  NC_000913/986468‑986662

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