Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,553,468 G→T 5.9% A82S (GCC→TCC)  hemF → coproporphyrinogen III oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,553,4680GT 5.9% 326.5 / 5.7 104A82S (GCC→TCC) hemFcoproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand):  ref base G (50/47);  new base T (3/4);  total (53/51)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 105 total aligned reads displayed.

TCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTT  >  NC_000913/2553319‑2553616
                                                                                                                                                     |                                                                                                                                                    
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)
        aGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGCAAGTTCCGGgcga                                                                                                                                              <  2:350702/150‑2 (MQ=255)
        aGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGgcgc                                                                                                                                              >  2:93717/1‑150 (MQ=255)
          ttGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGCAAGTTCCGGgcgatg                                                                                                                                            <  2:529404/150‑4 (MQ=255)
          ttGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGCAAGTTCCGGgcgatg                                                                                                                                            <  2:1114389/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)
                                          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)
                                                     cGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGAGATGTGCTTTTCACTGGTAGTGCATACGCATAAc                                                                                                 >  2:84257/1‑150 (MQ=255)
                                                       tAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAAccc                                                                                               <  1:1113157/150‑1 (MQ=255)
                                                        aaTGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTAGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACGGGTAGTGCATCCGCATAAcccc                                                                                              >  2:554663/1‑149 (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)
                                                                                   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)
                                                                                                        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)
                                                                                                             ggCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCggg                                         >  1:970027/1‑150 (MQ=255)
                                                                                                                 aTGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGcc                                     <  1:1581381/150‑1 (MQ=255)
                                                                                                                 aTGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGcc                                     <  1:1165445/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)
                                                                                                                      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:159982/1‑150 (MQ=255)
                                                                                                                             gCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCtgg                         <  2:1093668/150‑1 (MQ=255)
                                                                                                                                   gCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGttt                   >  2:409759/1‑150 (MQ=255)
                                                                                                                                   gCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGttt                   <  1:498100/150‑1 (MQ=255)
                                                                                                                                   gCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGttt                   <  1:504590/150‑1 (MQ=255)
                                                                                                                                   gCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGttt                   >  2:920392/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)
                                                                                                                                                    tGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACtt  >  2:1609148/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTT  >  NC_000913/2553319‑2553616

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: