Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,5660TC27.4% 399.5 / 190.0 261E139E (GAA→GAG) ompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base T (105/83);  new base C (39/32);  total (145/116)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.12e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 269 total aligned reads displayed.

GCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAAT  >  NC_000913/986467‑986665
                                                                                                   |                                                                                                   
gCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCGACACGACCAACGCAGAGGTCATCGCTGTATGCAGTATCACGACCAAAt                                                                                                     >  1:1540878/1‑100 (MQ=38)
gCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAAt                                                                                                     <  2:3633703/100‑1 (MQ=255)
gCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAAt                                                                                                     <  2:3761792/100‑1 (MQ=255)
  cATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTc                                                                                                   >  2:2207286/1‑100 (MQ=255)
   aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCt                                                                                                  <  2:2694457/100‑1 (MQ=255)
   aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTAGGCAGTATCACCACCAAATTCt                                                                                                  >  1:4115248/1‑100 (MQ=255)
    tCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTg                                                                                                 <  1:1249750/100‑1 (MQ=255)
      aaCCATACCAAGGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGGCATCGCTGTATGCAGTATCACCACCAAATTCTGgc                                                                                               <  2:396523/100‑1 (MQ=255)
      aaCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgc                                                                                               <  2:894054/100‑1 (MQ=255)
      aaCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgc                                                                                               <  2:83770/100‑1 (MQ=255)
      aaCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgc                                                                                               >  1:3908761/1‑100 (MQ=255)
      aaCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgc                                                                                               >  2:2566546/1‑100 (MQ=255)
       aCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgca                                                                                              >  2:1225594/1‑100 (MQ=255)
       aCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgca                                                                                              <  1:91356/100‑1 (MQ=255)
        ccAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGgcag                                                                                             <  2:482994/100‑1 (MQ=255)
        ccAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcag                                                                                             >  2:3384530/1‑100 (MQ=255)
         cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGgcagc                                                                                            <  2:1725507/100‑1 (MQ=255)
           gACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCat                                                                                          >  1:2775652/1‑100 (MQ=255)
           gACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGACATCGCTGTAAGCAGTATCACCACCAAATTCTGGCACCat                                                                                          >  2:1318839/1‑100 (MQ=255)
             ccAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCatat                                                                                        >  2:1446347/1‑100 (MQ=255)
              cAACGAAGTTGGCGCTACCATAGGTAGCAACGCCGCCAACCCGACCAACGAACAAGTCATCACTCTATGCGGTAACACCCCCAAATTCTGGCAGCatcaa                                                                                       >  2:4088883/1‑97 (MQ=255)
                aagaagTTGGAGTTACGATAGGTAGCAACGCCGGCAACACGACGAACGAAGAAGTCATCGCTGTATGCAGTATAACCACCAAATTCTGGCAGCATATTgg                                                                                     >  2:3586332/1‑100 (MQ=255)
                aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCgg                                                                                     <  2:2107207/100‑1 (MQ=255)
                aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCgg                                                                                     <  2:604220/100‑1 (MQ=255)
                aagaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCgg                                                                                     <  2:2043379/100‑1 (MQ=255)
                aagaagTGGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCgg                                                                                     <  2:3042634/100‑1 (MQ=255)
                  aaaattGGAGGTACGATAGGTAACAACGCCGCCAAGACGACCAACGAAGAAGTCAGGGCGGTATGCAGTATCACCACCAAATGCTGGCAGCAGATCtgtg                                                                                   <  2:1057067/96‑1 (MQ=255)
                   aagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGgtgt                                                                                  <  1:3458215/100‑1 (MQ=255)
                     gTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTaa                                                                                >  2:3101815/1‑100 (MQ=255)
                        ggAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAAccc                                                                             >  1:1705010/1‑100 (MQ=38)
                        ggAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAAccc                                                                             <  2:2036074/100‑1 (MQ=255)
                        ggAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAAccc                                                                             <  1:689785/100‑1 (MQ=255)
                         gAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCa                                                                            >  1:3623471/1‑100 (MQ=255)
                          aGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAg                                                                           <  1:2049210/100‑1 (MQ=255)
                           gTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGt                                                                          <  2:3007315/100‑1 (MQ=255)
                           gTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGt                                                                          >  2:4318845/1‑100 (MQ=255)
                            ttACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTg                                                                         <  1:2492369/100‑1 (MQ=255)
                               cGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcat                                                                      >  1:2998672/1‑100 (MQ=255)
                                gATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGcatc                                                                     <  2:788292/100‑1 (MQ=38)
                                gATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatc                                                                     >  2:3887665/1‑100 (MQ=255)
                                 aTAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGcatca                                                                    >  1:2026366/1‑100 (MQ=38)
                                 aTAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCAGCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCACTGcatca                                                                    >  2:2485505/1‑100 (MQ=255)
                                  tAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGcatcat                                                                   <  1:512433/100‑1 (MQ=255)
                                     gTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATaaa                                                                >  1:186087/1‑100 (MQ=255)
                                      tAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAc                                                               <  2:2507921/100‑1 (MQ=38)
                                      tAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAc                                                               >  2:79397/1‑100 (MQ=255)
                                       aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                              <  2:3772439/100‑1 (MQ=38)
                                       aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                              <  2:2228267/100‑1 (MQ=38)
                                       aGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcc                                                              <  1:3865682/100‑1 (MQ=255)
                                         cAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAAcctc                                                            <  2:237361/100‑3 (MQ=38)
                                           ccgccgcCAACACGACCAACGAAGAAGTCATCCCTGTAGGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCacac                                                          >  1:3920616/2‑100 (MQ=37)
                                           aCGCCGCCAACACGCCCATCTAAGAAGGCAGGGCTGTAGGCACTATCACCACCAACCTCCCGCAACATATCGGTTTAACCCAGTGCATCATAAAcccgac                                                          <  2:498971/100‑5 (MQ=255)
                                            cgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACAcc                                                         >  1:3389422/1‑100 (MQ=38)
                                            cgccgcCAACACGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACAcc                                                         >  1:833185/1‑100 (MQ=38)
                                             gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCg                                                        >  1:3737736/1‑100 (MQ=255)
                                              ccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGt                                                       >  1:3144017/1‑100 (MQ=255)
                                                 caaaCACAACAAACGAAGAAGTCATCACTGTATGCAGTATAACCAACCAATACTAACAGCAAAACAGTGTAACCCAGTGCATCATAAACCACACCCTata                                                    >  2:2617885/3‑98 (MQ=255)
                                                  cAACACGACCAACGAAGAAGTCATCGCTGTACGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGtt                                                   >  2:1804570/1‑100 (MQ=255)
                                                   aaCACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTa                                                  <  1:1957861/100‑1 (MQ=255)
                                                     cacGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACg                                                >  2:1115381/1‑100 (MQ=255)
                                                     cacGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACg                                                <  1:2603016/100‑1 (MQ=255)
                                                     cacGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACg                                                <  1:3646145/100‑1 (MQ=255)
                                                      gcggccAAAGGAGAAGTCATGCGTGCATGCCGTGTCACCGCCAAACTCAGGCGGCATTTCGGTGCAACGCAGTGCATCATAAACCACACCGTAGTTACgg                                               <  2:2542190/96‑1 (MQ=255)
                                                      acGACCAACGAAGAAGTCATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACgg                                               <  2:1996829/100‑1 (MQ=38)
                                                       cGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGc                                              >  1:2195086/1‑100 (MQ=255)
                                                         aCCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCg                                            >  1:3276252/1‑100 (MQ=255)
                                                         aCCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCg                                            <  1:1690529/100‑1 (MQ=255)
                                                          ncAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGt                                           >  1:1028911/2‑100 (MQ=255)
                                                           cAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTa                                          >  2:1664152/1‑100 (MQ=255)
                                                            aaCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTaa                                         >  1:1895015/1‑100 (MQ=255)
                                                             aCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAAt                                        >  1:3073941/1‑100 (MQ=255)
                                                                aagaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGa                                     <  2:3831079/100‑1 (MQ=255)
                                                                 agaagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGaa                                    <  2:1380466/100‑1 (MQ=255)
                                                                   aagTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAg                                  <  1:2685180/100‑1 (MQ=255)
                                                                      tCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAAc                               <  2:3231447/100‑1 (MQ=255)
                                                                       cATCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAAcc                              >  1:4407990/1‑100 (MQ=38)
                                                                        aTCGCTGTATGCCGTGTCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCa                             >  1:2402084/1‑100 (MQ=38)
                                                                        aTCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCa                             >  1:1845839/1‑100 (MQ=255)
                                                                         tCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCaa                            >  1:1625675/1‑100 (MQ=255)
                                                                           gCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACg                          <  1:651386/100‑1 (MQ=255)
                                                                           gCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACg                          <  1:2412869/100‑1 (MQ=255)
                                                                           gCTGTATGCAGTATCACCACCAAATTCTGGCAGCATACCGGTGTAACCCAGTGCATCATAAACCACACCCCAGTTACGGCCGTAATCGAAAGAACCAAcc                          >  2:4292884/1‑99 (MQ=255)
                                                                               tATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGTCGTAATCGAAAGAACCAACGTCAg                      >  2:3189999/1‑100 (MQ=255)
                                                                                 tGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACTCCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCg                    <  2:1516338/100‑1 (MQ=255)
                                                                                 tGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCg                    >  2:1889170/1‑100 (MQ=255)
                                                                                  gCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGt                   >  1:435770/1‑100 (MQ=255)
                                                                                   cAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTa                  >  1:3229656/1‑100 (MQ=255)
                                                                                   cAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTa                  <  2:4208432/100‑1 (MQ=255)
                                                                                     ntgtCACCGCCAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGGAGTTACGGCCTTAATCGAAAGAACCAACGTCAGCGTAtt                >  1:3531483/4‑100 (MQ=38)
                                                                                       aTCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTa              >  1:3717248/1‑100 (MQ=255)
                                                                                       aTCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTa              <  2:570647/100‑1 (MQ=255)
                                                                                          accgccAAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGa           >  1:2444704/1‑100 (MQ=255)
                                                                                            cgccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAcc         >  1:1416419/3‑100 (MQ=255)
                                                                                             gccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAccc        >  1:2571693/2‑100 (MQ=255)
                                                                                             gccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAccc        >  1:799825/2‑100 (MQ=255)
                                                                                             accaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAccc        >  1:3830821/1‑100 (MQ=255)
                                                                                             accaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGAccc        >  1:1316253/1‑100 (MQ=255)
                                                                                              ccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCcg       >  1:1157384/1‑100 (MQ=255)
                                                                                              ccaAACTCAGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCcg       <  2:2770367/100‑1 (MQ=255)
                                                                                                  aTTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATGCAAAGAACCAACGTCAGGGTATTTAAGAcgcggtaa   >  2:3700678/1‑96 (MQ=255)
                                                                                                   |                                                                                                   
GCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAAT  >  NC_000913/986467‑986665

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: ‑