Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,261 G→T 6.0% intergenic (+31/‑186) typA → / → yihL GTP‑binding protein/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,060,2610GT 6.0% 1826.4 / 64.0 499intergenic (+31/‑186)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (228/240);  new base T (18/13);  total (246/253)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 504 total aligned reads displayed.

TGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAA  >  NC_000913/4060112‑4060410
                                                                                                                                                     |                                                                                                                                                      
    gTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGAt                                                                                                                                                    <  2:7173403/150‑1 (MQ=255)
    gTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGAt                                                                                                                                                    >  2:5337676/1‑150 (MQ=255)
     ttCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATg                                                                                                                                                   >  1:1324705/1‑150 (MQ=255)
      tCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                                                                                  <  2:86797/150‑1 (MQ=255)
       cATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                                 >  1:7148999/1‑150 (MQ=255)
       cATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                                 >  2:1171009/1‑150 (MQ=255)
        atagatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTTACGGAAAAAGATCGACGCCGCGCAAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTTACAGCGATGGCa                                                                                                                                                >  1:4673182/4‑150 (MQ=255)
        atagatgatGAAGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAAAGTCAACTGACGGAAAACGATCGTCGCAGCGCCAACCGAGCACCGAAAGAAGATTAATTTCGTTTTTTAGTCATAAAAAACATGCCAGAGATggaa                                                                                                                                                >  1:6402438/4‑148 (MQ=255)
        aTCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCa                                                                                                                                                <  1:599717/150‑1 (MQ=255)
         tCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAg                                                                                                                                               <  1:613417/150‑1 (MQ=255)
         tCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAg                                                                                                                                               >  2:4511179/1‑150 (MQ=255)
         tCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAg                                                                                                                                               <  2:4637651/150‑1 (MQ=255)
         tCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAg                                                                                                                                               >  2:5198403/1‑150 (MQ=255)
           gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑t                                                                                                                                            <  2:3278727/150‑1 (MQ=255)
           gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑t                                                                                                                                            >  1:5915113/1‑150 (MQ=255)
             tgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttt                                                                                                                                          >  1:5584049/1‑150 (MQ=255)
                tgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tttttt                                                                                                                                       >  2:1043462/1‑150 (MQ=255)
                 gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttttttt                                                                                                                                      <  1:2771862/150‑1 (MQ=255)
                  aCGAACTGGTAGAAGTGACTCCGACCTCTCTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCACCCGCGCACCGAAAGACGATTAATATCGTTTTTTCGTCATAAAAAACCTGCCAGCGATGGCAGG‑tatttttt                                                                                                                                     <  2:4221320/150‑1 (MQ=255)
                  aCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttttttat                                                                                                                                     <  1:5791659/150‑3 (MQ=255)
                   cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttttttttt                                                                                                                                    <  1:3903548/150‑1 (MQ=255)
                   cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttttttttt                                                                                                                                    <  1:4904881/150‑1 (MQ=255)
                    gAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTg                                                                                                                                   <  1:4773846/150‑1 (MQ=255)
                       cTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCt                                                                                                                                <  2:2676429/150‑1 (MQ=255)
                          gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGc                                                                                                                                  <  2:3434870/145‑1 (MQ=255)
                          gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTc                                                                                                                             >  2:5534500/1‑150 (MQ=255)
                            aGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGt                                                                                                                           >  1:707059/1‑150 (MQ=255)
                            aGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGt                                                                                                                           <  2:1474933/150‑1 (MQ=255)
                                gTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAgc                                                                                                                       >  1:6304621/1‑150 (MQ=255)
                                     tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCGTGTCGTGTGCGcaca                                                                                                                  >  1:7234838/1‑150 (MQ=255)
                                      ccGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCacac                                                                                                                 <  2:5890695/150‑1 (MQ=255)
                                       cGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                                                                                <  1:686003/150‑1 (MQ=255)
                                        gACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTTTCGTGAGAGCACACCa                                                                                                               >  2:5652959/1‑150 (MQ=255)
                                        gACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCa                                                                                                               >  2:2134647/1‑150 (MQ=255)
                                           ctctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATg                                                                                                            >  1:5290017/1‑150 (MQ=255)
                                           ctctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATg                                                                                                            >  1:4016279/1‑150 (MQ=255)
                                           ctctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGAGTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATg                                                                                                            >  2:1972810/1‑150 (MQ=255)
                                            tctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGc                                                                                                           <  1:504121/150‑1 (MQ=255)
                                            tctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGc                                                                                                           >  2:6779665/1‑150 (MQ=255)
                                             ctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                                                                          <  1:1008317/150‑1 (MQ=255)
                                              tATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                         <  1:1770537/150‑1 (MQ=255)
                                              tATCCGCATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                         <  2:5248646/150‑1 (MQ=255)
                                                    tATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtct                                                                                                   >  2:2749309/1‑150 (MQ=255)
                                                    tATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtct                                                                                                   >  1:1559360/1‑150 (MQ=255)
                                                     aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                                                                  <  2:6600488/150‑1 (MQ=255)
                                                      ttCTTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCTGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                                                                                 <  1:810899/150‑1 (MQ=255)
                                                       tCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                                                                                >  2:513279/1‑150 (MQ=255)
                                                       tCGTAAACGTCACCGGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGACAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                                                                                <  2:3682896/150‑1 (MQ=255)
                                                        cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                                                                               <  2:4365810/150‑1 (MQ=255)
                                                          tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                                                                             >  1:3365910/1‑150 (MQ=255)
                                                           aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                                                                            <  1:694011/150‑1 (MQ=255)
                                                                   ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgt                                                                                    <  1:2925819/150‑1 (MQ=255)
                                                                    cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                                                                   <  1:456800/150‑1 (MQ=255)
                                                                    cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                                                                   <  1:3768378/150‑1 (MQ=255)
                                                                        cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                                                                               <  2:6222935/150‑1 (MQ=255)
                                                                         ggACAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAAGTTCGTTTGTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGGCCTGTCGTGAGCGCACACCCATGCGAGTCTCTCCCTTATAAGTCTGTGGTTAc                                                                              >  2:6984203/1‑150 (MQ=255)
                                                                         ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAAGGTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTCTCTCCCTTATAAGTCTGTGGTTAc                                                                              >  2:6678973/1‑150 (MQ=255)
                                                                           aaaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGt                                                                            <  1:2070846/150‑1 (MQ=255)
                                                                                gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                                                                       <  2:3579824/150‑1 (MQ=255)
                                                                                gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                                                                       >  2:2364197/1‑150 (MQ=255)
                                                                                  tcgtcgCCGCGCCAACCGCGCACCGCAATAAGTTTTATTTCGTTTTTTAGTCATAAAAAACCTGCCAGGGTTGGCATG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                                                                     <  1:6759443/150‑1 (MQ=255)
                                                                                    gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                                                                   <  2:1771163/150‑1 (MQ=255)
                                                                                     tcgCCGCGCCAACCGCGCACCGAACGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                                                                  <  1:1918579/150‑1 (MQ=255)
                                                                                     tcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCACCCTTATAAGTCTGTGGTTAAGTCAAAagagag                                                                  >  1:6379521/1‑150 (MQ=255)
                                                                                        ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTTTAGTGAGCGCACACCAATGCGAGTATCTCCCTTAAAAGTATGAGGTTACGTCAATAGAGAGCtt                                                               >  2:6433125/1‑150 (MQ=255)
                                                                                         cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt                                                              <  1:5655253/150‑1 (MQ=255)
                                                                                         cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ATTATTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt                                                              <  1:5735715/150‑1 (MQ=255)
                                                                                           cgcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                                                                              <  1:5206856/132‑1 (MQ=255)
                                                                                            gcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc                                                           >  1:4120435/1‑150 (MQ=255)
                                                                                            gcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc                                                           <  1:3151762/150‑1 (MQ=255)
                                                                                             ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc                                                           <  1:2617555/150‑1 (MQ=255)
                                                                                             ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca                                                          <  2:1087856/150‑1 (MQ=255)
                                                                                                aCCGCGCACCGAACGACTATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCCGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa                                                       <  1:145553/150‑1 (MQ=255)
                                                                                                       aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTg                                                <  1:4393749/150‑1 (MQ=255)
                                                                                                          gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCaa                                             >  1:3480908/1‑150 (MQ=255)
                                                                                                           aaaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCaaa                                            <  2:4982454/150‑1 (MQ=255)
                                                                                                            aaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAg                                           <  1:6493252/150‑1 (MQ=255)
                                                                                                               aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTc                                         <  2:410607/150‑1 (MQ=255)
                                                                                                               aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCa                                        <  1:124205/150‑1 (MQ=255)
                                                                                                                cGATTAATTTCTTTTTTTAGTCATAAAAAACCTGCCATCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTTCGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCAt                                       <  2:5732551/150‑1 (MQ=255)
                                                                                                                 gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCaca                                      <  2:4839221/150‑3 (MQ=255)
                                                                                                                 gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCATATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATa                                      >  1:1657010/1‑150 (MQ=255)
                                                                                                                 gATTAAATTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTTCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCATATAAGTCTGTGGTTACGTCAATAGAGATCTTTGTCACAATAATCTGCAAAGTCATa                                      >  2:1656436/1‑150 (MQ=255)
                                                                                                                  aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTATGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATAc                                     >  1:5383975/1‑150 (MQ=255)
                                                                                                                         tCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTa                              >  2:6106628/1‑150 (MQ=255)
                                                                                                                         tCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTa                              >  2:6928691/1‑150 (MQ=255)
                                                                                                                         tCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGGTa                              >  2:378422/1‑150 (MQ=255)
                                                                                                                         tCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTAGAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTa                              >  2:2509798/1‑150 (MQ=255)
                                                                                                                                ttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctt                       <  1:6678350/150‑1 (MQ=255)
                                                                                                                                        aaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTATTGCCTGTCGTGAGAGCACACCAATTCGAGTATCTCCCTAATAAGTCAGTGGTTACGTCAATAGAGAGCTTTGTCACAATTGTCTGCAAAGTCATACGCCGTGAATTGCTTTCtgtttt               >  2:2986744/1‑150 (MQ=255)
                                                                                                                                          aaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTgg             >  2:3136041/1‑150 (MQ=255)
                                                                                                                                          aaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTgg             >  1:3602102/1‑150 (MQ=255)
                                                                                                                                           aaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgc            >  1:1128038/1‑150 (MQ=255)
                                                                                                                                           aaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgc            >  1:2416157/1‑150 (MQ=255)
                                                                                                                                           aaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgc            >  2:3251376/1‑150 (MQ=255)
                                                                                                                                            aaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcg           >  1:2181797/1‑150 (MQ=255)
                                                                                                                                               cTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaa        <  2:5437874/150‑1 (MQ=255)
                                                                                                                                               cTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaa        <  1:4230005/150‑1 (MQ=255)
                                                                                                                                                    aGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaagcgta   <  1:3244923/150‑1 (MQ=255)
                                                                                                                                                     gCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaagcgtaa  >  2:3987778/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAA  >  NC_000913/4060112‑4060410

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