Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,553,466 T→G 5.5% L81R (CTT→CGT) ‡ hemF → coproporphyrinogen III oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,553,4660TG 5.5% 310.1 / 6.4 104L81R (CTT→CGT) ‡hemFcoproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand):  ref base T (50/48);  new base G (3/3);  total (53/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 106 total aligned reads displayed.

TGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTG  >  NC_000913/2553317‑2553606
                                                                                                                                                     |                                                                                                                                            
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: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)
               gCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGttt                                                                                                                               <  2:1224329/150‑1 (MQ=255)
               gCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTAGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGAATGTCCACGGTGAGGCGATGCCGGCTTCCGCCACCGCTCATCGCCCGGACCTTGACGGGAGAAGttt                                                                                                                               >  2:234109/1‑150 (MQ=255)
                cAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTc                                                                                                                              >  1:1509210/1‑150 (MQ=255)
                cAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTc                                                                                                                              >  2:1131952/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)
                                        cgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTa                                                                                                      >  2:1659875/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)
                                                             ggtggtGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTAt                                                                                 <  1:1444505/150‑1 (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:560794/1‑150 (MQ=255)
                                                                                                    gTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGcc                                          >  2:642508/1‑150 (MQ=255)
                                                                                                      ccACGGTGAGGCGATGCCGGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGa                                        <  2:867571/150‑1 (MQ=255)
                                                                                                          ggTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCCCGCGAATGTGCGGTTTTTTATTGCCGaaaaa                                    >  1:800760/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)
                                                                                                               ggCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCggg                               >  1:970027/1‑150 (MQ=255)
                                                                                                                   aTGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGcc                           <  1:1165445/150‑1 (MQ=255)
                                                                                                                   aTGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGcc                           <  1:1581381/150‑1 (MQ=255)
                                                                                                                   aTGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGcc                           <  1:1131952/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)
                                                                                                                         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:159982/1‑150 (MQ=255)
                                                                                                                               gCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCtgg               <  2:1093668/150‑1 (MQ=255)
                                                                                                                                     gCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGttt         <  1:504590/150‑1 (MQ=255)
                                                                                                                                     gCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGttt         <  1:498100/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)
                                                                                                                                       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)
                                                                                                                                                     |                                                                                                                                            
TGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTG  >  NC_000913/2553317‑2553606

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: