Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,506,0380AC41.6% 152.8 / 11.7 77F15V (TTT→GTT) ydcOBenE family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base A (40/5);  new base C (0/32);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 9.45e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATCTCGC  >  NC_000913/1505889‑1506185
                                                                                                                                                      |                                                                                                                                                    
ccATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGaaa                                                                                                                                                      <  2:951118/150‑1 (MQ=255)
   tAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑cccc                                                                                                                                                 >  1:1100968/1‑151 (MQ=255)
     aTGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGc                                                                                                                                               >  2:70938/1‑151 (MQ=255)
      tGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGcc                                                                                                                                              >  1:879880/1‑151 (MQ=255)
      tGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGcc                                                                                                                                              >  1:1079413/1‑151 (MQ=255)
       gtaagagTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGcc                                                                                                                                              >  2:863286/4‑150 (MQ=255)
       gTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGcc                                                                                                                                              <  2:243928/150‑1 (MQ=255)
              tCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCg                                                                                                                                       >  1:1085770/1‑150 (MQ=255)
                      cggccGCCCTTTGCCAGCCCCAGCCCCGTCACCCAGCCAGAGATTTGTGAAGGGGGGGCTCCGGCAACAACCGTTGCTTGCCAGATTATT‑CCGCTTGACCTGGGTTACCAAATTAATCCCGCAAAAAC‑CCCCGCCAGTAGCGTGGGTgga                                                                                                                               <  1:1054256/146‑1 (MQ=255)
                         aCGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTggagga                                                                                                                            >  2:280600/1‑150 (MQ=255)
                           gcccttggccccccccGGCGCCGTCACCCGGCCAGAGTTTGGTGGGGGGGTGCCTCGGGCGACAATCGTTGCTTGCCAGATTATGCCCGC‑TGAACGGGGGTAGCAAATAAATCCCGCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGaa                                                                                                                           <  2:706726/137‑1 (MQ=255)
                              cATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGa                                                                                                                       >  1:47029/1‑150 (MQ=255)
                                ttcccggCCCCAGCGCCGTCACCCAGCCAGAGATTTGTGCAGTGGTGGCCCCGGCAACAACCGCTGTTTGCCAGATTATTGCCGC‑TGAACGGGGGTAGCCAATTAATCCCGCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGaa                                                                                                                      <  1:435344/143‑1 (MQ=255)
                                   ccacccccgccccggTCACCCAGCGAGATTTTTGTGAGGGGGGGGCTCCGGGAAAAACCGCTGTTTCCAAGATTTTTCCGGC‑GGACCGGGGGTAGCAAATTAATACCCCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGaa                                                                                                                      <  2:1194497/132‑1 (MQ=255)
                                   ccAGCCCCGGCGCCCTCACCCAGCCAGATTTTTGTGCAGTGGGGGCTCCGGCGAAAACCGCTGTTTGCCAGATTATTGCCGT‑TGAACGGGGGTAGCAAATAAATCCCGCAAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACa                                                                                                                  <  2:1162878/150‑1 (MQ=255)
                                    cAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAAc                                                                                                                   >  1:700506/1‑148 (MQ=255)
                                    cAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAAc                                                                                                                   >  2:11148/1‑148 (MQ=255)
                                        cccAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGc                                                                                                             >  2:842001/1‑150 (MQ=255)
                                         ccggcgcCGTCACCCAGCAAGAGTTTTGGGCGGTGGTGGCTCCGGCGACATTCTTTGTTGGCCAGATTATTGCCGC‑TAAACGGGCGTACCAAATTAATCCCCCAAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCAt                                                                                                           <  1:843919/148‑1 (MQ=255)
                                             ccccGTCACCCAGCCAGAGTTTGGTGCAGGGGTGGCTCCGGCGAAAATCGCTGTTTCCAAAATTATGGCCGC‑GGAACTGGCGTAGCAAATAAACCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCAt                                                                                                           <  1:331259/145‑1 (MQ=255)
                                               ccttCACCCAGCCAGAGTTTTGTGCTGTGGGGGCTCGGGCAACATTCGCTGCTTGCCAGTTTTTTGCCGC‑TGAACGGGGGTAGCAAATAAATCCCGCAAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAg                                                                                                         <  2:627848/144‑1 (MQ=255)
                                               ccGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCAc                                                                                                     >  1:607803/1‑151 (MQ=255)
                                                  tCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTc                                                                                                  >  1:400363/1‑151 (MQ=255)
                                                    acccgcccagAATTTTGTGGAGTGGTGGCTCGGGCAAAAATCGCTTCTTGCAAGATTATTGCCGC‑TGAATGGGCGTACCAAATTAACCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTc                                                                                                  <  1:71443/143‑1 (MQ=255)
                                                    aTCCAGCCAGAGTTTTGTGCAGTGGGGGCTCCGGCAACAACCGCTGCTTGCCAGATTATTCCCGC‑GGAACGGGGGTAGCAAATAAATCCCGCAAAAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCt                                                                                                   <  1:366175/148‑1 (MQ=255)
                                                      ccggcaagggttttGTGGGGTGGTGGCTCGGGCAAAAATCGCTGCTTCCAATATTATTGCCGT‑GGAACGGGCGTAGCAAATTAATCCCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCg                                                                                                 <  2:950810/137‑1 (MQ=255)
                                                      ccagccagAGATTTGTGCAGGGGTGGCCCCGGCGACAATCGCTGTTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCAAATTAATCCCGCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGt                                                                                                <  1:202707/149‑1 (MQ=255)
                                                            agagTTTTGTGCAGTGGTGGCCCCGGCGAAAATCGCTGCTTGCCAAATTATTGCCGC‑GGAACGGGCGTACCCAATTAATCCCGCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtat                                                                                         <  1:174744/150‑1 (MQ=255)
                                                               gaTTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCg                                                                                     >  1:1045746/1‑151 (MQ=255)
                                                               gaTTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCg                                                                                     >  1:495176/1‑151 (MQ=255)
                                                                 tttGTGCAGGGGTGGCTCCGGCGAAAATCGCTGCTTGCAAGATTATTGCCGC‑TGAACTGGGGTAGCAAATTAATCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAg                                                                                       <  1:1037114/147‑1 (MQ=255)
                                                                    gtgCGGTGGGGGCTCCGGCGACAATCGCTGTTTGCCAGATTATTGCCGC‑GGAACTGGCGTAGCAAATAAATCCCGCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTAc                                                                                  <  1:796490/149‑1 (MQ=255)
                                                                      gcggtggtgGCCCCGGCGACAATCGCTGCTTGCCAGATTATGGCCGC‑TGACCGGGCGTAGCAAATAAATCCCGCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAgg                                                                               <  1:21441/147‑1 (MQ=255)
                                                                         gaggtggCCCCGGCAAAAACCGCGGCTCCCAAATTTATCGCCGA‑TGACTGGGCGTAGCAAATAAACCCCGCAAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTAc                                                                                  <  2:252858/142‑1 (MQ=255)
                                                                          tggtggCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCg                                                                          >  1:676532/1‑151 (MQ=255)
                                                                                 tCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAtt                                                                   >  1:695800/1‑151 (MQ=255)
                                                                                 tCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAtt                                                                   >  1:718521/1‑151 (MQ=255)
                                                                                    ggCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGt                                                                >  1:778302/1‑151 (MQ=255)
                                                                                    ggCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGt                                                                >  1:785840/1‑151 (MQ=255)
                                                                                    ggCAACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACGGGCGTAGCAAATTAATCCCGCAAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAt                                                                    <  1:1063661/147‑1 (MQ=255)
                                                                                     ggaaaaactggcgctTGCCATATTTTTGCCGG‑TAACCGGGGGGAGCAATTTAAACGCGCCAGACC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAttt                                                                  <  2:488826/136‑1 (MQ=255)
                                                                                      cGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTcc                                                              >  1:595212/1‑151 (MQ=255)
                                                                                        aaaaTCGCTGTTGCCCATATTATTGCCGC‑TGACCGGGCGTAGCAATTAAATCCCCCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTcc                                                              <  1:460225/147‑1 (MQ=255)
                                                                                        aCAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGt                                                            >  2:1005757/1‑151 (MQ=255)
                                                                                           aTCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTtat                                                         >  2:843703/1‑151 (MQ=255)
                                                                                              gctctttGCCAGATTTTTGCCGC‑GAAAGGGGCGTAGCAAATTAATCCCCCAAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATaa                                                       <  2:329462/145‑1 (MQ=255)
                                                                                                   ttGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcac                                                  >  1:302278/1‑150 (MQ=255)
                                                                                                   ttGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcac                                                  >  2:798982/1‑150 (MQ=255)
                                                                                                    tGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcaca                                                 >  1:67004/1‑150 (MQ=255)
                                                                                                       cAGTTTATTCCCGC‑TAACCGGGCGTAGCAAATAAATCCCGCAAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcaca                                                 <  1:773861/147‑1 (MQ=255)
                                                                                                        aGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGa                                            <  1:1191429/151‑1 (MQ=255)
                                                                                                           ttattGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATa                                          >  2:1074317/1‑150 (MQ=255)
                                                                                                           ttattCCCGC‑TGACCGGGCGGAACAAATAAATACCGCAAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGAt                                           <  1:369884/149‑1 (MQ=255)
                                                                                                            tattGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATaa                                         >  1:783115/1‑150 (MQ=255)
                                                                                                             attGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAAc                                       >  2:281737/1‑151 (MQ=255)
                                                                                                                 ccGC‑TGACCGGGCGTAGCAAATTAATCCCGCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTcc                                   <  1:1050542/151‑1 (MQ=255)
                                                                                                                   gC‑TGAACGGGGGTAGCAAATTAATCCCGCAAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCg                                  <  2:378414/150‑1 (MQ=255)
                                                                                                                    c‑TGAACGGGCGTAGCAAATAAATCCCGCAAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACtt                                     <  1:698924/146‑1 (MQ=255)
                                                                                                                       gAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtttt                              >  1:1168758/1‑151 (MQ=255)
                                                                                                                        accgggggtaccaaaTAAATCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtttt                              <  1:717342/137‑1 (MQ=255)
                                                                                                                        aacggggttAGCAAATAAACCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCg                                  <  2:123894/138‑1 (MQ=255)
                                                                                                                        aaCTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTg                             >  1:337426/1‑151 (MQ=255)
                                                                                                                        aaCTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTg                             >  1:667541/1‑151 (MQ=255)
                                                                                                                         aCTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGc                            >  1:906714/1‑151 (MQ=255)
                                                                                                                          cTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGcc                           >  2:505755/1‑151 (MQ=255)
                                                                                                                            gggttAGCAAATAAACCCCACCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGc                            <  2:667287/144‑1 (MQ=255)
                                                                                                                                 aGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCggagg                      >  2:1112873/1‑149 (MQ=255)
                                                                                                                                   caaaTTAATCCCGCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGt                   <  1:196195/148‑1 (MQ=255)
                                                                                                                                     aaTTAATCCCGCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCg                 <  1:317241/150‑1 (MQ=255)
                                                                                                                                     aaTTAATACCGCCAGAAA‑CCCCGCCAGAAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCg                 >  2:672622/1‑150 (MQ=255)
                                                                                                                                     aaTTAATACCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGt                   <  1:1032526/148‑1 (MQ=255)
                                                                                                                                         aaacccGCAAAAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCAt              <  2:450061/145‑1 (MQ=255)
                                                                                                                                          aacccccaaGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCAt              <  2:929617/140‑1 (MQ=255)
                                                                                                                                              ccccAGAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGaaa         <  2:943617/148‑1 (MQ=255)
                                                                                                                                              cGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAAt       >  2:842679/1‑151 (MQ=255)
                                                                                                                                               gCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAAtc      >  1:931686/1‑151 (MQ=255)
                                                                                                                                                 cAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAAtctc    >  1:367154/1‑151 (MQ=255)
                                                                                                                                                  aGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATctcg   >  1:542970/1‑151 (MQ=255)
                                                                                                                                                    aaa‑CCCC‑CCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATctcgc  <  1:810414/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGC‑TGAACTGGCGTAGCCAATTAATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAAATCTCGC  >  NC_000913/1505889‑1506185

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