Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,587,365 C→T S139S (AGC→AGT mrr → methylated adenine and cytosine restriction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,587,3650CT97.2% 1311.7 / 16.7 397S139S (AGC→AGTmrrmethylated adenine and cytosine restriction protein
Reads supporting (aligned to +/- strand):  ref base C (5/6);  new base T (192/193);  total (197/200)
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.36e-01

Only 100 of 406 total aligned reads displayed.

TTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGT  >  NC_000913/4587217‑4587512
                                                                                                                                                    |                                                                                                                                                   
ttGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGtcc                                                                                                                                                   >  5:53613/1‑151 (MQ=255)
 tGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCg                                                                                                                                                  >  5:470998/1‑151 (MQ=255)
   cGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGAGTTAGATCCGAAAGACGACTCTCAAGATCATGAAGAACTTGCAAAGAGCAGtcccg                                                                                                                                                 >  4:894008/1‑148 (MQ=255)
     tctcATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACg                                                                                                                                              >  8:1111164/1‑151 (MQ=255)
      ctcATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACg                                                                                                                                              >  2:218561/1‑150 (MQ=255)
        cATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACg                                                                                                                                              <  7:14375/148‑1 (MQ=255)
                    aTGACGGCGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAgg                                                                                                                               <  3:1223538/151‑1 (MQ=255)
                    aTGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAgg                                                                                                                               <  3:1446867/151‑1 (MQ=255)
                    aTGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAgg                                                                                                                               >  2:1634372/1‑151 (MQ=255)
                     tGACGGAGCAGGAGACGAACCATCTCGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGgc                                                                                                                              >  5:417945/1‑151 (MQ=255)
                          caGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAACCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTaa                                                                                                                         <  6:1204038/150‑1 (MQ=255)
                            ccAGGAGACGAACCATCTGCCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTaaaa                                                                                                                       <  1:1224024/150‑1 (MQ=255)
                             cAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTaaa                                                                                                                        >  8:212367/1‑149 (MQ=255)
                                 agaCGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCt                                                                                                                  >  1:399424/1‑151 (MQ=255)
                                   aCGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTc                                                                                                                >  8:173712/1‑151 (MQ=255)
                                     caaCCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGt                                                                                                              <  2:833660/150‑1 (MQ=255)
                                      aaCCATCTGGCCTTCGATTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTg                                                                                                             <  5:1528860/151‑1 (MQ=255)
                                        ccATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGAt                                                                                                           <  4:757597/151‑1 (MQ=255)
                                           tCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCg                                                                                                        >  3:416727/1‑151 (MQ=255)
                                            cTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCgg                                                                                                       >  2:45211/1‑151 (MQ=255)
                                               gCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTgg                                                                                                    >  4:168170/1‑151 (MQ=255)
                                               cccTTCGCTTTGGTGAATGTCAAACTTAAGGGACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGgggg                                                                                                    <  2:612517/150‑1 (MQ=255)
                                                 gttCGCTTTTGTGATTGTCAAAGTTAAGTCACGGCGGGATGCCGTCTATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCt                                                                                                  <  3:807174/150‑1 (MQ=255)
                                                   tCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCt                                                                                                  <  3:372925/149‑1 (MQ=255)
                                                    cGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCtgat                                                                                               <  8:541102/151‑1 (MQ=255)
                                                       ttttGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCtgatg                                                                                              >  5:310964/1‑149 (MQ=255)
                                                        tttGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAg                                                                                            <  8:382944/150‑1 (MQ=255)
                                                         ttGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGt                                                                                          >  5:916540/1‑151 (MQ=255)
                                                            tgAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCt                                                                                       >  8:1592356/1‑151 (MQ=255)
                                                               aTGTCAAACTTAAGGCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGa                                                                                    <  5:1355682/151‑1 (MQ=255)
                                                                 gTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATAAGGTTCTGGaa                                                                                   <  6:1000900/150‑1 (MQ=255)
                                                                  tCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTTGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGaaa                                                                                  >  8:727350/1‑150 (MQ=255)
                                                                  tCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGaaa                                                                                  >  3:237282/1‑150 (MQ=255)
                                                                   aaaaCTTAAGTCCCGGCCGGCTGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAAc                                                                                <  4:138/150‑1 (MQ=255)
                                                                    aaaCTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACt                                                                               <  8:1524975/151‑1 (MQ=255)
                                                                     aaCTTACGTCACGGCCGGATGCCGTCGATTTAGATCCTAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACtt                                                                              <  5:1233442/151‑1 (MQ=255)
                                                                     aaCTTAAGTCCCGGCCGGATGCCGTCTATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACtt                                                                              <  1:537487/151‑1 (MQ=255)
                                                                       cTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGAGCTAAAAGACCTTCGTGATCCGGTGGGTGATGAGGTTCTGGAAAACTTAt                                                                            >  3:1104556/1‑151 (MQ=255)
                                                                         tAAGTCCCGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTAtt                                                                           <  6:330138/150‑1 (MQ=255)
                                                                         tAAGTCACGGCCGGATGCGGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTAtt                                                                           <  7:382215/150‑1 (MQ=255)
                                                                         tAAGTCACCGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCGGATGAGGTTCTGGAAAACTTAtt                                                                           <  2:671814/150‑1 (MQ=255)
                                                                          aaGTCACGGCGGGATGCCGTCGATTTAGATCCGAAAGGCGACTCTCCCGATGATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTg                                                                          <  7:460455/150‑1 (MQ=255)
                                                                          aaGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGc                                                                         >  3:175201/1‑151 (MQ=255)
                                                                           cgTCACGGCCGGCTGCCTTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCa                                                                        <  2:1152505/150‑1 (MQ=255)
                                                                           aGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGc                                                                         >  7:93974/1‑150 (MQ=255)
                                                                           aGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGc                                                                         >  4:519156/1‑150 (MQ=255)
                                                                           aGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCa                                                                        <  7:565247/151‑1 (MQ=255)
                                                                             gccgggCCGGATGCGGTCTATTTAGAGCCCAACGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAgg                                                                      <  1:746434/147‑1 (MQ=255)
                                                                               aCGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGCCTCTCCCGATCATGAAGCACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTTAAAGAGCTTCGTGATGCGGTGGGTGATGAGGTTCTGGAACACTTATTGGAGGtt                                                                    >  6:1621756/1‑151 (MQ=255)
                                                                               aCGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGtt                                                                    >  6:1355895/1‑151 (MQ=255)
                                                                                cGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGttt                                                                   <  6:740913/151‑1 (MQ=255)
                                                                                cGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGttt                                                                   <  5:951863/151‑1 (MQ=255)
                                                                                cGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGttt                                                                   >  7:821198/1‑151 (MQ=255)
                                                                                cGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGttt                                                                   >  4:1116633/1‑151 (MQ=255)
                                                                                 ggCCGGATGCCGTCTATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTtc                                                                  <  4:604778/151‑1 (MQ=255)
                                                                                 ggCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTtc                                                                  >  7:1319818/1‑151 (MQ=255)
                                                                                 ggCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTtc                                                                  >  1:1523455/1‑151 (MQ=255)
                                                                                 ggCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTtc                                                                  >  3:1384056/1‑151 (MQ=255)
                                                                                  gCCGGATGCCGTCGATTTAGATCCGAACGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTtct                                                                 <  1:1683809/151‑1 (MQ=255)
                                                                                   ccGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTtctc                                                                <  2:274807/151‑1 (MQ=255)
                                                                                   ccGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTtctc                                                                <  3:565237/151‑1 (MQ=255)
                                                                                       aTGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGc                                                                                                   >  2:1432386/1‑112 (MQ=255)
                                                                                       aTGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTc                                                            <  5:758872/151‑1 (MQ=255)
                                                                                       aTGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTc                                                            <  3:977210/151‑1 (MQ=255)
                                                                                         gCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTcg                                                           <  2:320734/150‑1 (MQ=255)
                                                                                          ccGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCCAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTcgcg                                                         <  6:1110075/151‑1 (MQ=255)
                                                                                            gTCTATTAAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCt                                                       <  4:133630/151‑1 (MQ=255)
                                                                                               gATTTAGATCAGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGACCAGTACGCACGATAGGTTTGATAAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGTTTCTGGAAAACTTATTGCAGGTTTCTCATTcccgggttc                                                    >  6:565902/1‑146 (MQ=255)
                                                                                                aTTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTg                                                    <  1:472496/150‑1 (MQ=255)
                                                                                                    acaTCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGt                                                <  5:1009457/148‑1 (MQ=255)
                                                                                                     gATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTc                                               >  5:840498/1‑150 (MQ=255)
                                                                                                      aTCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTGAt                                             >  3:1047552/1‑151 (MQ=255)
                                                                                                      aTCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCAt                                             >  7:1146332/1‑151 (MQ=255)
                                                                                                      aTCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAAATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCAt                                             <  5:107552/151‑1 (MQ=255)
                                                                                                        cgtaacgCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTg                                           <  6:1493844/145‑1 (MQ=255)
                                                                                                         cGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGt                                          >  7:1080913/1‑151 (MQ=255)
                                                                                                                cGACTCTCCCGATCATGAAGAACTTACAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGAt                                   >  7:708681/1‑151 (MQ=255)
                                                                                                                    tctcCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGttt                               >  6:691600/1‑151 (MQ=255)
                                                                                                                     gtcCCGATCATGAAGAACTTGCAAAGAGGAGCGCGGACGATCGGTTAGATCAGGCGGTAAAAGATCTTCGTGATGCGGTGGCTGATGAGTTTCTGGATAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGttt                               <  4:548778/149‑1 (MQ=255)
                                                                                                                     ctcCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGtt                                <  5:207493/149‑1 (MQ=255)
                                                                                                                       cccGATCATGAAGCACTTGCAAAGAGCAGTCCGGACGATCGGGTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGtttt                              <  6:870332/149‑1 (MQ=255)
                                                                                                                       cccGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAATAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTGATTGTTCTGGATGttt                               <  6:197433/148‑1 (MQ=255)
                                                                                                                       cccGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGtttt                              <  4:453912/149‑1 (MQ=255)
                                                                                                                       cccGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGtttt                              <  1:988259/149‑1 (MQ=255)
                                                                                                                       cccGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTg                             <  4:282621/150‑1 (MQ=255)
                                                                                                                        ccGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCa                           <  7:167195/151‑1 (MQ=255)
                                                                                                                         cGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCAt                          >  7:1537677/1‑151 (MQ=255)
                                                                                                                            tCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGc                       >  2:136619/1‑151 (MQ=255)
                                                                                                                             cATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGcc                      >  3:1577820/1‑151 (MQ=255)
                                                                                                                              aTGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCt                     >  1:610770/1‑151 (MQ=255)
                                                                                                                               tGAAGAACTTGCAAAGAGCAGTCCGGACTATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTATGGATGTTTTGAATAGCCTg                    >  2:1487425/1‑151 (MQ=255)
                                                                                                                               tGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTg                    >  3:639969/1‑151 (MQ=255)
                                                                                                                               tGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATCCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTCTTCTGGATGTTTTCCATCGCCTg                    >  1:1442745/1‑151 (MQ=255)
                                                                                                                                gaagaaCTAGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTg                    <  2:905688/150‑1 (MQ=255)
                                                                                                                                 aagaaCTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTggg                  <  6:173465/151‑1 (MQ=255)
                                                                                                                                       ttGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATgg            >  1:109516/1‑151 (MQ=255)
                                                                                                                                        tGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATc                         >  1:885404/1‑137 (MQ=255)
                                                                                                                                        tGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATg             <  8:1310695/149‑1 (MQ=255)
                                                                                                                                              gagCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCAc     <  1:1007456/151‑1 (MQ=255)
                                                                                                                                               atcagtCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGGGGCTGATGAGGTTCTGGAAAAGTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCAcc    <  7:1192415/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGT  >  NC_000913/4587217‑4587512

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