Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,067,560 T→C 5.9% N174S (AAC→AGC)  yggE ← uncharacterized protein YggE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,067,5600TC 5.9% 210.0 / 3.0 68N174S (AAC→AGC) yggEuncharacterized protein YggE
Reads supporting (aligned to +/- strand):  ref base T (36/28);  new base C (2/2);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

ATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGACCCGCCTTCAGCGCGCCATCCAGCAAGGAATTCAG  >  NC_000913/3067412‑3067708
                                                                                                                                                    |                                                                                                                                                    
atcaAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGtt                                                                                                                                                     >  2:154444/1‑150 (MQ=255)
atcaAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGtt                                                                                                                                                     >  1:203374/1‑150 (MQ=255)
    aaCTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGcc                                                                                                                                                 <  1:323484/150‑1 (MQ=255)
     aCTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCa                                                                                                                                                <  1:1356632/150‑1 (MQ=255)
       tGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGt                                                                                                                                              >  2:712338/1‑150 (MQ=255)
        gAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGtt                                                                                                                                             >  2:497633/1‑150 (MQ=255)
         aaTAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTc                                                                                                                                            >  2:657343/1‑150 (MQ=255)
                    gCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCtgat                                                                                                                                 <  2:511225/150‑1 (MQ=255)
                    gCTCGTAAGGTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGTGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCtgat                                                                                                                                 <  2:1087459/150‑1 (MQ=255)
                     cTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCtgatg                                                                                                                                <  2:1000206/150‑1 (MQ=255)
                         tAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATc                                                                                                                            <  1:853945/150‑1 (MQ=255)
                            gTTTCCTGGGCGGACACAGGCGCGGCATCGGCTTTAATCATCCGCACCATTGGGCTGGGCTGATAGTCTGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATcgcg                                                                                                                         <  2:902641/150‑1 (MQ=255)
                                ccTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTAt                                                                                                                     >  1:1322367/1‑150 (MQ=255)
                                     gCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATc                                                                                                                <  1:4115/150‑1 (MQ=255)
                                      cGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCg                                                                                                               >  1:394877/1‑150 (MQ=255)
                                       ggACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGc                                                                                                              >  2:1292471/1‑150 (MQ=255)
                                            ccGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCt                                                                                                         >  2:753903/1‑150 (MQ=255)
                                            ccGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCt                                                                                                         <  1:657343/150‑1 (MQ=255)
                                             cGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCtt                                                                                                        <  2:249886/150‑1 (MQ=255)
                                                     cATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCtt                                                                                                <  2:624310/150‑1 (MQ=255)
                                                       tCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTg                                                                                              <  2:876214/150‑1 (MQ=255)
                                                         ggCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTc                                                                                            <  2:241254/150‑1 (MQ=255)
                                                           cTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCtt                                                                                          >  1:1128435/1‑150 (MQ=255)
                                                               catcatCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTtata                                                                                      <  2:1216714/150‑1 (MQ=255)
                                                                 tcatcCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCACACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAgg                                                                                    >  2:63190/1‑150 (MQ=255)
                                                                       gCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCg                                                                              >  2:226339/1‑150 (MQ=255)
                                                                             ttGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCCGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTgcg                                                                        <  1:497633/150‑1 (MQ=255)
                                                                               gggctgggctGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGcc                                                                      >  2:568581/1‑150 (MQ=255)
                                                                                ggctgggctgATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCa                                                                     >  2:111031/1‑150 (MQ=255)
                                                                                ggctgggctgATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCa                                                                     >  2:153618/1‑150 (MQ=255)
                                                                                    gggctgATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGcc                                                                 <  2:1391304/150‑1 (MQ=255)
                                                                                      gctgATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGCTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTg                                                                                              >  1:49935/1‑119 (MQ=255)
                                                                                      gctgATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGCTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTg                                                                                              <  2:55189/119‑1 (MQ=255)
                                                                                      gctgATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGCTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTg                                                                                              >  1:55189/1‑119 (MQ=255)
                                                                                      gctgATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGCTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTg                                                                                              <  2:49935/119‑1 (MQ=255)
                                                                                       ctgATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAg                                                              >  1:283764/1‑150 (MQ=255)
                                                                                       ctcaTAGTTGGAAACATGGTAGCGCACGCTATAAACCGCACCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCAGCGCCTGATGAATCGAGGTATCAATCGCTGCCTTACGCGCTATGTCTTTAAAGGCATCCGGATGCGCCACGCCCac                                                              >  1:758187/4‑149 (MQ=255)
                                                                                         gATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTGTAGGCATCCGGCTGCGCCACGCCCAGCg                                                            <  1:85611/150‑1 (MQ=255)
                                                                                         gATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCg                                                            <  1:154444/150‑1 (MQ=255)
                                                                                           tAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGac                                                          >  1:631330/1‑150 (MQ=255)
                                                                                            aGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGTTGAAAGCCGTTCGCCAGTTCCTGAGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGaca                                                         >  2:1163718/1‑150 (MQ=255)
                                                                                             gTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGacac                                                        >  1:178162/1‑150 (MQ=255)
                                                                                              ttGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGAcaca                                                       <  2:1395624/150‑1 (MQ=255)
                                                                                                 gAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGaa                                                    >  2:827775/1‑150 (MQ=255)
                                                                                                     cATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGaa                                                >  1:824663/1‑150 (MQ=255)
                                                                                                          tAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCg                                           <  1:753903/150‑1 (MQ=255)
                                                                                                           aGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTc                                            >  2:796100/1‑148 (MQ=255)
                                                                                                           aGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTc                                            <  1:796100/148‑1 (MQ=255)
                                                                                                            gcgcACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGtt                                         <  1:701002/150‑1 (MQ=255)
                                                                                                            gcgcACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGtt                                         <  2:283764/150‑1 (MQ=255)
                                                                                                              gcACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCCGCGACACAGAACGAATTTCGTTaa                                       >  2:1186385/1‑150 (MQ=255)
                                                                                                                aCGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGa                                     >  2:512947/1‑150 (MQ=255)
                                                                                                                 cGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGAc                                    >  2:1150012/1‑150 (MQ=255)
                                                                                                                  gCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGAcc                                   >  1:1054878/1‑150 (MQ=255)
                                                                                                                   cTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGAccc                                  <  1:768057/150‑1 (MQ=255)
                                                                                                                    tatataCCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGACCCg                                 >  2:1098517/1‑150 (MQ=255)
                                                                                                                      tataCCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGACCCGcc                               <  1:542420/150‑1 (MQ=255)
                                                                                                                      tataCCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGACCCGcc                               <  2:1322367/150‑1 (MQ=255)
                                                                                                                       ataCCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGACCCGCCt                              >  2:479123/1‑150 (MQ=255)
                                                                                                                       ataCCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGACCCGCCt                              >  1:1336447/1‑150 (MQ=255)
                                                                                                                         aCCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCAGGCCCAGCGACACAGAACGAATTTCGTTAAGACCCGCCTTc                            <  2:955010/150‑1 (MQ=255)
                                                                                                                         aCCGCCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGACCCGCCTTc                            >  2:337398/1‑150 (MQ=255)
                                                                                                                            ggCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGACCCGCCTTCAgc                         >  1:606704/1‑150 (MQ=255)
                                                                                                                                  aGTTAACGATGAAAGCCGTTCGCCAGTACCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGACCCGCCTTCAGCGCGCCa                   >  1:442361/1‑150 (MQ=255)
                                                                                                                                     ttACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGACCCGCCTTCAGCGCGCCATcc                <  1:63190/150‑1 (MQ=255)
                                                                                                                                     ttACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATAGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGACCCGCCTTCAGCGCGCCATcc                >  1:159441/1‑150 (MQ=255)
                                                                                                                                           tGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGACCCGCCTTCAGCGCGCCATCCAGCAAg          >  1:7666/1‑150 (MQ=255)
                                                                                                                                           tGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGACCCGCCTTCAGCGCGCCATCCAGCAAg          >  1:818479/1‑150 (MQ=255)
                                                                                                                                                  cGTTCGCCAGTTTCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGACCCGCCTTCAGCGCGCCATCCAGCAAGGAATTcc   >  2:974399/1‑149 (MQ=255)
                                                                                                                                                   gTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGACCCGCCTTCAGCGCGCCATCCAGCAAGGAATTCAg  >  2:577951/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGCTTTCATCATCCGCACCATTGGGCTGGGCTGATAGTTGGAAACATGGTAGCGCACGCTATATACCGGCCCCAGTTTACGATGAAAGCCGTTCGCCAGTTCCTGCGCCTGATGAATCGCGTTATCAATCGCTGCCTTACGCGCTTTGTCTTTATAGGCATCCGGCTGCGCCACGCCCAGCGACACAGAACGAATTTCGTTAAGACCCGCCTTCAGCGCGCCATCCAGCAAGGAATTCAG  >  NC_000913/3067412‑3067708

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