Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,584,054 A→C L903R (CTC→CGC)  hsdR ← endonuclease R Type I restriction enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,584,0540AC97.1% 233.9 / ‑3.7 69L903R (CTC→CGC) hsdRendonuclease R Type I restriction enzyme
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base C (32/35);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-01

TGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGAGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTC  >  NC_000913/4583907‑4584199
                                                                                                                                                   |                                                                                                                                                 
tGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCt                                                                                                                                                >  2:440182/1‑151 (MQ=255)
  cGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTc                                                                                                                                              <  1:251843/151‑1 (MQ=255)
   ggCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCg                                                                                                                                             >  1:163127/1‑151 (MQ=255)
    gCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGc                                                                                                                                            >  1:619465/1‑151 (MQ=255)
       cTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGccac                                                                                                                                         >  1:306102/1‑151 (MQ=255)
         ttttCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGccacc                                                                                                                                        >  1:904313/1‑150 (MQ=255)
            tCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAgt                                                                                                                                     <  1:493314/150‑1 (MQ=255)
                 cAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTc                                                                                                                                >  2:723962/1‑150 (MQ=255)
                  aGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCAt                                                                                                                              >  2:97411/1‑151 (MQ=255)
                   ggCGCGAGGCGAAGCCGTTAACGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCGAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGATCGATGATACGCTGGCGCTTCGCAACCAGTTGTTCAt                                                                                                                              >  2:15397/1‑150 (MQ=255)
                   ggCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATg                                                                                                                             <  1:449915/151‑1 (MQ=255)
                        gAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCtgt                                                                                                                        <  1:701645/151‑1 (MQ=255)
                                    ttAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAAc                                                                                                            >  1:764197/1‑151 (MQ=255)
                                      aaaGTTCACGCCCGCCTCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACt                                                                                                           <  1:363266/150‑1 (MQ=255)
                                            cACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCa                                                                                                    <  1:618949/151‑1 (MQ=255)
                                              cgcccgccGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCAt                                                                                                   <  2:253246/150‑1 (MQ=255)
                                               gcccgccGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATcc                                                                                                 <  2:907875/151‑1 (MQ=255)
                                                cccgccGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCg                                                                                                <  2:644788/151‑1 (MQ=255)
                                                cccgccGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCg                                                                                                >  1:838133/1‑151 (MQ=255)
                                                    ccgcGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTc                                                                                            <  1:364424/151‑1 (MQ=255)
                                                      gcgTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCgg                                                                                          <  2:809516/151‑1 (MQ=255)
                                                       cgTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGt                                                                                         >  2:683486/1‑151 (MQ=255)
                                                          ccTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGAt                                                                                      <  1:184518/151‑1 (MQ=255)
                                                             ggCATAGCTCATCCAGACGTCGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGGTGTTCATGGCTGTGCTCGGCAAATCTGCGGCCAGCCGCTTCGGTGAtttt                                                                                   <  2:899935/151‑1 (MQ=255)
                                                             ggCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGAtttt                                                                                   <  2:528666/151‑1 (MQ=255)
                                                             ggCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGAGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGAtttt                                                                                   <  1:235860/151‑1 (MQ=255)
                                                              gCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGAttt                                                                                    <  1:508576/149‑1 (MQ=255)
                                                                aTAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTtat                                                                                 >  1:197200/1‑150 (MQ=255)
                                                                  aGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAg                                                                               <  1:300322/150‑1 (MQ=255)
                                                                     tCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGttt                                                                           >  1:161078/1‑151 (MQ=255)
                                                                      cATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGtt                                                                            <  1:866344/149‑1 (MQ=255)
                                                                       aTCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGCGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCt                                                                         >  2:495441/1‑151 (MQ=255)
                                                                         ccAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTg                                                                        >  1:136567/1‑150 (MQ=255)
                                                                         ccAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGa                                                                       <  1:234050/151‑1 (MQ=255)
                                                                         ccAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGa                                                                       <  1:539391/151‑1 (MQ=255)
                                                                         ccAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGa                                                                       <  2:378329/151‑1 (MQ=255)
                                                                          cAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGaa                                                                      <  2:735638/151‑1 (MQ=255)
                                                                            gacgacGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATc                                                                    >  1:145297/1‑151 (MQ=255)
                                                                              cgacgCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCgg                                                                  <  1:391440/151‑1 (MQ=255)
                                                                                  gCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAAt                                                              >  2:759736/1‑151 (MQ=255)
                                                                                        gTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCAtt                                                        <  2:608379/151‑1 (MQ=255)
                                                                                        gTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCAtt                                                        <  1:727132/151‑1 (MQ=255)
                                                                                          ttATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGAGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTg                                                       <  1:418696/150‑1 (MQ=255)
                                                                                                atcGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCGATCCGCTTAGCTGATTTTTTAGGTTTCTGAATCGGTAATTTGATTGACCAg                                                  >  2:835371/1‑149 (MQ=255)
                                                                                                atcGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCa                                                   <  2:245623/148‑1 (MQ=255)
                                                                                                atcGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCa                                                   <  2:203803/148‑1 (MQ=255)
                                                                                                atcGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCa                                                   <  1:754638/148‑1 (MQ=255)
                                                                                                atcGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGc                                                 <  2:601945/150‑1 (MQ=255)
                                                                                                   gTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGttt                                             >  2:1462/1‑151 (MQ=255)
                                                                                                   gTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGttt                                             >  1:763338/1‑151 (MQ=255)
                                                                                                     ttCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGttt                                             >  2:14551/1‑149 (MQ=255)
                                                                                                     ttCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGc                                           >  2:133110/1‑151 (MQ=255)
                                                                                                        gCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGt                                        >  1:129984/1‑151 (MQ=255)
                                                                                                          tGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTc                                      >  2:173745/1‑151 (MQ=255)
                                                                                                           gCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGtt                                       <  1:387849/149‑1 (MQ=255)
                                                                                                           gCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTc                                      <  2:689328/150‑1 (MQ=255)
                                                                                                                cACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCtt                                >  2:238744/1‑151 (MQ=255)
                                                                                                                cACGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCtt                                <  2:837734/151‑1 (MQ=255)
                                                                                                                 aCGGTTAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTc                               >  2:460070/1‑151 (MQ=255)
                                                                                                                      tAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGgcg                          >  2:448881/1‑151 (MQ=255)
                                                                                                                      tAAACGTGGCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGgcg                          >  1:771931/1‑151 (MQ=255)
                                                                                                                             ggCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGcaccacc                   >  2:776266/1‑151 (MQ=255)
                                                                                                                              gCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGcaccacca                  <  1:448891/151‑1 (MQ=255)
                                                                                                                              gCCAGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGcaccacca                  <  1:400468/151‑1 (MQ=255)
                                                                                                                                 aGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCg               >  2:82084/1‑151 (MQ=255)
                                                                                                                                 aGACCGATGATACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGAACCACCACCg               >  2:770294/1‑151 (MQ=255)
                                                                                                                                        tgatACGCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCa         >  2:797318/1‑150 (MQ=255)
                                                                                                                                              gCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTc  <  2:163127/151‑1 (MQ=255)
                                                                                                                                              gCTGGCGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTc  <  1:83111/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGATGATACGCTGGAGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTC  >  NC_000913/4583907‑4584199

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: