Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,582,272 C→T G169E (GGG→GAG)  ypfG ← DUF1176 domain‑containing protein YpfG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,582,2720CT94.2% 173.7 / ‑1.8 71G169E (GGG→GAG) ypfGDUF1176 domain‑containing protein YpfG
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (29/36);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CAAATCAATGGTGTTATAGGCCCCTGCCTCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCCCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCATCGGTTACTAATAACCATGGTGAAATCCGCCACTTGTCACCACTTAACGCCAACGGCTCGCCATCTAACAGCAGCCGT  >  NC_000913/2582025‑2582517
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          
cAAATCAATGGTGTTATAGGCCCCTGCCTCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCtt                                                                                                                                                                                                                                                         >  1:495211/1‑251 (MQ=255)
   accaaaggTGTTATAGGCCCCTGCCTCACAGCTAACCATCATCAGCGCTTTATCATCAGCCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCCGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCTGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCt                                                                                                                                                                                                                                                      <  2:418992/245‑1 (MQ=255)
       aTGGTGTTATAGGCCCCTGCCTCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGTCGTCGGCGTTGGGTTACCCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCATCCCCTTTCtttat                                                                                                                                                                                                                                                  >  2:50147/1‑248 (MQ=255)
         gttgtTATAGGCCCCTGCATCACAGCTAATAATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCACGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCCTAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGGACTCTGAGAGGCGGTTCGGCCACTTTTTTTATcc                                                                                                                                                                                                                                                >  2:415228/3‑251 (MQ=255)
          gtgtTATAGGCCCCTGCCTCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCa                                                                                                                                                                                                                                               <  1:446455/251‑1 (MQ=255)
           tgtTATAGGCCCCTGCCTCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTGGTCCTCTTTCTTTATcccc                                                                                                                                                                                                                                              >  2:93430/1‑249 (MQ=255)
           tgtTATAGGCCCCTGCCTCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCtctttctttgtccac                                                                                                                                                                                                                                              >  1:345335/1‑244 (MQ=255)
            tttATAGGCCCCTGCCTCACAGCTAATCAGCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCAc                                                                                                                                                                                                                                              <  1:574894/249‑1 (MQ=255)
            gtTATAGGCCCCTGCCTCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCAc                                                                                                                                                                                                                                              <  1:414140/250‑1 (MQ=255)
              tataGGCCCCTGCTCCACAGCTAATCAACACCAGCACTTTATCATCAGTCAGAGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGTGCAGCGCAGCCCATTCATCCGCCAGTTTCCATAATCCAGCAAACCGTGGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGTCCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTGCTTGATCCAcgc                                                                                                                                                                                                                                            <  2:563320/250‑1 (MQ=255)
                   gCCCCTGCCTCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCATCCTCCTTCTTTTTCCACCCGGTTTc                                                                                                                                                                                                                                      >  1:370629/1‑251 (MQ=255)
                         gCCTCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTg                                                                                                                                                                                                                                  <  1:562711/249‑1 (MQ=255)
                            ccacaGCTAACCATCATCAGCGCTTTAGCATCAGGCAGCGCAGTGACATTCACCTCACGACGCAATGGACCAAGCGGGCAGCGCAGACCATGCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTGCACTGCCAAc                                                                                                                                                                                                                             <  2:309382/250‑1 (MQ=255)
                             cacaGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCCCCGCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGAGGGGTCACTGACAAcg                                                                                                                                                                                                                            <  1:586977/251‑1 (MQ=255)
                                 gCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGGTTCACTGCCAACGCGCt                                                                                                                                                                                                                        >  1:641484/1‑251 (MQ=255)
                                   tAATCATCATCAGCGCTTTATCATCAGTAAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTGCGTCCTCTTTTTTTATCCAAGCGGGTTCCCTGCCAAccctcttt                                                                                                                                                                                                                      >  2:556948/1‑244 (MQ=255)
                                    aaTCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCAGTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCtttt                                                                                                                                                                                                                     <  2:275019/251‑1 (MQ=255)
                                     atcatcatcAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTTCCAACGCGCttttt                                                                                                                                                                                                                    >  1:195264/1‑250 (MQ=255)
                                           atcaGCGCTTTATCACCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGc                                                                                                                                                                                                              <  1:354876/251‑1 (MQ=255)
                                                                        gACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGc                                                                                                                                                                               <  1:50147/251‑1 (MQ=255)
                                                                        gACAGTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGc                                                                                                                                                                               <  1:320511/251‑1 (MQ=255)
                                                                             tCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATACTCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTAGTCCTCTTTCTTGATCCACGCGGTTTCACCGCCAGCGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTc                                                                                                                                                                          >  1:326023/1‑251 (MQ=255)
                                                                             tCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACt                                                                                                                                                                           <  1:643339/250‑1 (MQ=255)
                                                                                    acgacgCAATGGAACAAGCGAGCAGCGCAGTCCATTCTTCCGCCAGTTTCCATAATCCAGCAAATCGTGGCGTTCTTCGAGTGAGAGTGGCGGCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCCCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTGTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGaaa                                                                                                                                                                   <  2:351352/251‑1 (MQ=255)
                                                                                       acgCAATGGATAAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCAATAATCCAGCAAATCGTTGCGTTCTTCGAGTAAGAGTGACGTAGGCGTTTGGTTAACCACCGCGCCCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCTGTTCGTCCTCCTTCTTGTTCCCCCCGGTGTCACTGCCAACGCGCTTTTGCTCAGCATCAAATAAACACAACACCGCTTTTaaaacaccaatgcaaaagg                                                                                                                                                                >  2:68812/1‑234 (MQ=255)
                                                                                        cgCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACATCAACGCTGCTTTTAAGCCACTCAGAGAAATGGt                                                                                                                                                               <  1:167938/251‑1 (MQ=255)
                                                                                            aTGGATCAAGCGAGCAGGGCAGACCATTCATCCGCCAGGTTCCATAATCCAGCAAAACGTTGCGTTCTCCGAGTGAGAGTGGCGTCGGCTTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGATCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCCCTCAGAGAAATTGTCTGa                                                                                                                                                           <  2:555852/251‑1 (MQ=255)
                                                                                                            gagcggccctttCATCCGCCAGTTCCCAAAACCAGCCAAATCGTGCCGTTTTACGAGTGAGAGTGGGGTGGGCGTTGGGTTAACCACCGCGCCCTCTTCCGGCGCAGGCGCGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCGTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGt                                                                                                                                         <  1:53975/241‑1 (MQ=255)
                                                                                                              gcAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTg                                                                                                                                        <  2:666472/251‑1 (MQ=255)
                                                                                                              gcAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGATCTTCGAGTGAGAATGGCATCGGCGTTATGATAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCCTCATCTTTGTTGATCCACCAGGTTTCACGGCCAGCGCACTTATGGACAGAATCAATACACGACCACACTGCTTTTTAGCCACCCAGAGCAATGGtattatacccctacccaaaggtg                                                                                                                                         >  2:350460/1‑229 (MQ=255)
                                                                                                                aGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTTATCCCCGCGGGTTCACTGCCAACGCGCTTTTGCTGTGCCTCAATAAACAACCACGCCGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGAt                                                                                                                                      >  1:595389/1‑251 (MQ=255)
                                                                                                                           ccGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAAAAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGccttccctt                                                                                                                           >  1:173876/1‑251 (MQ=255)
                                                                                                                              ccaggtttcATAATTCAGCAAATCGTTGAGTTCTTCGAGTGAGAGGGGCGTCGGCGTTGGGTTAACCACCGCGACTTTTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCTGTGGTTCGTCCTTTTTCTTGATCCCCTCGGTTTCACTGTCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAATCCACTCAGAGAAATGGTCTGATCGTCATCGCGTAAGGTGATTGCCTTCCCTTCCt                                                                                                                        <  1:430082/243‑1 (MQ=255)
                                                                                                                              ccAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCAACGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCCGTTTTACTGCCCACGCGCTTTTGCTTGGCATCAATTAACAACAACGCTACTTTTAAGCCACTCAGAGGAATGGGCTGATCGCCATCGCGTAGGGTGATTGccttacctttct                                                                                                                        >  2:185115/1‑248 (MQ=255)
                                                                                                                                 gTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCCTCCTCTTTCTTGAGCCACGCGGGTTCACTGCCAACGCGCCTTTGCTGCGCATCCATAAAACAACACCCTGCTTTTCAGACACTCCGAGAAATGGTCTTATCGCCATCGCGGAAGGGGACTGCCTTCCCTTCCttg                                                                                                                      >  2:460010/1‑248 (MQ=255)
                                                                                                                                     ccATAATCCAGCAAATCTTTGCTGTCTTCGAGTGAGAGTGGCGTCGGCGGTGGGTTAACCCCCGCTACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGGTCGTCCTCTTTCGTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAatcat                                                                                                                 <  2:495211/251‑1 (MQ=255)
                                                                                                                                     ccATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGGGATTGCCTTCCCTTCCTGAatcat                                                                                                                 >  1:581759/1‑251 (MQ=255)
                                                                                                                                         aaTCCAGCAAATCGGTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTGc                                                                                                             <  1:350460/251‑1 (MQ=255)
                                                                                                                                           tCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCaa                                                                                                                                                                                                                              >  1:639109/1‑139 (MQ=255)
                                                                                                                                           gccAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCaa                                                                                                                                                                                                                              <  2:639109/138‑1 (MQ=255)
                                                                                                                                                  aaaTCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCCACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCACAAAcgc                                                                                                    >  1:350957/1‑251 (MQ=255)
                                                                                                                                                   aaTCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTACCCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGGTTCACTGCCAACGCGCTTTTGCTTGGGATCAATAAACAAAAACGCTGGTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGGTATTGCCTTTCCTTCCTGAATCATCTGCaaaaaag                                                                                                     >  1:267418/1‑247 (MQ=255)
                                                                                                                                                        ttGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCGGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACgct                                                                                                                                                                                         <  1:655608/162‑1 (MQ=255)
                                                                                                                                                           cGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGAGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGAACGATGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGGTTCAATGCCAAAGCGCTTTTGGTGAGCATCAATAAACAACAAGGCGGCATTTAAGCCACACAGAGGAAAGGTCTTATCGCCATCGCGGAAGGTGATTGCCTTCCCTTCCTGAATAATCTGAAACAAAGCGGttattat                                                                                           >  2:16117/1‑245 (MQ=255)
                                                                                                                                                                  cGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAAAGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAATGTGATTGCCTTCCATTCCTGAATACTCTGCAAAAACGCGGTAAGGGGTGCCGt                                                                                     >  1:494705/1‑250 (MQ=255)
                                                                                                                                                                   gAGTTAGAGTGGCGTCGGCGTTAGGTTAACCACCGCGAACTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGATCCGCCACCTTCTTGATACCCACCGTTACACTGCCCACGCGACTTTGGTGAGCCTCAATAAACAACAACACTGCTTTAAAACCACACCGAGGAATGGTCTGATCGCCTTCCAGTTAGGTGGTTGACTGTCCTACCTGAATCATCGTCCAAAACGCggggactgttgcccaac                                                                                   >  2:453138/1‑236 (MQ=255)
                                                                                                                                                                     gTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCCCTTTCTTGATCCACGCGGGTTCCCTGCCCACCCGCTTTTGCTGAGCATCAATAAAACACACCGCTGCTTTTAAGCCACCCCGAGCAAAAGTCTGAGCGCCATAGTGTAAGGTGATTCCCTTCCCTTCCTGAATCCTCTACAAAAACGCGGTGATTGTGTCCGTatcc                                                                                 >  2:505058/1‑250 (MQ=255)
                                                                                                                                                                     gTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGAGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGGTTCAATTCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGGCTGATCGCCATCGCGTAAGGGGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTatc                                                                                  >  1:569998/1‑250 (MQ=255)
                                                                                                                                                                     gTGAGAGGGGGGGCGGTGTTGGGTTTACCCCCGCGGCCTCTTTCAGCGCAGGCGCGGGTGGTACGCGGAGCGGCGGGTTGTACTCTTTCTTGATCCACGCGGTGTCACTGCAAACGCTCGTTTGCTGAGCATCAATTAACAACCAAGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGAGCGCAATCGCGTAAGGTGATTGCCTGCCCTTCCTGAATCATCGGCAGAAACTCGGTGATGGTTGTCGTagca                                                                                 <  2:50076/251‑4 (MQ=255)
                                                                                                                                                                       gagaGTGGCGTCGGCGTTGGGTTAACCACCGCGCCCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACCCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACCCAGAGAAATGGTCTGATCGCCCTCGCGTAAGGAGATTGCCTTCCCTTCCTGAAACATCTGACACAACGCGGGGATGGTTGCCGTatcatc                                                                               >  1:492064/1‑251 (MQ=255)
                                                                                                                                                                       gagaGGGGCGGCGGCGTTGGGTTAACCCCCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTGAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCGGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTatcatc                                                                               <  2:350957/251‑1 (MQ=255)
                                                                                                                                                                          agTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCATCGGt                                                                            <  2:641484/251‑1 (MQ=255)
                                                                                                                                                                           gTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCAGCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCATCGGtt                                                                           <  1:96344/251‑1 (MQ=255)
                                                                                                                                                                              gCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCATCGGTTACt                                                                        <  1:415228/251‑1 (MQ=255)
                                                                                                                                                                                 tCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCCTCGCGTAAGGGGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCATCGGTTACtaat                                                                    >  1:165872/1‑251 (MQ=255)
                                                                                                                                                                                  cGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATTGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCGTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCATCGGTTACtaat                                                                    <  1:622454/250‑1 (MQ=255)
                                                                                                                                                                                         gggTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCATCGGTTACTAATAACCATg                                                             <  1:568918/250‑1 (MQ=255)
                                                                                                                                                                                            ttAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGGTCCTCCTCTTTCTTAGTCCACGCGGTTTCACTGCCCAAGCGCTTTTTCTGAGCATCAATAAACAACACCGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATAGCCATCGAGTAAAGAGATTGCCTTCCCTTCCTGGATCATCTGCCAAAAACCGGAGATGGTTGCCGTATCCTCGGTTACTAATAACCAtaggg                                                          >  2:174962/1‑246 (MQ=255)
                                                                                                                                                                                             tAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCCTCCTCTTTCCTGATCCCCGCGGGTTCACTGCCAACCCGCTTTTGCTGCGCATCAATAAACAACAACGCTGCTTTTAAGCCAATCAGAGAAATGGTCTGCTCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGGGATGGTTGCCGCAACATAGGTTACTAAAAACCATgggga                                                         >  2:79473/1‑247 (MQ=255)
                                                                                                                                                                                              aaCCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCATCGGTTACTAATAACCATGGTGaa                                                        <  1:92210/250‑1 (MQ=255)
                                                                                                                                                                                               accaccGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAATCGCGGTGATGGTTGCCGTATCATCGGTTACTAATAACCATGGTGAAAt                                                      <  1:213187/251‑1 (MQ=255)
                                                                                                                                                                                                       gACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCCCTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGCATCATCTGCAAAAACGCGGGGATGGTTGCCCTATCATCGGTTACTAATAACCATGGTGAAATCCGCCACt                                              >  1:414517/1‑251 (MQ=255)
                                                                                                                                                                                                       gACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAatcat                                                                                                                 <  1:95269/186‑1 (MQ=255)
                                                                                                                                                                                                       gACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGGTTCCCTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACCACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCCTCGCGTAAGGTGATTGCCTTCCCTTCCTGAatcat                                                                                                                 >  2:95269/1‑186 (MQ=255)
                                                                                                                                                                                                             tttCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTTTTGATCCACGCGGGTTCACTGCCCACGCGCCTTTGCTGAGCATCAATAAAAAACAACGCTGCTTTTTAGCCACCCAGAGAAATGGTCTGATACCCCACGCGTAAGGGGATTGCCTTCCCCTCCTGAATCATCTGCCACAACGCGGGGGTGGGTGCCGGATCCACGGTTACTTAATAACCATGTGAAAttccccccctctcac                                        >  2:253202/1‑237 (MQ=255)
                                                                                                                                                                                                                 aGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGGGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCATCGGTTACTAATAACCATGGTGAAATCCGCCACTTGTCACCACt                                    >  1:89483/1‑251 (MQ=255)
                                                                                                                                                                                                                           cgggcggTACGCTGGGCGGCGGTTAGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCTTCGCGGAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCATCGGTTACTAATAACCATGGTGAAATCCGCCACGTGTCACCACTTAACTCCAAc                          <  2:173876/251‑1 (MQ=255)
                                                                                                                                                                                                                                ggTACGCTGATGGGCGGGTCGGCCTTTTTCTTGATCCACGCGGTTTCCCTGGCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTTATGGTTGCCGTATCATCGGTTACTAATAACCATGGTGAAATCCGCCACTTGTCACCACTTAACGCCAACGGCTc                     <  1:216540/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCAGCGGTTACTAATAACCATGGTGAAATCCGCCACTTGTCACCACTTAACGCCAACGGCTCg                    <  1:556948/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gTACGCTGAGCGGCGGTTCGTCCTCTTTATTTATCCACTCGTTTGCACTGCAAACGCGATTTTGTGGAGAATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGGAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGTTGATGGTTGCCGTATCATCGGTTCCTAATAACCATGGTGAAATCCGCCACTTGTCACCACTTAACGCCACCGGCTCg                    <  2:492064/251‑1 (MQ=255)
                                                                                                                                                                                                                                             gcggTTCGTCCCCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAATGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCATCGGTTACTAATAACCATGGTGAAATCCGCCACTTGTCACCACTTAACGCCAACGGCTCGCCATCTAAcagc        >  1:26523/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   cGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCATCGGTTACTAATAACCATGGTGAAATCCGCCACTTGTCACCACTTAACGCTAACGGCTCACCATCTAAcagccaccgg  >  2:472069/1‑245 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          
CAAATCAATGGTGTTATAGGCCCCTGCCTCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCCCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCATCGGTTACTAATAACCATGGTGAAATCCGCCACTTGTCACCACTTAACGCCAACGGCTCGCCATCTAACAGCAGCCGT  >  NC_000913/2582025‑2582517

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