Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,147,196 C→T P135S (CCG→TCG)  yceD → DUF177 domain‑containing protein YceD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,147,1960CT96.2% 154.7 / ‑1.4 52P135S (CCG→TCG) yceDDUF177 domain‑containing protein YceD
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (23/27);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.74e-01

GTTAAACGGCGATGCGAAGGTGACGGTAACGCTCGAGTGTCAGCGTTGCGGGAAGCCGTTTACTCATCAGGTCTACACAACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTCCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTCCATGGCCGTACAACAGAATAAACCAACCCGTTCCAAACGTGGCATGCGTCGTTCCCATGACGCGCTGACCGC  >  NC_000913/1146973‑1147437
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 
gtTAAACGGCGATGCGAAGGTGACGGTAACGCTCGAGTGTCAGCGTTGCGGGAAGCCGTTTACTCATCAGGTCTACACAACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGaa                                                                                                                                                                                                                         <  1:535386/250‑1 (MQ=255)
gatAAACGGCGATGCGAAGGAGACGGGAACGATCGAGTGTCAGCGTTGCGGGAAGACGTTTATTCAGCAGGTCTACACAACGTATTGTTTTAGTACTGTGCGTTCAGACGAACAGGCTGAGGCACTGCTGGAAACGTATGAACCGATTGAGGTTAACGGATTCGGTGAAATCGATCTGCTTGCAATGGTTGAGGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGCTTCTGAACACTGTGa                                                                                                                                                                                                                          >  2:407182/3‑249 (MQ=255)
        gcgagtcGAAGGTGACGGGAACGATCGAGTTGCAGCGGTGCGGGAAGCCGTTTACTCATCAGGTCTACATAACGTATTGTTTTAGTCCTGTGCGTTTAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGaa                                                                                                                                                                                                                <  2:94667/245‑1 (MQ=255)
        gCGATGCGAAGGTGACGGTAACGCTCGAGTGTCAGCGTTGCGGGAAGCCGTTTACTCATCAGGTCTACACAACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGTTGCATGATTCTGAACACTGTGAAGTGTCCGaa                                                                                                                                                                                                                <  1:417520/251‑1 (MQ=255)
                     gaccgTAACGCTCGAGTGTCAGCGTTGCTGGAAGCCGTTTACTCATCAGGTCTACACAACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGTTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCt                                                                                                                                                                                                   <  1:274603/247‑1 (MQ=255)
                      aCGGTAACGCTCGAGTGTCAGCGTTGCGGGAAGCCGTTTACTCATCAGGTCTACACAACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCtt                                                                                                                                                                                                  >  1:80122/1‑251 (MQ=255)
                       cGGTAACGCCCGAGTGTCAGCGTTGCGGGAAGCCGTTTACTCATCAGGTCTACACAACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCttt                                                                                                                                                                                                 <  2:189789/251‑1 (MQ=255)
                          tAACGCTCGAGTGTCAGCGTTGCGGGAAGACGTTTACTCATCAGGTCTACACAACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGCAGATGAAATCATCCTCGCCTTGCCAGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGt                                                                                                                                                                                              >  1:146383/1‑251 (MQ=255)
                                gcGAGTGTCAGCGTTGCGGGAAGTCGTTTTCTCATCAGGGTTACACAACGTATTGTTTTAGTCCTGTGCTGTCAGACGAACAGGCTTAAGACCTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGAGCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACTCTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACt                                                                                                                                                                                         <  2:423580/249‑1 (MQ=255)
                                 cGAGTGTCAGCGTTGCGGGAAGCCGTTTACTCATCAGGTCTACACAACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGc                                                                                                                                                                                       >  1:112694/1‑251 (MQ=255)
                                 cGAGTGTCAGCGTTGCGGGAAGCCGTTTACTCATCAGGTCTACACAACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGc                                                                                                                                                                                       >  1:482730/1‑251 (MQ=255)
                                   agtttCAGTGTTGCGGGATGTCGTTTACTCTACAGTTTTACTCAACGTAGTGTTTTAGTCCTGGGCGTCTAGTCGAACAGCCTGAAGCACTGCCGGAACCGTAGGTACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGACATCATCCTCGCCTGGCCGGTAGTTTCGGTGCATGATTCTTAACACTGTGAATTGTCCGAAGCGGACATGGTCTTTGTTCAACTGCCt                                                                                                                                                                                     <  2:208311/247‑1 (MQ=255)
                                           cGTTGCGGGAAGCCGTTTACTCATCAGGTCTACACAACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAgc                                                                                                                                                                             >  2:127830/1‑251 (MQ=255)
                                            gTTGCTGGAAGCCGTTTACTCATCAGGTCTACACAACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGTCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCTCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTTTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAgcg                                                                                                                                                                            <  1:10711/251‑1 (MQ=255)
                                                           ttACTCATCAGGTCTACACAACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCTCTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAAccc                                                                                                                                                              <  1:196951/250‑1 (MQ=255)
                                                             aCTCATCAGGTCTACACAACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCAtt                                                                                                                                                           <  1:400236/251‑1 (MQ=255)
                                                                catcaGGTCTACACAACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGc                                                                                                                                                        >  1:132595/1‑251 (MQ=255)
                                                                catcaGGCCTACACAACGTATTGTTTTCGTCCTGGGCGTTCAGACGAACAGGCTGAAGTACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGACATCATCCTCGCCTTGCCTGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCCTTGGTGAACTGCCTTACGAAGCGCAAAAGCCAAGCCCATTTGc                                                                                                                                                        <  1:461119/251‑1 (MQ=255)
                                                                 atcaGGTCTACACAACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTCCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGcc                                                                                                                                                       <  1:210816/251‑1 (MQ=255)
                                                                            acaACGTTTGGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTTAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAg                                                                                                                                            <  2:159107/251‑1 (MQ=255)
                                                                                   aTTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCAGTATTAGCCAGCTTaaa                                                                                                                                      <  1:375823/250‑1 (MQ=255)
                                                                                    ttGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGACTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCAAGCTTAAAGc                                                                                                                                    >  1:13450/1‑251 (MQ=255)
                                                                                         ttAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCgta                                                                                                                                 >  1:269700/1‑249 (MQ=255)
                                                                                                 gtgCGTTCAGACTAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGGGATCGGCCCGCAGAAGCGCAATAGCCAAACCAATTTGCCGTATGAGCAAGCTTAAAGCGAAAGTAATTggg                                                                                                                       >  2:92091/1‑250 (MQ=255)
                                                                                                 gtgCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGt                                                                                                                       <  1:127830/251‑1 (MQ=255)
                                                                                                          gACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGt                                                                                                              >  1:475475/1‑251 (MQ=255)
                                                                                                                cAGGCTGAAGTGATGCCGGAAGCGTTTTAACCGATTAAGGTTAACGATTGCTGTGCAATCGATCTGTTGTCAATGGATGAAGATGAAACCATCCTCGCGTTGCCGGTAGTGTCGGTGCTTGGTTCGGAACACTGTGAAGTGTCCGAAGCGTACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCTAACCCATTTGCCGTATTAGCCAGCATAAAGCGTAAGTAATTGGTGCTCACCGTTGTATc                                                                                                        <  2:468179/251‑1 (MQ=255)
                                                                                                                  ggCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAg                                                                                                                                     >  1:154725/1‑220 (MQ=255)
                                                                                                                  ggCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAg                                                                                                                                     <  2:154725/220‑1 (MQ=255)
                                                                                                                          cattggctgaatcGTATGAACCGATGGAGGTTAACGAATTCGGGGAAATCGTGCTTCTTGCAATGTGTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGAGTCTGAGCACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGTGGATAAAc                                                                                              <  2:414946/239‑1 (MQ=255)
                                                                                                                                        tATGAACCGATTGATGTTAACGAATTAGGTGAAATCGCTCTGCTTGCAATGGTTGAAGATGACATCATCCTCGCCTTGCCGGTAGTTTCGGAGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGACTTTGGTGAACTGCATGAAGAAGCGCAAAAACCAAACCCATTTGCCGTATTAGCCAGCTTAAAGAGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATTAACTGTAATTGAGGAGt                                                                                >  2:331764/1‑251 (MQ=255)
                                                                                                                                               cGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTcc                                                                         <  1:237042/251‑1 (MQ=255)
                                                                                                                                                   tGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATTGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCGCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTCCATg                                                                      <  2:265020/250‑1 (MQ=255)
                                                                                                                                                          aatgAATTCTGTGAAAGCGCTCTGCGTGCCATGGTTTAATGTGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATTATTCTGAACACTGTGAAGTGTCCGAAGCGGACATTGTCTTTAGTGAACTGCCTGAAGAAGCGCAAAAGCCATGACCAGTTGCCGTATTAGCCAGCGTAAAGCGGAAGTAATTGGTGCTCCCCGTTGGATCGGGGATTAACCGTAATTGAGGTGTAAGGTCCATGGCCGTaca                                                              <  2:400415/248‑1 (MQ=255)
                                                                                                                                                                    gTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTCCATGGCCGTACAACAGAATaa                                                     >  1:301535/1‑250 (MQ=255)
                                                                                                                                                                      gAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCTGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTCCATGGCCGTACAACAGAATAaac                                                   <  1:385698/250‑1 (MQ=255)
                                                                                                                                                                       aaaTCGATCTTCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGGGTCCGAAGCGGACATGGTCTTTGGTTAACTGCCTGAAGAAGCGCAAAAGCCAAACCCACTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGATGGATCGGGGATAAACCGTAATTGAGGAGTAAGTTCCATGGCCGTACAACAGAATAaacca                                                 <  1:362583/251‑1 (MQ=255)
                                                                                                                                                                           cGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGCAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTCCATGGCCGTACAACAGAATAAACCAAccc                                             >  2:409294/1‑251 (MQ=255)
                                                                                                                                                                              tCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGCGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTcc                                                                         <  1:252668/220‑1 (MQ=255)
                                                                                                                                                                              tCTGCTTGAAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGCGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTcc                                                                         >  2:252668/1‑220 (MQ=255)
                                                                                                                                                                                 tcTTGCAATGGTTGAAGATGAAATCAGCCTCGCCTTGCCTGTTGTTTCTGTTCATGATGCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTCCATGGCCGTACAACAGAATAAACCAACCCGtt                                          <  2:301535/247‑1 (MQ=255)
                                                                                                                                                                                 gCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTCCATGGCCGTACAACAGAATAAACCAACCCGTTc                                         <  2:35625/249‑1 (MQ=255)
                                                                                                                                                                                                aaGATGAAATCATCCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTACATGGCCGTACAACAGAATAAAACAACCCGATCCAAACGTGGCaagcat                         >  2:454971/1‑246 (MQ=255)
                                                                                                                                                                                                          catcCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTCCATGGCCGTACAACAGAATAAACCAACCCGTTCCAAACGTGGCATGCGTCGTTCCCATg               <  1:407182/251‑1 (MQ=255)
                                                                                                                                                                                                           atcCTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTCCATGGCCGTACAACAGAATAAACCAACCCGTTCCAAACGTGGCATGCGCCGTTCCCATg               >  2:195913/1‑250 (MQ=255)
                                                                                                                                                                                                             ccTCGCCTTGCCGGTAGTTCCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTCCATGGCCGTACAACAGAATAAACCAACCCGTTCCAAACGCGGCATGCGTCGTTCCCATGAcg            >  1:456901/1‑251 (MQ=255)
                                                                                                                                                                                                              cTCGCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTCCATGGCCGTACAACAGAATAAACCAACCCGTTCCAAACGTGGCATGCGTCGTTCCCATGAcgc           >  2:111996/1‑251 (MQ=255)
                                                                                                                                                                                                                 gCCTTGCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTCCATGGCCGTACAACAGAATAAACCAACCCGTTCCAAACGTGGCATGCGTCGTTCCCATGACGCGCt        >  1:473400/1‑251 (MQ=255)
                                                                                                                                                                                                                      gCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGTTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTCCATGGCCGTACAACAGAATAAACCAACCCGTTCCAAACGTGGCATGCGTCGTTCCCATg               <  2:210394/239‑1 (MQ=255)
                                                                                                                                                                                                                      gCCGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTCCATGGCCGTACAACAGAATAAACCAACCCGTTCCAAACGTGGCATGCGTCGTTCCCATg               >  1:210394/1‑239 (MQ=255)
                                                                                                                                                                                                                        cGGTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTCCATGGCCGTACAACAGAATAAACCAACCCGTTCCAAACGTGGCATGCGTCGTTCCCATGACGCGCTGACCGc  >  1:100104/1‑250 (MQ=255)
                                                                                                                                                                                                                         ggTAGTTTCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTCCATGGCCGTACAACAGAATAAACCAACCCGTTCCAAACGTGGCATGCGTCGTTCCCATGACGCGCTGACCGc  >  1:472096/1‑249 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 
GTTAAACGGCGATGCGAAGGTGACGGTAACGCTCGAGTGTCAGCGTTGCGGGAAGCCGTTTACTCATCAGGTCTACACAACGTATTGTTTTAGTCCTGTGCGTTCAGACGAACAGGCTGAAGCACTGCCGGAAGCGTATGAACCGATTGAGGTTAACGAATTCGGTGAAATCGATCTGCTTGCAATGGTTGAAGATGAAATCATCCTCGCCTTGCCGGTAGTTCCGGTGCATGATTCTGAACACTGTGAAGTGTCCGAAGCGGACATGGTCTTTGGTGAACTGCCTGAAGAAGCGCAAAAGCCAAACCCATTTGCCGTATTAGCCAGCTTAAAGCGTAAGTAATTGGTGCTCCCCGTTGGATCGGGGATAAACCGTAATTGAGGAGTAAGGTCCATGGCCGTACAACAGAATAAACCAACCCGTTCCAAACGTGGCATGCGTCGTTCCCATGACGCGCTGACCGC  >  NC_000913/1146973‑1147437

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