Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,264 A→C 5.2% intergenic (+34/‑183) 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,2640AC 5.2% 348.7 / 6.2 122intergenic (+34/‑183)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (63/50);  new base C (4/3);  total (67/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 123 total aligned reads displayed.

CATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGAT  >  NC_000913/4060119‑4060413
                                                                                                                                                  |                                                                                                                                                       
cATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                                     >  2:1014727/1‑150 (MQ=255)
 aTCGATGATGACGAACTGGGAGCAGTGACTCCGACCTCTATCCGTATGCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCa                                                                                                                                                    <  1:899135/150‑1 (MQ=255)
   cGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTAATCGTAAACGTCACCTGACGGACAACGATCGTCGCCGCGACACCCGCGCACCGAAAGCCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGAGATGGCAgg                                                                                                                                                  >  2:1853635/1‑150 (MQ=255)
    gatgatGACGAACTGGTAGAAGTGACTCCTTCCTCTATCCGTTTTCGTAAACGGCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGGAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑t                                                                                                                                               <  2:1688315/150‑1 (MQ=255)
    gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑‑t                                                                                                                                               <  1:204450/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:1568950/1‑150 (MQ=255)
           aCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑tttttttt                                                                                                                                        <  1:711949/150‑1 (MQ=255)
            cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑ttttttttt                                                                                                                                       >  1:182504/1‑150 (MQ=255)
            cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATATCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑ttttttttt                                                                                                                                       >  2:494737/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: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:254744/150‑1 (MQ=255)
                              tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGcaca                                                                                                                     >  2:1346449/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)
                                            gTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                                                                       <  1:1448138/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)
                                                   tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                                                                                >  2:592801/1‑150 (MQ=255)
                                                      aCGTCACCAGACGGAAATAGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                                                                             <  2:1150119/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)
                                                              tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGCGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                                     >  1:1227174/1‑150 (MQ=255)
                                                              tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                                     <  2:1616075/150‑1 (MQ=255)
                                                              tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                                     >  2:1492672/1‑150 (MQ=255)
                                                              tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                                     <  2:1220492/150‑1 (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)
                                                                   gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                                                                                >  1:627314/1‑150 (MQ=255)
                                                                   gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTAAGTCATAAAAAACCTGCCAGCGATGGCAGC‑‑TGATTTTTAGCCTGTAGTGAGCGCACACCAAGGCGACTATCTCCCTTATAAGTCTGTGGGTACg                                                                                >  2:256732/1‑150 (MQ=255)
                                                                      aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGACAGAGATGGCAGG‑‑TTTTTTTTTGCATGTCGTGAGCGCACACCAATGCGAGTATATCCCTTATAAGTCTGTGGTTACGTCa                                                                             >  2:1054158/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)
                                                                                 ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGACTGTCGTGCGCGCACACCAATGCGAGTATCTACCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt                                                                  >  1:349851/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)
                                                                                          ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATTAGTCTGTGGTGACGTCAATAGAGAGCTTTGTCACAAt                                                         >  2:1041431/1‑150 (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)
                                                                                                   gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCa                                                 >  1:828242/1‑150 (MQ=255)
                                                                                                       gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTc                                            <  2:1128937/150‑1 (MQ=255)
                                                                                                           aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATAc                                        <  1:570221/150‑1 (MQ=255)
                                                                                                              aaTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGcc                                     >  1:1329533/1‑150 (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)
                                                                                                                                tAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCttttt                   >  1:1201401/1‑150 (MQ=255)
                                                                                                                                 aaaaaaaCCTGCCATCGCTGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCttttt                   >  2:1545287/2‑150 (MQ=255)
                                                                                                                                  aaaaaaCCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGGTTTCtttttt                  >  1:173814/1‑150 (MQ=255)
                                                                                                                                   aaaaaCCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTg                 <  2:1473544/150‑1 (MQ=255)
                                                                                                                                   aaaaaCCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTg                 <  2:572512/150‑1 (MQ=255)
                                                                                                                                   aaaaaCCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTg                 <  1:1178418/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)
                                                                                                                                    aaaaCCTGCCAGCGAGGTCAGG‑‑TTTTTTTTTGCCTGTCGTGAGAGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTgg                <  2:903945/150‑1 (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)
                                                                                                                                               gCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaagcgtaa     >  1:1316612/1‑150 (MQ=255)
                                                                                                                                                 gATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGa   <  2:646328/150‑1 (MQ=255)
                                                                                                                                                  aTGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGAt  <  2:577011/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                       
CATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGAT  >  NC_000913/4060119‑4060413

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: