Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 885,118 A→G V38A (GTT→GCT)  ybjH ← uncharacterized protein YbjH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913885,1180AG98.2% 192.2 / NA 57V38A (GTT→GCT) ybjHuncharacterized protein YbjH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (30/25);  minor base T (0/1);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

AACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAACCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGCCGGGCGTGCTGCCGGTTTTATTGGTCAAATGGCGCTTCGTGTTTAAGAACTTTATCGATCACA  >  NC_000913/884877‑885361
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
aaCAGATGGGATATTCTCCCGCAGGGTAAAAAACTGCCTGAATGCTTTGTGCGATCAGGCATTCTCGAATTAATGGAGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGAGGAAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGGGGAACGCTGACTGAATCAGGagccacaaccc                                                                                                                                                                                                                                            >  2:283837/1‑250 (MQ=255)
  cAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACt                                                                                                                                                                                                                                          >  1:114571/1‑251 (MQ=255)
               ctcACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCACTCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCg                                                                                                                                                                                                                             >  1:210663/1‑251 (MQ=255)
                    gAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGcac                                                                                                                                                                                                                         >  1:5124/1‑250 (MQ=255)
                      aGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCCCATc                                                                                                                                                                                                                      >  2:314155/1‑251 (MQ=255)
                        ggTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTAGCCTTAATGGTGATGGTCGACAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGAAAGATGGCGAACATAAGGGTGCGGGAGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGAGAAGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCAACAACCACTTTGTATCCACGCGGCACATCg                                                                                                                                                                                                                     >  1:419824/1‑250 (MQ=255)
                         gTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGCACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGAc                                                                                                                                                                                                                   >  2:134081/1‑251 (MQ=255)
                             aaaaaTGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAAGTAATGGTGATGGTCGTCAAGCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTg                                                                                                                                                                                                               <  2:114284/251‑1 (MQ=255)
                              aaaaTGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCACTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACAACTTTGTATCCACTCGGCACATCGACGGTgc                                                                                                                                                                                                              >  2:383373/1‑250 (MQ=255)
                               aaaTGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGac                                                                                                                                                                                                             >  1:40109/1‑251 (MQ=255)
                                    ccTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTGGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACATCCTCTTTGTATCCACTCGACACATCGACGGTGACACttt                                                                                                                                                                                                        <  1:154812/251‑1 (MQ=255)
                                         cTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGcc                                                                                                                                                                                                   >  1:16036/1‑251 (MQ=255)
                                                    gATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTGTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCCCAcc                                                                                                                                                                                        >  2:298650/1‑251 (MQ=255)
                                                       cAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGt                                                                                                                                                                                     <  1:314348/251‑1 (MQ=255)
                                                        cggCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGCGACACTTTGCGCCATCGCCACACCAGt                                                                                                                                                                                     >  2:236433/2‑250 (MQ=255)
                                                        aGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGt                                                                                                                                                                                     >  2:321297/1‑250 (MQ=255)
                                                        aGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGt                                                                                                                                                                                     >  2:172180/1‑250 (MQ=255)
                                                           cATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTaa                                                                                                                                                                                   >  1:429180/1‑249 (MQ=255)
                                                            aTTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGc                                                                                                                                                                                >  1:83717/1‑251 (MQ=255)
                                                              tctcGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCACCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGcc                                                                                                                                                                               <  2:47662/250‑1 (MQ=255)
                                                                  gAATTACTGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCAtt                                                                                                                                                                           >  2:135804/1‑250 (MQ=255)
                                                                   aaTTAATGGTGATGGTCGGCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGTGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACGTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTaa                                                                                                                                                                         <  2:330284/251‑1 (MQ=255)
                                                                     ttAATGGCGATGGTCGTCACTATGGTGTTCGATAACCATACCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGAAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGCCACTTTGCGCCATCGCCACACCAGTAAAGCCCATTaaa                                                                                                                                                                        >  1:124741/1‑250 (MQ=255)
                                                                             gATGGTAGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGATGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCCTCGCAACACCAGTAAAGCCCATTAAAAGCGCGccc                                                                                                                                                               >  2:373120/1‑251 (MQ=255)
                                                                                ggTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAAc                                                                                                                                                            <  1:134081/251‑1 (MQ=255)
                                                                                   cgtcAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGCCCGCTTGTGGAACGCTGACTGAATAAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCCTCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGt                                                                                                                                                         >  2:91055/1‑251 (MQ=255)
                                                                                    ggcAATCGGGTGTTCGATAACCATCCCGTCTCCTACGCTGGCAAGATGGCGAACAGAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGGCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCTCAACCACTTTGTATCCGCTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTa                                                                                                                                                        <  1:374987/249‑1 (MQ=255)
                                                                                       aaTCTGGTGTTCGATAACCATCCCTTCACCTACGCGGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGAc                                                                                                                                                     <  2:255393/251‑1 (MQ=255)
                                                                                          cTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAAt                                                                                                                                                  >  1:22189/1‑251 (MQ=255)
                                                                                            ggTGTTCGATAACCATCCCTTCCCCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACAGAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCGTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACTGTGACACTTTGCGCCATCGTCACACCAGAAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAAttt                                                                                                                                                <  1:381845/251‑1 (MQ=255)
                                                                                            ggTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACAGAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACAGAGACGGTGTGCGGCACAGTCGCGTCGCTGTCCGCTGGTGTAACGCTGACTGAATCAGGATCCACAACCACTTGGTATCCACTCGGCACATCGACGGTGACACTGTGCGCCATCGCCACACCAGTAAAGCCCAGGAAAAGCGCGCCCAGCAGTAGACAAttt                                                                                                                                                <  2:225812/251‑1 (MQ=255)
                                                                                              ggtTCGATAACCATCCCTTCACCTACGCTGGCAAGAGGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACTCCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAAttttt                                                                                                                                              <  2:5124/250‑1 (MQ=255)
                                                                                                    aTAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAAt                                                                                                                                        <  1:167296/251‑1 (MQ=255)
                                                                                                      aaCCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCa                                                                                                                                      <  1:236433/251‑1 (MQ=255)
                                                                                                       aCCATCCCTTCACCTACGCTGGCAACATGGCGAACATCAGGATGAGAGCAGTAAGCCGGAGCTGGCGCAGGACACAAATAGACGGTCTGCGGCACAGTCGCAACGCCGACCGCTTGTGGAACACTCGCTGAACCAGGAGCCACACCCACTTTGAATCAACTCGGCACATCGACGGTGACACTTTCCGCCATCGACACACCAGTAACGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCAt                                                                                                                                     >  1:406982/1‑251 (MQ=255)
                                                                                                        ccATCCCTTCACCTACGCTGGGAAGAGGGGGAACATAGGGATGCGGGCGGTAAGCCGGATCTGGCGCAGTGGCCACATAGACGGTGTGCGGCACAGACGCGACGCTGCCCGCTTGTGGAACGCTGACTGAATCAGGAGCCACATCGACTTTGTATCCACTCGGCACATCTACGGTGACACTTTGCGCCATCGCCACACCAGGAAAGCCCATGAAAAGCGCGCCCAACAGTAGACAATTGTTCATAATCAtt                                                                                                                                    <  2:24069/251‑1 (MQ=255)
                                                                                                             ccTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGcc                                                                                                                               <  1:338151/251‑1 (MQ=255)
                                                                                                                                             aGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAAt                                                                                               >  1:408012/1‑251 (MQ=255)
                                                                                                                                                    gggcggTAAGCCGGAGCTGGCTCAGGAGCCACATAGACGTTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACGCCATTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTTAAATAATGCAGCGGGTTAAGTTCCGCGAATTATGCAg                                                                                        <  1:389697/251‑1 (MQ=255)
                                                                                                                                                               cGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCCCAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAg                                                                             <  1:82512/251‑1 (MQ=255)
                                                                                                                                                                           aGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGTATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCAAACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCTTTCATTCAGGCTTTGTGTTGc                                                                 >  2:19139/1‑251 (MQ=255)
                                                                                                                                                                            ggggCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGcc                                                                <  2:83717/248‑1 (MQ=255)
                                                                                                                                                                             gAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTGTCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGCCg                                                               <  1:52760/251‑1 (MQ=255)
                                                                                                                                                                              aGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGCCgg                                                              >  2:139031/1‑251 (MQ=255)
                                                                                                                                                                                     aaGACGGTTTGTGGCCCAGTCGCGACGCTAACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCCCTTTGTAACCACTCGGCACATCGACGGTGACTCTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACATTAGACAATTTTTCATAATCATTGCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGCCGGGCGtgct                                                       <  2:114571/250‑1 (MQ=255)
                                                                                                                                                                                      aGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGACCGCTGACTGAATCAGGAGCCACAACACCTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATAGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGCCGGGCGtgctg                                                      >  1:422103/1‑251 (MQ=255)
                                                                                                                                                                                        aCGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGCCGGGCGTGCTGcc                                                    <  1:283837/251‑1 (MQ=255)
                                                                                                                                                                                              tGCGGCACAGTCGCTGCGCTGCCCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGCCGGTCGTGCTGCCGGttt                                               <  2:408012/250‑1 (MQ=255)
                                                                                                                                                                                                cGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGCCGGGCGTGCTGCCGGTTTTAt                                            >  2:79142/1‑251 (MQ=255)
                                                                                                                                                                                                     caGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTGGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGCCGGGCGTGCTGCCGGTTTTATTGGTc                                       <  1:298650/251‑1 (MQ=255)
                                                                                                                                                                                                       gTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGCCGGGCGTGCTGCCGGTTTTATTGGTCa                                      >  1:349307/1‑250 (MQ=255)
                                                                                                                                                                                                               gCTGACCGCTTGTGGAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGCCGGGCGTGCTGCCGGTTTTATTGGTCAAATGGCGCt                             >  2:180652/1‑251 (MQ=255)
                                                                                                                                                                                                                            ggAACGCTGACTGAATCAGGAGCCACAACAACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGCCGGGCGTGCTGCCGGTTTTATTGGTCAAATGGCGCTTCGTGTTTAAGa                 >  1:352976/1‑250 (MQ=255)
                                                                                                                                                                                                                             gAACGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGCCGGGCGTGCTGCCGGTTTTATTGGTCAAATGGCGCTTCGTGTTTAAGAAc               <  1:373120/251‑1 (MQ=255)
                                                                                                                                                                                                                               aCGCTGACTGAATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGGAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTGCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTGCATGCAGGCATTGTGTTGCCGGGCGGGCTGCCGGGTTTATTGGTCAAATTGCGCTTCGTGTTTAAGAACtt             <  2:210663/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gactgaATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGCCGGGCGTGCTGCCGGTTTTATTGGTCAAATGGCGCTTCGTGTTTAAGAACTTTAt          <  1:371179/249‑1 (MQ=255)
                                                                                                                                                                                                                                    gactgaATCAGGAGCCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGCCGGGCGTGCTGCCGGTTTTATTGGTCAAATGGCGCTTCGTGTTTAAGAACTTTATCg        >  2:42809/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aTCAGGACCAACAACCACTTTGTATCCACTCGGCACATCCACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGCAGACAATTTTTCATAATCATTTCGCCTGAAATACTGCAGCGGGTTAATTTCCGCGAATTATGCCGCCTTCATTCAGGCATTGTGTTGCCGGGCGTGCTGCCGGTTTTATTGGTCAAATGGCGCTTCGCGTTTAATGACTTTATCGTTcacc  >  2:377333/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
AACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACGCTGACCGCTTGTGGAACGCTGACTGAATCAGGAACCACAACCACTTTGTATCCACTCGGCACATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGCCGGGCGTGCTGCCGGTTTTATTGGTCAAATGGCGCTTCGTGTTTAAGAACTTTATCGATCACA  >  NC_000913/884877‑885361

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: