Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,261 G→T 7.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 7.0% 248.9 / 7.8 108intergenic (+31/‑186)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (58/42);  new base T (2/6);  total (60/48)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.25e-01

Only 100 of 110 total aligned reads displayed.

GGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTA  >  NC_000913/4060113‑4060409
                                                                                                                                                    |                                                                                                                                                    
ggAGTGCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAg                                                                                                                                                      <  1:979487/149‑1 (MQ=255)
   gTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGa                                                                                                                                                   >  2:554382/1‑149 (MQ=255)
       aTCGATGATGACGAACTGGTAGAAGTGACTCCTCCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAAAGATCGTCGCCGCGCCAACCGCGCCCCGAAAGACGATTAATTACGTTTTTTATTCATAAAAAACCAGCCAGCGATGGc                                                                                                                                               >  2:745812/1‑149 (MQ=255)
       aTCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACTGAAAACGATCGTCGCCGCTCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                               <  2:372953/149‑1 (MQ=255)
       aTCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCACCCTCGCACCGAAAGACGATTAATTTCGTTTTTTAGACATAAAAAACCTGACAGCGATGGc                                                                                                                                               >  1:696828/1‑149 (MQ=255)
       aTCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                               >  1:679638/1‑149 (MQ=255)
       aTCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                               >  1:182091/1‑149 (MQ=255)
        tctatgatgACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTAATAAAAAACCTGACAGCGATGGCa                                                                                                                                              >  1:4091/4‑149 (MQ=255)
         cGATGAAGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAg                                                                                                                                             <  2:256982/149‑1 (MQ=255)
          gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAgg                                                                                                                                            <  1:238665/149‑1 (MQ=255)
          gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATGTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAgg                                                                                                                                            >  2:928876/1‑149 (MQ=255)
          gatgatGACGAACTGGTAGAAGTGACTCAGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAgg                                                                                                                                            <  1:454886/149‑1 (MQ=255)
            tgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                                                                          >  1:980042/1‑149 (MQ=255)
            tgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATAAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                                                                          >  2:789847/1‑149 (MQ=255)
              atgaCGAACTGGTAGAAGTGACTCCGACCTCTATCAGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttt                                                                                                                                        <  2:8540/149‑1 (MQ=255)
               tgaCGAACTGGTAGAAGTGAATCCGACCTCTATCCGTATTCGTAAACGTCACCTGAAGGAAAACGATCGTCGGCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttt                                                                                                                                       <  2:881873/149‑1 (MQ=255)
                gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCACCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATTAAAAACCTGCCAGCGATGGCAGGtttttt                                                                                                                                      >  1:666749/1‑149 (MQ=255)
                gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGGCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttt                                                                                                                                      <  2:72038/149‑1 (MQ=255)
                gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGGCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttt                                                                                                                                      <  2:342543/149‑1 (MQ=255)
                  cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttat                                                                                                                                    <  2:378487/149‑3 (MQ=255)
                      cTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGAGATGGCAGGTTTTTTTTTGcc                                                                                                                                >  1:654407/1‑149 (MQ=255)
                         ggaGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTgg                                                                                                                             >  2:343749/3‑148 (MQ=255)
                         gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGCCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGt                                                                                                                             >  1:1006034/1‑149 (MQ=255)
                         gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTgg                                                                                                                             >  2:843869/1‑148 (MQ=255)
                         gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGt                                                                                                                             >  1:1006895/1‑149 (MQ=255)
                          tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTc                                                                                                                            <  1:395703/149‑1 (MQ=255)
                             aaGTGACTCCGACCTCTATCCGTATTCGTATACGTCACCTGCCGGAAAACGATCGTCGCCGCGCCAACCGCGCGCCGAAAGACGATTAATTTCGTTATTTAGTCATAAAAACCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTg                                                                                                                         >  2:813259/1‑149 (MQ=255)
                               gTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTAGCCGCGCCAACCGCGCACCGAAAGACGATTAATGTCGTTTTTTAGTCATAAAAAACCTGCCAGCGAAGGCACGTTATTATTTGCCTGTCGTGAg                                                                                                                       >  1:743379/1‑149 (MQ=255)
                                    tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCAGCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGcac                                                                                                                  >  2:917311/1‑149 (MQ=255)
                                       gACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTTTCGTGAGCGCACAcc                                                                                                               >  1:203800/1‑149 (MQ=255)
                                         ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAgcg                                                                                                                     <  1:517901/141‑1 (MQ=255)
                                         ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAgcg                                                                                                                     >  2:517901/1‑141 (MQ=255)
                                           tctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCGACCGCGCCCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATg                                                                                                           >  2:79163/1‑149 (MQ=255)
                                            ctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGc                                                                                                          >  2:154214/1‑149 (MQ=255)
                                              aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATTGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                        >  2:849480/1‑149 (MQ=255)
                                              aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCCCACCAATGCGa                                                                                                        >  1:110846/1‑149 (MQ=255)
                                              aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                        >  2:264538/1‑149 (MQ=255)
                                              aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                        >  2:241263/1‑149 (MQ=255)
                                                 cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                                                                     >  2:65931/1‑149 (MQ=255)
                                                 cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACAAATGCGAGTa                                                                                                     >  1:237446/1‑149 (MQ=255)
                                                    aTTAGTAAACGTAACCTGACGGAAAACGATGGTCGCCGCGACAACCGCGCACCGAAAGACGTATAATTTCGTTGTTTAGACATAAAAAACCTGCCAGCGAGGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtct                                                                                                  <  1:663320/149‑1 (MQ=255)
                                                     tgagTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                                                                 <  1:459898/146‑1 (MQ=255)
                                                       cGTAAACGTCAGCTGACGGAAAAAGATCGTATCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGGGATGTCAGGTTTTTTTTTGACTGTTGTGAGCGCACACCAATGCGAGTATCTccc                                                                                               <  2:812249/149‑1 (MQ=255)
                                                         tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGCTGGCAGGTTTTGTTTTGCCTGTCGTGAGCGCACACAAATGCGAGTATCTCCCtt                                                                                             >  1:1009663/1‑149 (MQ=255)
                                                              gTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                                                                        >  2:301200/1‑149 (MQ=255)
                                                               tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGt                                                                                       <  1:360752/149‑1 (MQ=255)
                                                                  ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGAGGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACCCCAATGCGAGTATCTCCCGTATATGTcgg                                                                                    >  2:618439/1‑147 (MQ=255)
                                                                   cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgt                                                                                   >  2:1026038/1‑149 (MQ=255)
                                                                    tGACGGAAAACGTTAGTCGCCGTGCCAACCGGGGACCGTAAGTCTATTAATTTCGTTTTTTAGTCATAAAATACCTGCCAGCGTTGGCATGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGAGAGTATCTCCCTTATAAGTCtgtg                                                                                  <  2:517659/149‑1 (MQ=255)
                                                                    tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGTCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                                                                  <  2:756276/149‑1 (MQ=255)
                                                                    tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                                                                  >  1:12150/1‑149 (MQ=255)
                                                                    tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                                                                  <  1:568278/149‑1 (MQ=255)
                                                                     gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                                 <  2:803289/149‑1 (MQ=255)
                                                                       cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGtt                                                                               <  1:660838/149‑1 (MQ=255)
                                                                        ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTTTCGTGAGCGCACACCAATGCTAGTATCTCCCTTATAAGTCTGTGGTTa                                                                              >  2:856053/1‑149 (MQ=255)
                                                                        ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACCCCAAGGCGAGTATCTCCCTTATAAGTCTGTGGTTa                                                                              >  2:730330/1‑149 (MQ=255)
                                                                        ggAAAACGATCGTCGACGCGCCAACCGCGCACCGAAAGACGTTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATTGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                                                                              <  1:677056/149‑1 (MQ=255)
                                                                               taTCGTAGCAGCGCAAAACGAGCACCGAAAGAAGATTAATTTAGTTTTTTAGTAATAAATAACCTGCCAGCGATGGCAGGTTTTTTTTTGACTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                                                                       <  1:986581/148‑1 (MQ=255)
                                                                               gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTGTGTGGTTACGTCAAt                                                                       >  2:224921/1‑149 (MQ=255)
                                                                                 tcgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                                                                     >  2:277734/1‑149 (MQ=255)
                                                                                   gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCTTGAGCGCACACCAATGCGATTATCTCCCTTATAAGTATGTGGTTACGTCAATagtg                                                                   >  2:957474/1‑147 (MQ=255)
                                                                                   gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                                                                   >  2:296829/1‑149 (MQ=255)
                                                                                    tcgCCGCGCCAACCGCGCACCGAAAGACGATTAATGTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                                                                  >  1:375777/1‑149 (MQ=255)
                                                                                      gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                                                                >  2:533492/1‑149 (MQ=255)
                                                                                       ccgcGCCCACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCt                                                               >  1:133360/1‑149 (MQ=255)
                                                                                           gcCAACCGCGCACCGAAAGACGATTAATTTAGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt                                                           <  1:277742/149‑1 (MQ=255)
                                                                                            ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTTAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc                                                          <  2:797931/149‑1 (MQ=255)
                                                                                              aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac                                                        >  1:275860/1‑149 (MQ=255)
                                                                                                ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGCTGGCAGGTTTTTTTTTGCCTGTCGGGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa                                                      <  2:204262/149‑1 (MQ=255)
                                                                                                      aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCt                                                <  1:366790/149‑1 (MQ=255)
                                                                                                      aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCt                                                <  1:104041/149‑1 (MQ=255)
                                                                                                      aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCt                                                >  1:101648/1‑149 (MQ=255)
                                                                                                          aaaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACGAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTAACTGCaa                                            >  2:425762/1‑149 (MQ=255)
                                                                                                             gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGt                                         >  2:686553/1‑149 (MQ=255)
                                                                                                               cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCa                                       >  2:885094/1‑149 (MQ=255)
                                                                                                               cGATTAATTTCGTATTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGATCGCACACCAATGCGAGTATCTCCCTTATAACTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCa                                       <  2:920616/149‑1 (MQ=255)
                                                                                                                 aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTGTTTGCCTGTCGTGCGCGCACACCAATGCGAGGATCTCCCTTATAAGTCTGTGGTTACGTGAATAGAGCGCGTTGTCACAATTATCTTCAAAGTCATa                                     >  2:462191/1‑149 (MQ=255)
                                                                                                                        tCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTGTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTTTCTGCCTTTTAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTGACAATTATGTGGAAAGTCATACGCCGtt                              >  2:492989/1‑149 (MQ=255)
                                                                                                                            tttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAAtt                          <  1:65931/149‑1 (MQ=255)
                                                                                                                               ttattaaTACAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCATATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGct                       <  2:207851/143‑1 (MQ=255)
                                                                                                                                 aGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGcttt                     <  1:798824/149‑1 (MQ=255)
                                                                                                                                  gcaaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACAAATGCGAGTATCTCCCTTTTAAGTCTGTGGTTACGTCAAAAGAGAGCTTTGTCAGAATTATCTGCAAAGTCATACGCCTTTTAGTGctttg                    >  1:1017518/6‑148 (MQ=38)
                                                                                                                                   tCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctttct                   >  2:186623/1‑149 (MQ=255)
                                                                                                                                     aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAAGGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctttcttt                 >  2:473363/1‑149 (MQ=255)
                                                                                                                                       aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGGGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCttatt               >  2:213887/1‑149 (MQ=255)
                                                                                                                                        aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttttt              <  1:533492/149‑1 (MQ=255)
                                                                                                                                        aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttttt              <  2:730129/149‑1 (MQ=255)
                                                                                                                                        aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGACCGCACACCAATTCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGACCTTAGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttttt              >  2:377533/1‑149 (MQ=255)
                                                                                                                                         aacaCCTCCCAGCGATGGCAGGGTTTTTTTTGCCTGTCGCTAGCGCACACCAATGCGAGTATCTCCCTAATAAGTCTGGGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTg             <  2:666749/146‑1 (MQ=255)
                                                                                                                                         aaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTg             <  1:928876/149‑1 (MQ=255)
                                                                                                                                          aaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTgg            >  2:552874/1‑149 (MQ=255)
                                                                                                                                           aaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACATTTATCGGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgc           <  1:745812/149‑1 (MQ=255)
                                                                                                                                           aaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgc           <  2:392829/149‑1 (MQ=255)
                                                                                                                                           aaCCTGCCAGCGATGGCAGGTTTTTTATTGCCTGTCGAGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgc           <  1:849480/149‑1 (MQ=255)
                                                                                                                                                gCCAGCGATGGCAGGTTTTTTTTTGCCTGGCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaag      <  1:354481/149‑1 (MQ=255)
                                                                                                                                                gCCAGCGAGGGCAGGTTTTTTTTTGCATGTCGTGAGCGCACACCAATGCGAGTAGCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaag      <  1:902358/149‑1 (MQ=255)
                                                                                                                                                 ccAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaagc     <  1:1026038/149‑1 (MQ=255)
                                                                                                                                                 ccAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATAATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaagc     <  1:447359/149‑1 (MQ=255)
                                                                                                                                                   aGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaagcgt   <  1:241263/149‑1 (MQ=255)
                                                                                                                                                    gCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaagcgta  <  1:264538/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTA  >  NC_000913/4060113‑4060409

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