Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,261 G→T 5.8% 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 5.8% 356.6 / 9.2 120intergenic (+31/‑186)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (62/51);  new base T (4/3);  total (66/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

Only 100 of 121 total aligned reads displayed.

CATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAA  >  NC_000913/4060119‑4060410
                                                                                                                                               |                                                                                                                                                       
cATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                                  >  2:1014727/1‑150 (MQ=255)
cATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                                  >  1:1376206/1‑150 (MQ=255)
 aTCGATGATGACGAACTGGGAGCAGTGACTCCGACCTCTATCCGTATGCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCa                                                                                                                                                 <  1:899135/150‑1 (MQ=255)
     atgatgACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑tt                                                                                                                                           >  2:1187764/1‑150 (MQ=255)
      tgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑ttt                                                                                                                                          >  1:129765/1‑150 (MQ=255)
         tgaCGAACTGGTAGAAGTGCCTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑tgtttt                                                                                                                                       >  1:1097724/1‑150 (MQ=255)
          taCGAACTGGGAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑ttttttt                                                                                                                                      <  1:1101096/149‑1 (MQ=255)
          gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑ttttttt                                                                                                                                      >  2:1305563/1‑150 (MQ=255)
           ccGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑tttttttt                                                                                                                                     <  2:1788984/149‑1 (MQ=255)
           aCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAAACCTGCCATCGCTGGCAGG‑‑tttttta                                                                                                                                      <  2:992807/150‑2 (MQ=255)
           aCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑tttttttt                                                                                                                                     <  1:711949/150‑1 (MQ=255)
           aCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑tttttttt                                                                                                                                     >  1:1568950/1‑150 (MQ=255)
            cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑‑ttttttatg                                                                                                                                    <  1:628329/150‑4 (MQ=255)
            cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑ttttttttt                                                                                                                                    >  1:182504/1‑150 (MQ=255)
              atcTGTTAGAAGTGAATCCGACCTCTATCCGTATTCGTAATCGTCACCTGACGTAAAACGATCGTCTCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGc                                                                                                                                  <  1:434587/148‑1 (MQ=255)
               aCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGcc                                                                                                                                 <  1:1014727/150‑1 (MQ=255)
                 tGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGCCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGCCAGG‑‑TTTTTTTTAGCCTg                                                                                                                               >  2:1192061/1‑150 (MQ=255)
                   gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTc                                                                                                                             >  2:973834/1‑150 (MQ=255)
                    tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCg                                                                                                                            >  2:532494/1‑150 (MQ=255)
                    tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCg                                                                                                                            <  2:807575/150‑1 (MQ=255)
                     aGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGACTGTcct                                                                                                                           >  2:1769876/1‑148 (MQ=255)
                        aGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAg                                                                                                                        >  1:1839981/1‑150 (MQ=255)
                             cTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAATCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGcac                                                                                                                   >  2:1657749/1‑150 (MQ=255)
                             cTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTTGCCTGTCGTGAgcgc                                                                                                                     >  2:874756/1‑150 (MQ=255)
                              tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGcaca                                                                                                                  >  2:1346449/1‑150 (MQ=255)
                              tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGcaca                                                                                                                  <  2:254744/150‑1 (MQ=255)
                               ccGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCGTGTCGTGAGCGCacac                                                                                                                 >  2:1819327/1‑150 (MQ=255)
                               ccGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCacac                                                                                                                 >  1:1729164/1‑150 (MQ=255)
                                cGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                                                                                <  2:1020380/150‑1 (MQ=255)
                                   ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAAt                                                                                                             >  1:756012/1‑150 (MQ=255)
                                     tctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGc                                                                                                           >  2:1041282/1‑150 (MQ=255)
                                      ctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACCCCAATGc                                                                                                           >  2:1110107/1‑150 (MQ=255)
                                       tATCCGTATTCGTAAACGTCACCTGACGGAAAACGATGGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                         <  2:602162/150‑1 (MQ=255)
                                        aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                                                                        <  1:1177809/150‑1 (MQ=255)
                                        aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                                                                        <  1:1339726/150‑1 (MQ=255)
                                                 cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccat                                                                                               >  2:1486708/1‑148 (MQ=255)
                                                   tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                                                                             <  2:523918/150‑1 (MQ=255)
                                                   tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                                                                             <  2:1564043/150‑1 (MQ=255)
                                                      aCGTCACCAGACGGAAATAGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                                                                          <  2:1150119/150‑1 (MQ=255)
                                                       cGTCACCTGACGGAAATCGATCGTCGCAGCGCCAACGGCGGACCGAAAGACGATTATTTTCGTTTTTTAGTCATAAAAAAACTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACGAATGCGAGTATCTCCCTTATAAg                                                                                         <  2:1642772/150‑1 (MQ=255)
                                                       cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCGCACCAATGCGTGTATCGCCCTTATAAg                                                                                         >  2:118093/1‑150 (MQ=255)
                                                       cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGt                                                                                        >  2:532092/1‑150 (MQ=255)
                                                         tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGtg                                                                                       >  1:1817213/1‑149 (MQ=255)
                                                          cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCt                                                                                      >  1:1512009/1‑150 (MQ=255)
                                                           aCCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGCGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtg                                                                                     <  2:1113018/150‑1 (MQ=255)
                                                              tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCTTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                                  <  2:1138581/150‑1 (MQ=255)
                                                              tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGCGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                                  >  1:1227174/1‑150 (MQ=255)
                                                              tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                                  <  2:1616075/150‑1 (MQ=255)
                                                              tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                                  <  2:1220492/150‑1 (MQ=255)
                                                              tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                                  >  2:1492672/1‑150 (MQ=255)
                                                                 cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGTTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCTGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACTCCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                                                                               >  1:302456/1‑150 (MQ=255)
                                                                 cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                                                                               <  1:26722/150‑1 (MQ=255)
                                                                  ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑GTTTTTTTTGCCTGTCGTGAGCGCACACCAATTCGAGTTTGTGCCGTATAAGTAGTTGGTTAc                                                                              >  2:1128142/1‑150 (MQ=255)
                                                                  ggAAAACGATCGTCGCCACGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                                                                              >  1:438046/1‑150 (MQ=255)
                                                                   gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                                                                             >  1:627314/1‑150 (MQ=255)
                                                                   gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTAAGTCATAAAAAACCTGCCAGCGATGGCAGC‑‑TGATTTTTAGCCTGTAGTGAGCGCACACCAAGGCGACTATCTCCCTTATAAGTCTGTGGGTACg                                                                             >  2:256732/1‑150 (MQ=255)
                                                                      aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGACAGAGATGGCAGG‑‑TTTTTTTTTGCATGTCGTGAGCGCACACCAATGCGAGTATATCCCTTATAAGTCTGTGGTTACGTCa                                                                          >  2:1054158/1‑150 (MQ=255)
                                                                            cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                                                                    >  2:608901/1‑150 (MQ=255)
                                                                             gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagt                                                                   >  2:1688537/1‑149 (MQ=255)
                                                                             gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                                                                   >  1:1048999/1‑150 (MQ=255)
                                                                               cgccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                                                                 >  1:460051/1‑150 (MQ=255)
                                                                                 ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt                                                               >  1:1172555/1‑150 (MQ=255)
                                                                                 ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCATTATAAGTCTGTGGTTACGTCAATAGAGAGCtt                                                               >  1:1353560/1‑150 (MQ=255)
                                                                                  cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg                                                             <  1:532092/150‑1 (MQ=255)
                                                                                   gcgcCACCCGCTCACCGAAAGTCGATTAATTTCGTTTTTTAGTCATAAAAAACCTTCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg                                                             <  2:1211605/150‑1 (MQ=255)
                                                                                   gcgcCAACCGCGCACCTAAAGACTTTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCTTGGGAAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg                                                             <  2:1376206/150‑1 (MQ=255)
                                                                                   gcgcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTTTCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg                                                             >  1:1445970/1‑150 (MQ=255)
                                                                                   gcgcCAACCGCGCACCGAAAGACCATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGAAGGCAGG‑‑TTTTTTTTTGCCTGTAGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGATACGTCAATAGAGAGCTTTg                                                             <  2:1217707/150‑1 (MQ=255)
                                                                                    cgcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCGGCCAGCGTTGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt                                                            <  2:451599/150‑1 (MQ=255)
                                                                                     gcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGta                                                           >  1:1087204/1‑149 (MQ=255)
                                                                                      ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca                                                          <  2:688119/150‑1 (MQ=255)
                                                                                       cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTTGTTACGTCAATAGAGAGCTTTGTcac                                                         <  2:129765/150‑1 (MQ=255)
                                                                                         aCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa                                                       >  1:967363/1‑150 (MQ=255)
                                                                                          ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt                                                      <  1:2875/150‑1 (MQ=255)
                                                                                            gcgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTa                                                    >  1:1021633/1‑150 (MQ=255)
                                                                                              gcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATc                                                  <  1:19792/150‑1 (MQ=255)
                                                                                              gaactgAATGACGATTATTTTCTTTTTTTAGTCTTAAAAAACCTGCCAGCGATGGAAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCGGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATc                                                  <  1:1269319/145‑1 (MQ=255)
                                                                                                 ccGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCTAGCGATGGCAGG‑‑TTTTTTTTTGCCTTTCGTGAGCGCAAACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGCGCTTAGTCACAATTATCTGc                                               >  2:289892/1‑150 (MQ=255)
                                                                                                  cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCa                                              >  1:1855275/1‑150 (MQ=255)
                                                                                                   gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCaa                                             <  1:1305563/150‑1 (MQ=255)
                                                                                                    aaaGACGATAAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTATGTGGTTACGTCAATAGATAGCTTTGTCACAATTATCTGCaaa                                            >  1:1186850/1‑150 (MQ=255)
                                                                                                              aaTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGcc                                  >  1:1329533/1‑150 (MQ=255)
                                                                                                              aaTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCGCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGcc                                  <  2:104191/150‑1 (MQ=255)
                                                                                                               aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCg                                 <  2:131165/150‑1 (MQ=255)
                                                                                                                  tCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTa                              >  2:718025/1‑150 (MQ=255)
                                                                                                                              cATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCAAAATTATCTGCAAAGTCATACGCCGTTAATTGctttcttt                  >  2:1214506/1‑150 (MQ=255)
                                                                                                                               aaaaaaaaaCCTGCCATCGCTGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctttcttt                  >  2:591382/4‑150 (MQ=255)
                                                                                                                                   aaaaaCCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTg              <  1:1178418/150‑1 (MQ=255)
                                                                                                                                   aaaaaCCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTg              <  2:572512/150‑1 (MQ=255)
                                                                                                                                   aaaaaCCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTg              <  2:1473544/150‑1 (MQ=255)
                                                                                                                                    aaaaCCTGCCATCGCTGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATTAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTgg             >  1:220246/1‑150 (MQ=255)
                                                                                                                                    aaaaCCTGCCATCGCTGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTgg             >  1:1299123/1‑150 (MQ=255)
                                                                                                                                    aaaaCCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTgg             <  2:1340852/150‑1 (MQ=255)
                                                                                                                                    aaaaCCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTgg             >  2:1502654/1‑150 (MQ=255)
                                                                                                                                    aaaaCCTGCCAGCGAGGTCAGG‑‑TTTTTTTTTGCCTGTCGTGAGAGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTgg             <  2:903945/150‑1 (MQ=255)
                                                                                                                                      aaCCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcg           >  2:387518/1‑150 (MQ=255)
                                                                                                                                          tggcAGCGATTGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATCAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACTCCGTTAATTGCTTTCTTTTTTGgcgtaag       <  1:305696/147‑1 (MQ=255)
                                                                                                                                          tGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaag       >  1:531094/1‑150 (MQ=255)
                                                                                                                                              aGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaagcgta   <  2:1568950/150‑1 (MQ=255)
                                                                                                                                              aGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaagcgta   <  1:1657749/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                       
CATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAA  >  NC_000913/4060119‑4060410

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: