Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,349,1450AC40.3% 120.8 / 16.9 67F110V (TTC→GTC) inaAacid‑inducible Kdo/WaaP family putative kinase
Reads supporting (aligned to +/- strand):  ref base A (35/5);  new base C (0/27);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 4.91e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATATGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCG  >  NC_000913/2348995‑2349286
                                                                                                                                                        |                                                                                                                                             
aGATATGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGaa                                                                                                                                               <  2:1028390/151‑1 (MQ=255)
  atatGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAAcc                                                                                                                                             <  2:1203130/151‑1 (MQ=255)
    atGGCGAACATAACAACAGCTATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAAcccc                                                                                                                                           >  1:255568/1‑151 (MQ=255)
    atGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAAcccc                                                                                                                                           >  1:894994/1‑151 (MQ=255)
    atGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAAcccc                                                                                                                                           <  1:788924/151‑1 (MQ=255)
    atGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAAccc                                                                                                                                            <  2:351439/150‑1 (MQ=255)
     tGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCg                                                                                                                                          >  2:257347/1‑151 (MQ=255)
        cGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGcc                                                                                                                                        >  1:415819/1‑150 (MQ=255)
         gAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCat                                                                                                                                      >  1:399317/1‑151 (MQ=255)
          aaCATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCat                                                                                                                                      >  2:417309/1‑150 (MQ=255)
                  aacaGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGt                                                                                                                             >  2:334058/1‑151 (MQ=255)
                      gCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCAcc                                                                                                                         >  2:1253877/1‑151 (MQ=255)
                       ccATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCAcc                                                                                                                         >  2:657045/1‑150 (MQ=255)
                        cATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCAcc                                                                                                                         >  1:944781/1‑149 (MQ=255)
                         aTGCTGACTATTAATGTTATGCTTTTTCTAAACCGCCAGCGCCTCTGCTTTAAACATGGCTTCCCGTTTTTCGTCAGAAAAGGGCGATACTGCATGCGGGGCATCCCAGTCAGCAATGCTGATAACCCCCGCCATATCTTCAGTCACCa                                                                                                                        <  2:1290176/149‑1 (MQ=255)
                          tGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGc                                                                                                                      >  2:581573/1‑150 (MQ=255)
                           gCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCa                                                                                                                     >  2:942837/1‑150 (MQ=255)
                           gCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCa                                                                                                                     >  2:1264358/1‑150 (MQ=255)
                            cTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCa                                                                                                                               >  2:35969/1‑139 (MQ=255)
                            cTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAAc                                                                                                                   >  2:650487/1‑151 (MQ=255)
                              gACTATTATTTCTCTGCTTTTTTTTAAACGCCAGCGCCCTTGCTTTCAACATGGCTTGCCGTATTTCGTCAAAAAAAGGCGTTTCTGCCTGCTGGGCATCCTAGCCAGCAATGCTGAAGACCCCCGCCATATCTTCAGTCACCAGCAAc                                                                                                                   <  2:927405/149‑1 (MQ=255)
                               aCGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCt                                                                                                                >  2:316448/1‑151 (MQ=255)
                                cgtttAATGCTTTGTTTTTTCTTAACCGCCGCCCCCACTTTTTCCAACATGGCTTGCCGTTCTTCGTAAGAAAAAGGCGATACGGCATGCTGGGCATCCCACTCAGCAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCt                                                                                                                <  2:108380/147‑1 (MQ=255)
                                cGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTc                                                                                                               >  1:749311/1‑151 (MQ=255)
                                cGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTc                                                                                                               >  1:772911/1‑151 (MQ=255)
                                  aTTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGc                                                                                                             >  2:152624/1‑151 (MQ=255)
                                   ttaattttatcttttttCTAAAACGCCGGCCCCATTGCTTACAAAATGGTTGGCCGTCCTTCGTAAAAATAGGGCGATCCGGCATGCTGGGCATCCCAGTCAGAAATGTTGATGACCCCCGCCATATCTTCAGTCACCAGCAACg                                                                                                                  <  1:1282408/132‑1 (MQ=255)
                                   ttAATGCTCTGTTTTTTCTAAACGGCCGGCGCCGCTGCTTACATCTTGCTTGGCCGTGGTGCGCAAAAACAAGCCGATACTGCGTGCTGGGCAGACCACTCAAAATTGCGGAGGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCg                                                                                                              <  2:390205/149‑1 (MQ=255)
                                     aaTGTTATGCTTTTTCTAAACCGCCGGCCCCACTTTTTTAAACATGTCTTCCCTTCCTTCGTAAAAAAAGGGCGATCCTGCATGCGGGCCATCCCAGTCAGAAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCa                                                                                                                     <  1:102637/140‑1 (MQ=255)
                                        gCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTc                                                                                                        >  1:653308/1‑150 (MQ=255)
                                          tATCCATTTTCTAAAACCCCGGCGCCACTGCTTTAAACATGGCTGGCCTTCCTTCTTCAGAAAAAGGAGATACTGCTTGCTGGGCATACCAGTCAGAAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTc                                                                                                        <  2:291314/148‑1 (MQ=255)
                                                 tttCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAAtt                                                                                              >  1:180414/1‑151 (MQ=255)
                                                  ttCTAAAACGCCGGCCCCACTGCTTTAAACAGGGCTTGCCGTACTTCGACAAAATAAGGCGATACTGCAGGCGGGCCATACCAGTCAGCAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAAtt                                                                                              <  1:1175065/150‑1 (MQ=255)
                                                     ttAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTc                                                                                           >  1:1275619/1‑150 (MQ=255)
                                                     ttAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTc                                                                                           >  1:786350/1‑150 (MQ=255)
                                                       aaaCGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCAcc                                                                                        >  2:270830/1‑151 (MQ=255)
                                                        acccccAGCCCCACTTCTTTCAACATGGCTTGCCGTACTTCGTAAGAATAAGGCGATCCTGCAGGCTGGGCATCCCAGTCAGCAATGCTAATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCg                                                                                       <  1:166238/147‑1 (MQ=255)
                                                           gCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCtt                                                                                    >  1:776546/1‑151 (MQ=255)
                                                            ccggccccACTGTTTACAAAAGGGTTTGCCTTCCTGCGAAAAAAAAAGGCGATACGGCAGGGGGGGACCACCAGTCAGCATGCCTGAGCACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCg                                                                                       <  1:125806/141‑1 (MQ=255)
                                                            ccAGCCCCTCTGTTTTCAACATGGCTTGCCGTTCTTCTTAAAAATAGGGCGATACTGCATGCGGGGCATCCCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTc                                                                                   <  2:257735/151‑1 (MQ=255)
                                                                cgcCCCTGTTTTCAACATGGCTTGCCGTCCTCCGTCAAAAAAAGGCGATACTGCATGCGGGGAATCCCAGTCAGCATTGCGGATAAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCAc                                                                                 <  2:603107/149‑1 (MQ=255)
                                                                cgcCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCAcc                                                                                >  2:424865/1‑150 (MQ=255)
                                                                   cccctgtttCAACCATGGCTTGCCGTACTTCGTAAAAATAAGGCATTCCTGCTTGCGGGGCACCCCAGCCAACATTGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCaa                                                                              <  1:672779/144‑1 (MQ=255)
                                                                   cccctcTTTCAACCATGGTTGCCCGTACTTCCTAAAAAAAGGGCGATCCGCCGGCTGGGCCAACCCAGACACCAATGCTGATCACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCaaa                                                                             <  1:672404/146‑1 (MQ=255)
                                                                   cactttttacAACATGGTTTGCCTTCCTTCGTCAAAAAAGGGCGATCCTGCATGCGGGCCATCCCAGCCAGCAATGGTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCaaaa                                                                            <  2:779042/142‑1 (MQ=255)
                                                                       gCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACg                                                                         >  1:236731/1‑150 (MQ=255)
                                                                             aaCATGGTTTCCCTTACTTCTTCAAAATAGGGCGATACTGCTTGCTGGGAATCCCAGTCAGCAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTc                                                                   <  2:613318/150‑1 (MQ=255)
                                                                             aaCATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCg                                                                  >  2:464545/1‑151 (MQ=255)
                                                                             aaCATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCg                                                                  >  2:490403/1‑151 (MQ=255)
                                                                              accATGGTTTGCCTTCCTTCGTAAAAAAAGGGCATTCCTGCATGCGGGCAATCCCATCCAGCAATGTTAATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCg                                                                  <  1:762289/149‑1 (MQ=255)
                                                                                 aTGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGc                                                                >  1:341314/1‑150 (MQ=255)
                                                                                   gtcttcccttccTTCGAAAGAAAAAGGCGATCCTGCATGGTGGCCACACCAGTCAGCAATGTGAAGGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGAt                                                                       <  2:207827/130‑1 (MQ=255)
                                                                                      tggCCGTCCTTCGTAAAAAAAGGGCGATCCTGCATGCTGGCCATCCCAGCCAGAAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGc                                                                <  2:1034276/143‑1 (MQ=255)
                                                                                        ccgggaCTGCGTGAAAATAAGGCGATCCTGCAGGCGCGCCATACAAGAAAGAATCGCGGATCACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCGTCGGCACGATg                                                          <  1:558885/144‑1 (MQ=255)
                                                                                               ttCGTAAAAAAAGGGCGATCCTCCATGCGGGGCATCCCAGTCAGCAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGc                                                  <  1:711187/150‑1 (MQ=255)
                                                                                                 cgtaaaaaTAGGGCTATCCTGCATGCTGGGCAACCCAGTCAGAAAGGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGcc                                                 <  2:122603/143‑1 (MQ=255)
                                                                                                     aaaaaaaGGCGATACTGCATGCGGGCCATCCCAGTCAGCAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGtt                                             <  2:738479/144‑1 (MQ=255)
                                                                                                          aaGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTc                                      >  1:616580/1‑151 (MQ=255)
                                                                                                                aTACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAAt                                >  1:461095/1‑151 (MQ=255)
                                                                                                                        tGCTGGGCATACCAGTCAGCAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCa                         <  2:1060490/150‑1 (MQ=255)
                                                                                                                         gCTGGGCATCCCATACAACAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCttt                                   <  1:847945/139‑1 (MQ=255)
                                                                                                                           tGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCt                     >  2:1118767/1‑151 (MQ=255)
                                                                                                                                    ccAGTCAGCATTGCTAATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCAc                  <  1:1038490/146‑1 (MQ=255)
                                                                                                                                      agacacccctcgtgacgaCCCCCGCCAGATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGAc              <  1:1100100/136‑1 (MQ=255)
                                                                                                                                             aatcttgataaCCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACGATCGtt       <  2:1295704/143‑1 (MQ=255)
                                                                                                                                             aaTGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGt    >  2:289139/1‑151 (MQ=255)
                                                                                                                                               tGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCg  >  1:760418/1‑151 (MQ=255)
                                                                                                                                                 gtaatgaCCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCg  <  1:722681/146‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                             
AGATATGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCG  >  NC_000913/2348995‑2349286

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: