Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,553,465 C→A 5.3% L81I (CTT→ATT) ‡ hemF → coproporphyrinogen III oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,553,4650CA 5.3% 332.4 / 4.5 106L81I (CTT→ATT) ‡hemFcoproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand):  ref base C (52/48);  new base A (2/4);  total (54/52)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

TTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTG  >  NC_000913/2553316‑2553606
                                                                                                                                                     |                                                                                                                                             
ttGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAAc                                                                                                                                               >  2:1238761/1‑150 (MQ=255)
 tGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACt                                                                                                                                              <  2:1384286/150‑1 (MQ=255)
  gTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACtt                                                                                                                                             >  1:1145399/1‑150 (MQ=255)
   tCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTg                                                                                                                                            >  2:219784/1‑150 (MQ=255)
     gAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGACCTTGcc                                                                                                                                          >  1:1594284/1‑150 (MQ=255)
     gAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACGTGcc                                                                                                                                          >  1:1668888/1‑150 (MQ=255)
     gAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCAGGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGcc                                                                                                                                          >  2:94457/1‑150 (MQ=255)
      aaGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCg                                                                                                                                         <  1:1612497/150‑1 (MQ=255)
           aGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGCAAGTTCCGGgcga                                                                                                                                    <  2:350702/150‑2 (MQ=255)
           aGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGgcgc                                                                                                                                    >  2:93717/1‑150 (MQ=255)
             ttGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGCAAGTTCCGGgcgatg                                                                                                                                  <  2:1114389/150‑4 (MQ=255)
             ttGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGCAAGTTCCGGgcgatg                                                                                                                                  <  2:529404/150‑4 (MQ=255)
              tGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGt                                                                                                                                 <  1:1309407/150‑1 (MQ=255)
              tGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGGTGCGTAATGGGGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGGCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGt                                                                                                                                 <  2:1104791/150‑1 (MQ=255)
               ggCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGGAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGtt                                                                                                                                <  1:964185/150‑1 (MQ=255)
                gCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTAGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGAATGTCCACGGTGAGGCGATGCCGGCTTCCGCCACCGCTCATCGCCCGGACCTTGACGGGAGAAGttt                                                                                                                               >  2:234109/1‑150 (MQ=255)
                 cAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTc                                                                                                                              >  2:1131952/1‑150 (MQ=255)
                 cAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTc                                                                                                                              >  1:1509210/1‑150 (MQ=255)
                    cgcgAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAg                                                                                                                           <  2:727752/150‑1 (MQ=255)
                           cTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATgg                                                                                                                    <  2:1413868/150‑1 (MQ=255)
                             ggcggcggcGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGc                                                                                                                  >  2:743899/1‑150 (MQ=255)
                                 gcggcggGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGttt                                                                                                              >  2:1619619/1‑150 (MQ=255)
                                     cgggcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACt                                                                                                          <  2:107298/150‑1 (MQ=255)
                                      gggcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGACCAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTg                                                                                                         <  1:1092817/150‑1 (MQ=255)
                                      gggcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTg                                                                                                         <  1:1605993/150‑1 (MQ=255)
                                       ggcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTgg                                                                                                        <  1:1169720/150‑1 (MQ=255)
                                        gcgTAGTCGGGTGTGGCGTAATGGGGGGGTTTTCGACCAGGCAGGCGGCAACTTTTCGCATGTCCACGGTGAGGCGAGGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGt                                                                                                       <  1:899033/150‑1 (MQ=255)
                                             gTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCTGTTTCGAGGCGATGGGCGGTTCACTGGTAGTGc                                                                                                  >  2:1599012/1‑150 (MQ=255)
                                              tCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCa                                                                                                 <  2:1171724/150‑1 (MQ=255)
                                              tCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCa                                                                                                 <  2:29223/150‑1 (MQ=255)
                                                  gtgtTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCg                                                                                             >  2:1165445/1‑150 (MQ=255)
                                                       gCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATaa                                                                                        >  1:78685/1‑150 (MQ=255)
                                                       gCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGATCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATaa                                                                                        <  2:1275785/150‑1 (MQ=255)
                                                        cGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGAGATGTGCTTTTCACTGGTAGTGCATACGCATAAc                                                                                       >  2:84257/1‑150 (MQ=255)
                                                          tAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAAccc                                                                                     <  1:1113157/150‑1 (MQ=255)
                                                           aaTGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTAGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACGGGTAGTGCATCCGCATAAcccc                                                                                    >  2:554663/1‑149 (MQ=255)
                                                            aTGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGt                                                                                   <  1:93717/150‑1 (MQ=255)
                                                              ggtggtGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCCCTGGTAGTGCATCCGCATAACCCGTAt                                                                                 >  2:102693/1‑150 (MQ=255)
                                                               gtggtgTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATg                                                                                <  2:1509210/150‑1 (MQ=255)
                                                                    gtTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTccc                                                                           <  2:628355/150‑1 (MQ=255)
                                                                     ttttCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCGCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCcca                                                                          <  2:1141988/150‑1 (MQ=255)
                                                                      tttCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCTACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCccac                                                                         >  2:498100/1‑150 (MQ=255)
                                                                      tttCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCTACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCccac                                                                         >  2:504590/1‑150 (MQ=255)
                                                                      tttCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGGTTCACTGGTAGTGCCTCCGCATACCCCGTATGTTCccac                                                                         >  2:151958/1‑150 (MQ=255)
                                                                       ttCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCccacc                                                                        <  2:78685/150‑1 (MQ=255)
                                                                         cGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAg                                                                      >  2:1522510/1‑150 (MQ=255)
                                                                          gAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGc                                                                     >  2:1234934/1‑150 (MQ=255)
                                                                                gcaggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCAcgcg                                                               <  2:19024/150‑1 (MQ=255)
                                                                                gcaggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCAcgcg                                                               >  2:647173/1‑150 (MQ=255)
                                                                                  aggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGaa                                                             >  1:1672448/1‑150 (MQ=255)
                                                                                    gcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAAtg                                                           <  2:462257/150‑1 (MQ=255)
                                                                                      gTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAAtgtg                                                         <  1:1659875/150‑1 (MQ=255)
                                                                                        cAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCg                                                       >  2:237781/1‑150 (MQ=255)
                                                                                          aCTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGt                                                     >  1:364543/1‑150 (MQ=255)
                                                                                           cTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGtt                                                    >  1:1198848/1‑150 (MQ=255)
                                                                                           cTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCAGGAACTTGCCGGGCGCAGGTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGtt                                                    <  2:723726/150‑1 (MQ=255)
                                                                                                cGCATGTCCACGGTGAGGCGATGCCTGCTTCCGACACCGCTCATCGCCCGGTACTTGCCGGGCGCAGTTTCGAGGGGAGGGGCGGTTCACTGGTAGTGCATCCGCCTAACCCGGATGTTCCCACCAGCCACGCGATAGTGCGGttttttt                                               >  1:1071943/1‑149 (MQ=255)
                                                                                                 gCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTAt                                              >  1:541114/1‑150 (MQ=255)
                                                                                                  cATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTAtt                                             >  1:439062/1‑150 (MQ=255)
                                                                                                   aTGTCCACGTTGAGTCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGACCTTGCCGGGCGCAGTTGCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTg                                            <  1:1214726/150‑1 (MQ=255)
                                                                                                    tGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGc                                           <  2:1657227/150‑1 (MQ=255)
                                                                                                     gTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGcc                                          >  2:642508/1‑150 (MQ=255)
                                                                                                       ccACGGTGAGGCGATGCCGGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGa                                        <  2:867571/150‑1 (MQ=255)
                                                                                                           ggTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCCCGCGAATGTGCGGTTTTTTATTGCCGaaaaa                                    >  1:800760/1‑150 (MQ=255)
                                                                                                           ggTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGaaaaa                                    >  1:1539819/1‑150 (MQ=255)
                                                                                                             tGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAAcc                                  >  1:1093668/1‑150 (MQ=255)
                                                                                                               aGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTTCCCGGCGCAGTTTCTAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCgg                                <  1:1238761/150‑1 (MQ=255)
                                                                                                               aGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCgg                                <  1:1361728/150‑1 (MQ=255)
                                                                                                               aGGCGATGCCTGCTTCCGCCACCGCTCATCGCCAGGAAATTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCgg                                <  2:1027086/150‑1 (MQ=255)
                                                                                                                    aTGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGcc                           <  1:1165445/150‑1 (MQ=255)
                                                                                                                    aTGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGcc                           <  1:1131952/150‑1 (MQ=255)
                                                                                                                    aTGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGcc                           <  1:1581381/150‑1 (MQ=255)
                                                                                                                    aTGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTAGAGGCGATGGGCGTTTCACTGGTAGAGCATCCGCATAACCCGTATGTTCCCACCAGCCCCGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGcc                           >  2:1148107/1‑150 (MQ=255)
                                                                                                                       ccTGCTTCCGCCACCGCTCATCGCCCGGTACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCTTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGAt                        <  1:19803/150‑1 (MQ=255)
                                                                                                                       ccTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCGCTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGAt                        >  1:1043369/1‑150 (MQ=255)
                                                                                                                       ccTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGAt                        >  1:789327/1‑150 (MQ=255)
                                                                                                                        cTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATc                       >  2:1199570/1‑150 (MQ=255)
                                                                                                                        cTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATc                       <  1:886890/150‑1 (MQ=255)
                                                                                                                        cTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTACTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATc                       <  1:94457/150‑1 (MQ=255)
                                                                                                                        cTGCTTCCGCCACCGATCATCGCACGGAACTTGCAGCGATAAGTTTCGAGGAGATGGTAGTTTCACTGGTAGTGCATACGCTTAAACCGTATGTTCACAAAAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTTCCCATc                       >  2:1199572/1‑150 (MQ=255)
                                                                                                                         tGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATcc                      >  1:971309/1‑150 (MQ=255)
                                                                                                                         tGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATcc                      <  2:541941/150‑1 (MQ=255)
                                                                                                                          gCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATccc                     <  1:647173/150‑1 (MQ=255)
                                                                                                                             tCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTc                  >  1:1118867/1‑150 (MQ=255)
                                                                                                                              ccGCCCCCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCCTCCGCAGACCCCGGAGGTTCCCACCAGCCACGCGATTGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGCTCGCGGCt                 >  1:88986/1‑150 (MQ=255)
                                                                                                                               cGCCACCGCTCATCGCCCGGAACTTTCCTGGCGAAGTTTAGAGGCGATGGGCGGTTCACAGGTAGTGCATCAGCAAAACACGTATGTTCCCACAAGCCAAGCGTATGTTAGGTTTTTTATTGCAGAAAAACATGGTGCATATGCCGtgtg                >  2:619972/1‑150 (MQ=255)
                                                                                                                               cGCCACCGCGCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCtg                <  2:568258/150‑1 (MQ=255)
                                                                                                                                gCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCtgg               <  2:1093668/150‑1 (MQ=255)
                                                                                                                                gCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCtgg               >  2:159982/1‑150 (MQ=255)
                                                                                                                                      gCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGttt         <  1:504590/150‑1 (MQ=255)
                                                                                                                                      gCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGttt         >  2:920392/1‑150 (MQ=255)
                                                                                                                                      gCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGttt         >  2:409759/1‑150 (MQ=255)
                                                                                                                                      gCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGAGGTTTTATATTGCCGAAAAACCGGGTGCCGTTCCCGTCTGGTGGttt         >  2:112579/1‑150 (MQ=255)
                                                                                                                                       cTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTg        <  1:743899/150‑1 (MQ=255)
                                                                                                                                        actgcGCCCGGCAAGTTCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTgg       >  1:641563/5‑150 (MQ=255)
                                                                                                                                         ctgcGCCCGGCAAGTTCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCCCCAGCCCCGCGAATGTGCGGTTTTGTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGCTGGTATGGc      >  2:1509048/4‑150 (MQ=255)
                                                                                                                                         cATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGc      <  2:1616416/150‑1 (MQ=255)
                                                                                                                                            cGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGt   <  1:409759/150‑1 (MQ=255)
                                                                                                                                                 ggAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTg  <  1:278767/146‑1 (MQ=255)
                                                                                                                                                 ggAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTg  >  2:278767/1‑146 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
TTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTG  >  NC_000913/2553316‑2553606

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: