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.

CATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAA  >  NC_000913/4060119‑4060410
                                                                                                                                               |                                                                                                                                                       
cATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGc                                                                                                                                                  >  2:1014727/1‑150 (MQ=255)
cATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGc                                                                                                                                                  >  1:1376206/1‑150 (MQ=255)
 aTCGATGATGACGAACTGGGAGCAGTGACTCCGACCTCTATCCGTATGCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCa                                                                                                                                                 <  1:899135/150‑1 (MQ=255)
     atgatgACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑tt                                                                                                                                           >  2:1187764/1‑150 (MQ=255)
      tgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑ttt                                                                                                                                          >  1:129765/1‑150 (MQ=255)
         tgaCGAACTGGTAGAAGTGCCTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑tgtttt                                                                                                                                       >  1:1097724/1‑150 (MQ=255)
          taCGAACTGGGAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑ttttttt                                                                                                                                      <  1:1101096/149‑1 (MQ=255)
          gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑ttttttt                                                                                                                                      >  2:1305563/1‑150 (MQ=255)
           ccGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑tttttttt                                                                                                                                     <  2:1788984/149‑1 (MQ=255)
           aCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAAACCTGCCATCGCTGGCAGG‑‑tttttta                                                                                                                                      <  2:992807/150‑2 (MQ=255)
           aCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑tttttttt                                                                                                                                     <  1:711949/150‑1 (MQ=255)
           aCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑tttttttt                                                                                                                                     >  1:1568950/1‑150 (MQ=255)
            cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCATCGCTGGCAGG‑‑ttttttatg                                                                                                                                    <  1:628329/150‑4 (MQ=255)
            cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑ttttttttt                                                                                                                                    >  1:182504/1‑150 (MQ=255)
              atcTGTTAGAAGTGAATCCGACCTCTATCCGTATTCGTAATCGTCACCTGACGTAAAACGATCGTCTCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGc                                                                                                                                  <  1:434587/148‑1 (MQ=255)
               aCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGcc                                                                                                                                 <  1:1014727/150‑1 (MQ=255)
                 tGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGCCGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGCCAGG‑‑TTTTTTTTAGCCTg                                                                                                                               >  2:1192061/1‑150 (MQ=255)
                   gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTc                                                                                                                             >  2:973834/1‑150 (MQ=255)
                    tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCg                                                                                                                            >  2:532494/1‑150 (MQ=255)
                    tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCg                                                                                                                            <  2:807575/150‑1 (MQ=255)
                     aGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGACTGTcct                                                                                                                           >  2:1769876/1‑148 (MQ=255)
                        aGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAg                                                                                                                        >  1:1839981/1‑150 (MQ=255)
                             cTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAATCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGcac                                                                                                                   >  2:1657749/1‑150 (MQ=255)
                             cTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTTGCCTGTCGTGAgcgc                                                                                                                     >  2:874756/1‑150 (MQ=255)
                              tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGcaca                                                                                                                  >  2:1346449/1‑150 (MQ=255)
                              tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGcaca                                                                                                                  <  2:254744/150‑1 (MQ=255)
                               ccGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCGTGTCGTGAGCGCacac                                                                                                                 >  2:1819327/1‑150 (MQ=255)
                               ccGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCacac                                                                                                                 >  1:1729164/1‑150 (MQ=255)
                                cGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                                                                                <  2:1020380/150‑1 (MQ=255)
                                   ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAAt                                                                                                             >  1:756012/1‑150 (MQ=255)
                                     tctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGc                                                                                                           >  2:1041282/1‑150 (MQ=255)
                                      ctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGT‑TTTTTTTTTGCCTGTCGTGAGCGCACCCCAATGc                                                                                                           >  2:1110107/1‑150 (MQ=255)
                                       tATCCGTATTCGTAAACGTCACCTGACGGAAAACGATGGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                         <  2:602162/150‑1 (MQ=255)
                                        aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                                                                        <  1:1177809/150‑1 (MQ=255)
                                        aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                                                                        <  1:1339726/150‑1 (MQ=255)
                                                 cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccat                                                                                               >  2:1486708/1‑148 (MQ=255)
                                                   tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                                                                             <  2:523918/150‑1 (MQ=255)
                                                   tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                                                                             <  2:1564043/150‑1 (MQ=255)
                                                      aCGTCACCAGACGGAAATAGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                                                                          <  2:1150119/150‑1 (MQ=255)
                                                       cGTCACCTGACGGAAATCGATCGTCGCAGCGCCAACGGCGGACCGAAAGACGATTATTTTCGTTTTTTAGTCAT‑AAAAAAACTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACGAATGCGAGTATCTCCCTTATAAg                                                                                         <  2:1642772/150‑1 (MQ=255)
                                                       cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCGCACCAATGCGTGTATCGCCCTTATAAg                                                                                         >  2:118093/1‑150 (MQ=255)
                                                       cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGt                                                                                        >  2:532092/1‑150 (MQ=255)
                                                         tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGtg                                                                                       >  1:1817213/1‑149 (MQ=255)
                                                          cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCt                                                                                      >  1:1512009/1‑150 (MQ=255)
                                                           aCCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGCGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtg                                                                                     <  2:1113018/150‑1 (MQ=255)
                                                              tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCTTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                                  <  2:1138581/150‑1 (MQ=255)
                                                              tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGCGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                                  >  1:1227174/1‑150 (MQ=255)
                                                              tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                                  <  2:1616075/150‑1 (MQ=255)
                                                              tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                                  <  2:1220492/150‑1 (MQ=255)
                                                              tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                                  >  2:1492672/1‑150 (MQ=255)
                                                                 cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGTTTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCTGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACTCCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                                                                               >  1:302456/1‑150 (MQ=255)
                                                                 cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                                                                               <  1:26722/150‑1 (MQ=255)
                                                                  ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑GTTTTTTTTGCCTGTCGTGAGCGCACACCAATTCGAGTTTGTGCCGTATAAGTAGTTGGTTAc                                                                              >  2:1128142/1‑150 (MQ=255)
                                                                  ggAAAACGATCGTCGCCACGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                                                                              >  1:438046/1‑150 (MQ=255)
                                                                   gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                                                                             >  1:627314/1‑150 (MQ=255)
                                                                   gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTAAGTCAT‑AAAAAACCTGCCAGCGATGGCAGC‑‑TGATTTTTAGCCTGTAGTGAGCGCACACCAAGGCGACTATCTCCCTTATAAGTCTGTGGGTACg                                                                             >  2:256732/1‑150 (MQ=255)
                                                                      aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGACAGAGATGGCAGG‑‑TTTTTTTTTGCATGTCGTGAGCGCACACCAATGCGAGTATATCCCTTATAAGTCTGTGGTTACGTCa                                                                          >  2:1054158/1‑150 (MQ=255)
                                                                            cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                                                                    >  2:608901/1‑150 (MQ=255)
                                                                             gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagt                                                                   >  2:1688537/1‑149 (MQ=255)
                                                                             gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                                                                   >  1:1048999/1‑150 (MQ=255)
                                                                               cgccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                                                                 >  1:460051/1‑150 (MQ=255)
                                                                                 ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt                                                               >  1:1172555/1‑150 (MQ=255)
                                                                                 ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCATTATAAGTCTGTGGTTACGTCAATAGAGAGCtt                                                               >  1:1353560/1‑150 (MQ=255)
                                                                                  cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg                                                             <  1:532092/150‑1 (MQ=255)
                                                                                   gcgcCACCCGCTCACCGAAAGTCGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTTCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg                                                             <  2:1211605/150‑1 (MQ=255)
                                                                                   gcgcCAACCGCGCACCTAAAGACTTTTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCTTGGGAAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg                                                             <  2:1376206/150‑1 (MQ=255)
                                                                                   gcgcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTTTCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg                                                             >  1:1445970/1‑150 (MQ=255)
                                                                                   gcgcCAACCGCGCACCGAAAGACCATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGAAGGCAGG‑‑TTTTTTTTTGCCTGTAGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGATACGTCAATAGAGAGCTTTg                                                             <  2:1217707/150‑1 (MQ=255)
                                                                                    cgcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCGGCCAGCGTTGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt                                                            <  2:451599/150‑1 (MQ=255)
                                                                                     gcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGta                                                           >  1:1087204/1‑149 (MQ=255)
                                                                                      ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca                                                          <  2:688119/150‑1 (MQ=255)
                                                                                       cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTTGTTACGTCAATAGAGAGCTTTGTcac                                                         <  2:129765/150‑1 (MQ=255)
                                                                                         aCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa                                                       >  1:967363/1‑150 (MQ=255)
                                                                                          ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt                                                      <  1:2875/150‑1 (MQ=255)
                                                                                            gcgcACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTa                                                    >  1:1021633/1‑150 (MQ=255)
                                                                                              gcACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATc                                                  <  1:19792/150‑1 (MQ=255)
                                                                                              gaactgAATGACGATTATTTTCTTTTTTTAGTCTT‑AAAAAACCTGCCAGCGATGGAAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCGGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATc                                                  <  1:1269319/145‑1 (MQ=255)
                                                                                                 ccGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCTAGCGATGGCAGG‑‑TTTTTTTTTGCCTTTCGTGAGCGCAAACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGCGCTTAGTCACAATTATCTGc                                               >  2:289892/1‑150 (MQ=255)
                                                                                                  cGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCa                                              >  1:1855275/1‑150 (MQ=255)
                                                                                                   gAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCaa                                             <  1:1305563/150‑1 (MQ=255)
                                                                                                    aaaGACGATAAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTATGTGGTTACGTCAATAGATAGCTTTGTCACAATTATCTGCaaa                                            >  1:1186850/1‑150 (MQ=255)
                                                                                                              aaTTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGcc                                  >  1:1329533/1‑150 (MQ=255)
                                                                                                              aaTTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCGCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGcc                                  <  2:104191/150‑1 (MQ=255)
                                                                                                               aTTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCg                                 <  2:131165/150‑1 (MQ=255)
                                                                                                                  tCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTa                              >  2:718025/1‑150 (MQ=255)
                                                                                                                              cAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑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)
                                                                                                                                               |                                                                                                                                                       
CATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑‑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: ‑