Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 302,186 C→T E130K (GAG→AAG)  paoA ← aldehyde dehydrogenase, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913302,1860CT96.5% 175.0 / ‑3.9 57E130K (GAG→AAG) paoAaldehyde dehydrogenase, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (23/32);  total (24/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

GACAGATGTTGCCGCTCATACGTTCACGGATCTCATCGGCAGTTGTTTCTGGAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCACCTTCAGTGTCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTG  >  NC_000913/301949‑302428
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
gACAGATGTTGCCGCTCATACGTTCACGGATCTCATCGGCAGTTGTTTCTGGAGCGGAAACCAAATCGCCCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAACGGCCGCCTGCATGGGGGGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCtgctgc                                                                                                                                                                                                                                       >  2:16941/1‑251 (MQ=255)
                    cGTTCACGGATCTCATCGGCAGTTGTTTCTGGAGCGGAAACCAAATAGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTcaaga                                                                                                                                                                                                                    >  1:135460/1‑247 (MQ=255)
                     gTTCCCGGATCTCATCGGCAGTTGTTTCTGGAGCGGAAACCCAATCTACCGCCTCCTGACTGGGAATGCCGTCCTGGATCTCTTTTAGCACCGCTACCGCTGCGCAAATTTGCCCGGAGGTTCAGATGCCGCACTTGAAGCCATCATGCTTGATAAAAGCCGCCTCCATGTGGTGCAGATTATCTGGCGAACCCAGGCCTTCACTGGTGGGGATCTTGGCCCCCTGCTGCATGACTGCacacgtccggccg                                                                                                                                                                                                                  >  2:178042/1‑239 (MQ=255)
                        cACGGATCTCATCGGCAGTTGTTTCTGGAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCaggcaggc                                                                                                                                                                                                                <  1:219041/250‑1 (MQ=255)
                         aCGGATCTCATCGGCAGTTGTTTCTGGAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCCATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCAt                                                                                                                                                                                                              >  2:225747/1‑251 (MQ=255)
                         aCGGATCTCATCGGCAGTTGTTTCTGGAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCAt                                                                                                                                                                                                              >  1:234805/1‑251 (MQ=255)
                                 cATCGGCAGTTGTTTCTGGAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTaa                                                                                                                                                                                                           <  1:225747/246‑1 (MQ=255)
                                 cATCGGCAGTTGTTTCTGGAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTg                                                                                                                                                                                                      >  1:261892/1‑251 (MQ=255)
                                    cGGCAGTTGTTTCTGGAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGGGGTGATCTTGGCCCCCTGATGCATGACTGAAAGCGTCAGGCAGGCATTAAGCCTGCgc                                                                                                                                                                                                   >  2:351666/1‑250 (MQ=255)
                                                  ggAGCGGAAACCAATTCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCTCCGCTACTGATGAGCAAATTTGCCCTGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAg                                                                                                                                                                                        <  2:350874/248‑1 (MQ=255)
                                                          aaccaagtcttccggtacgtgacTGGGAATGCCGTCCTGAATCTCTTGAAGCACCGCGACTGAGGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCATCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCCTGCATTAAGCCTGCGACCATTGACCAGTACGGTACAgg                                                                                                                                                                             <  2:321255/235‑1 (MQ=255)
                                                              aaaTCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTcc                                                                                                                                                                         >  1:35426/1‑251 (MQ=255)
                                                                 acgaccggtcggggacTGGGAATTCCGTCCTGAATCGCTTTTAGCACCGCTACTGATGAGCAAATTTGCCAGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCgta                                                                                                                                                                      <  2:12016/238‑3 (MQ=255)
                                                                    aCCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGAGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTg                                                                                                                                                                   >  1:234059/1‑251 (MQ=255)
                                                                           cgggacTGGGAATGCCGGCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGTTGCAGTAGCCGCACTGGAATCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCACCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGTACTGTCCGTgt                                                                                                                                                            <  2:224406/248‑2 (MQ=255)
                                                                                    gAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCt                                                                                                                                                   <  2:305911/251‑1 (MQ=255)
                                                                                        gCCGTCCTGAATCTCTTTTAGCACCGCGACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCAGTGACCAGCACGGTACAGTCTCCGCACTGTCCGTGATCGCAACCTTTc                                                                                                                                                <  1:118011/250‑1 (MQ=255)
                                                                                            ccctctagctctTTTAGCCCCGCTCCTGATGAGCAAGTTTCCCCGGTGTTGCCGGAGCCGCACTGGAAGCCATCATTCTTGATAAAGGCCCCCTGCCTGGGGTGAAGATTATCTGGCGATCCCCGGCCTTCAATTGTGGTCATCTTGGCCCCCTGATGCATGCCTGCAAGCGTCAGGCTTGCATTAAGCCTGCGACCATTGACCTGCTCGGTACTGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGt                                                                                                                                           <  2:35426/243‑1 (MQ=255)
                                                                                                gtaTCTCTTTTAGCACCGCGACTGATGAGCAAATTTGCCCGGATGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTTAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTAc                                                                                                                                         <  2:146129/247‑1 (MQ=255)
                                                                                                 aaTCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGa                                                                                                                                      <  1:224978/251‑1 (MQ=255)
                                                                                                        tttAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAAt                                                                                                                                <  1:177562/250‑1 (MQ=255)
                                                                                                        tttAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATg                                                                                                                               >  2:223367/1‑251 (MQ=255)
                                                                                                             atccGCTACGGATGAGCAATTTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGGATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCATGCATTAAGCCTGCGACCAGTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGAt                                                                                                                          <  2:128129/249‑1 (MQ=255)
                                                                                                                    tctgatgatcaaatgtctacGGAGTTGCAGTTGCCGCTCTGGAAGCTAGCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGAATATCTGGCGAACCCAGGACTGCAATGGTGGTGATCTTGGCCTCCTGATGCATGACTGCAAGCGTGAGGTAGGCATTAAGCCTGCGACCATTGTCCAGCTCGATACAGGCTTCGCACTATCCGTGATCGCATCCTTGGTTGGTACAGAGCAAAGGCAGAtttaca                                                                                                                     <  2:299044/230‑4 (MQ=255)
                                                                                                                            gCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCTGCACTGGCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAgtgt                                                                                                           <  1:138298/251‑1 (MQ=255)
                                                                                                                            gCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAgtgt                                                                                                           >  1:84066/1‑251 (MQ=255)
                                                                                                                             cAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTc                                                                                                          >  1:257639/1‑251 (MQ=255)
                                                                                                                                aTTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAg                                                                                                       <  2:257639/251‑1 (MQ=255)
                                                                                                                                    gCCCGGAGGTGCAGTAGCCGCCCTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTTAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGGCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGc                                                                                                   <  2:206329/251‑1 (MQ=255)
                                                                                                                                                    gCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTAGCTGGCGAACCCAGGACTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGGCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCa                                                                                   <  1:323898/251‑1 (MQ=255)
                                                                                                                                                    gCCGCACTGGAAGCCAACATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGACTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATGGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCa                                                                                   <  1:323916/251‑1 (MQ=255)
                                                                                                                                                       gCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAgg                                                                                                                                                                             <  1:40013/158‑1 (MQ=255)
                                                                                                                                                       gCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAgg                                                                                                                                                                             >  2:40013/1‑158 (MQ=255)
                                                                                                                                                              aaGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAgctg                                                                         <  1:345867/251‑1 (MQ=255)
                                                                                                                                                               aGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAgctgc                                                                        >  2:248292/1‑251 (MQ=255)
                                                                                                                                                               aGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAgctgc                                                                        >  2:285138/1‑251 (MQ=255)
                                                                                                                                                                  catcatGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCg                                                                     <  2:261892/251‑1 (MQ=255)
                                                                                                                                                                   atcatGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCgg                                                                    >  1:97973/1‑251 (MQ=255)
                                                                                                                                                                                  ggCCGCCTGACTGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATACCTACAAGCGCCAGGCAGGCATTAAGCCTGCGACCATTGACCCGCACGGTACAGGCTCCGCACTGTCCGTGATCGCACCCTTTCTTGGTACCGATCACATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCACCt                                                     >  2:39971/1‑251 (MQ=255)
                                                                                                                                                                                   gccgccTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCACCtt                                                    <  1:223367/251‑1 (MQ=255)
                                                                                                                                                                                         ggCATGGGTTGAAGATGATCGGGCGAACCCAGGCCGTCAATGGTGGTGATCTTGGCCCCCGGCTGCATGACAGCAAGCGGCAGGCAAGCATTAAGCCGGCTACCATTGATCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCTAGGAGCGTGGTTCGGGTATCCACCTCAAGCTTCTCGGTTTTGCCGTTCACCTTCAgtgt                                              <  2:349984/250‑1 (MQ=255)
                                                                                                                                                                                           cATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGCGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCACCTTCAGTGTc                                             >  2:107666/1‑250 (MQ=255)
                                                                                                                                                                                                   gAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCACCTTCAGTGTCAGGGGCAtt                                    <  2:196344/251‑1 (MQ=255)
                                                                                                                                                                                                    aaGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGTCGTTCCCCTTCAGTGTCAGGGGCAttt                                   >  2:234772/1‑250 (MQ=255)
                                                                                                                                                                                                    aaGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCACCTTCAGAGTCAGGGGCATTa                                   <  1:248292/251‑1 (MQ=255)
                                                                                                                                                                                                     aGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCCCCTTCAGTGTCAGGGGCATTAt                                  >  2:180806/1‑251 (MQ=255)
                                                                                                                                                                                                              ggCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCACCTTCAGTGTCAGGGGCATTATCTCTGgcgc                         >  2:223815/1‑251 (MQ=255)
                                                                                                                                                                                                              ggCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCACCTTCAGTGTCAGGGGCATTATCTCTGgcg                          <  1:180806/250‑1 (MQ=255)
                                                                                                                                                                                                                  aaCCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCACCTTCAGTGTCAGGGGCATTATCTCTGGCGCGGgt                     <  1:86192/251‑1 (MQ=255)
                                                                                                                                                                                                                           cTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGt                                                                                                                                           <  1:4415/124‑1 (MQ=255)
                                                                                                                                                                                                                           cTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGt                                                                                                                                           >  2:4415/1‑124 (MQ=255)
                                                                                                                                                                                                                            ttCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCACCTTCAGTGTCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCg           >  1:183181/1‑251 (MQ=255)
                                                                                                                                                                                                                                aTGGTTGTGAGCTTGGCCCCCGGATGCATGGCTGCAAGCGCTAGGCAGGCATTAAGCCTGCGCACATTGCTCAGCACGGTACATGCTCTGACCTGTCCGGGATCGCAACCTTTCTTGGTACCGATCAAATGCAGAGTTGCACGCAAAGTGTCCAGTAGTGTGGTTCGGGTATCCGCCCCAAGCTGCTCGGTGTTGTCGGTCGCCTGAAGGGTCACGGGCCTTACCTCTGCCGCGGGTGTAGTTGCCGGAAc       <  2:179988/251‑1 (MQ=255)
                                                                                                                                                                                                                                aTGGTGGTGATCTTTGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCAGTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGAACGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCACCTTCAGTGTCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGatt       <  2:27165/251‑3 (MQ=255)
                                                                                                                                                                                                                                aTGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCGGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCGCGGTTTTGCCGTTCACCTTCAGTGTCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAAc       <  2:90347/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gtggtgATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCACCTTCAGTGTCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCt    >  2:38328/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ggtgATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCACCTTCAGTGTCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTg  <  1:90247/251‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
GACAGATGTTGCCGCTCATACGTTCACGGATCTCATCGGCAGTTGTTTCTGGAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCACCTTCAGTGTCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTG  >  NC_000913/301949‑302428

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: