Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,366,433 G→A G295D (GGC→GAC)  puuE → 4‑aminobutyrate aminotransferase PuuE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,366,4330GA94.7% 100.5 / ‑3.8 38G295D (GGC→GAC) puuE4‑aminobutyrate aminotransferase PuuE
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (16/20);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

ATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGGCACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGTACGCGGCCTGGGGTCGATGATTGCGGTAGAGTTTAACGATCCGCAAACGGGCGAGCCGTCAGCGGCGATTGCA  >  NC_000913/1366219‑1366662
                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      
aTTCGCCGCCGGTGCGACGAGCCAGGTAGTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTGGCGCGTACCGGTAAGCTGTTTTCCATGGATCATGACGCCGATAAGCCGGATTTAATGTCGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGAAACCTACGCCGGTAACCCGCTGGCGGTGGCTGcc                                                                                                                                                                                                      <  1:221773/249‑1 (MQ=255)
    gccgccTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCAGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCAgcgccc                                                                                                                                                                                                >  2:135469/1‑249 (MQ=255)
    gccgccTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAAGGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCAc                                                                                                                                                                                                <  2:397558/251‑1 (MQ=255)
    gccgccTGTGCGACGAGCACGGTATTGGGATGATTGCTGATGAAGCGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGATGATGGCGAAAAGCCTCGCGGGCGGGATGCTGCTTTCGGGCGTGGTCGTTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCAc                                                                                                                                                                                                <  1:419/251‑1 (MQ=255)
                acgaGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGCTGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGAGGCTGCCCCGCACGCGGTGCTcca                                                                                                                                                                                     >  1:406342/1‑248 (MQ=255)
                 cgaGCACGGTATTGTGCTGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCTGTAAGCTGTTTTCCATGGATCATTACGCCGATAAGCCGGAGTTAATTACGAGGGCGAAAAGCCGAGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTAGGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACa                                                                                                                                                                                   <  2:102369/251‑1 (MQ=255)
                     cACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACCTTAt                                                                                                                                                                               >  1:231432/1‑251 (MQ=255)
                              gtgATTATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGGTTGCCATGGAGCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCGTTCTGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTGGGCGACACCTACGCCGGTAACCCGCTGGCGGGGGCTGCCGCGCACGCGGTGCTCAACATTACCGACAAAg                                                                                                                                                                       <  2:74766/250‑1 (MQ=255)
                                gatgatTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCCACATTAAAGGCAAAGAAt                                                                                                                                                                    >  2:175520/1‑251 (MQ=255)
                                           atgaAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAAAGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGTCGACACCTACGCCGGTAACCCGCTGGCGGGGGCTGCCGCGCAAGCGGTGCTCAACAGTATCGGCAAATAATCACTCTGCGa                                                                                                                                                          >  1:82565/1‑250 (MQ=255)
                                            tgacgTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCAATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGGAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGCCAAAGAATCACTCTGcccac                                                                                                                                                        >  2:3567/5‑247 (MQ=255)
                                              aaGTGCAAAGCGGCTGTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAAcgc                                                                                                                                                      <  2:89272/251‑1 (MQ=255)
                                                   cAAAGCGGCCTGGCTCATACCGTTAATCTGTTTTCCATGGATCATGCCGCCGATAAGCCGGACTTAATGACGATTGCGAAAAGCATCGCGGTCGGGATGCCGCTGTCGGGCGTGGTCGGGAACGCGAATATTATTGACGCACCCGCGCCGGGCGGGCTTTGCGACACCTGCGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCACCAGTAGCGACAAAGATTCACTCTGCGTACGCGCGaa                                                                                                                                                 <  2:413213/251‑1 (MQ=255)
                                                             ttGCGCGTACCGGTAAGCTGTTTGCCATTGATCATTACGCCGATAAGCCGGATTTAATGACGCTGGCGAAAAGCCGCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTgg                                                                                                                                         <  2:65362/249‑1 (MQ=255)
                                                                cgcgTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTAGCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAg                                                                                                                                    <  2:41170/251‑1 (MQ=255)
                                                                     aCCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGGGCTCAACATTATCGACa                                                                                                                                                                          >  1:59970/1‑208 (MQ=255)
                                                                     aCCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATAGTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACa                                                                                                                                                                          <  2:59970/208‑1 (MQ=255)
                                                                                gTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGTGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGt                                                                                                                   >  1:328306/1‑251 (MQ=255)
                                                                                           aTCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGCATATTATGGGCGCACCCGCGCCGGTCGGGCTTGGCGACACCTACGCCGGGAACCCCCTGGCGGTAGCGTCCGCGCACGCGGTGCTCCACATTCTCTGCCACGCATCACTCTGCCCACGCGCGAATCAACTGGGCACGCGGCTCAAAAACACGtttgtttatgcc                                                                                                        >  2:108624/1‑241 (MQ=255)
                                                                                               ttACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCGCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCaaa                                                                                                     <  2:182466/250‑1 (MQ=255)
                                                                                                tACGCCGATAAGCCTGATTTAAGGACGAGGGCGAAAAGCCTCGCGGGTGGGAGGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCGCAACATTAGCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCaaag                                                                                                    <  2:378581/250‑1 (MQ=255)
                                                                                                  cGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCaaagaaa                                                                                                 <  1:1443/251‑1 (MQ=255)
                                                                                                       aTAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCCCGCGGTGCTCAACATTATCGCCAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCCGCGTCTCAAAAACACGTTGATTGATGCCCCAGAAAGCGtt                                                                                            >  2:111470/1‑251 (MQ=255)
                                                                                                       aTAAGCCGGATTTAAGGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCGCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGtt                                                                                            <  1:402688/251‑1 (MQ=255)
                                                                                                         aaGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTTATGCCAAAGAAAGCGTTc                                                                                           <  1:108624/250‑1 (MQ=255)
                                                                                                         aaGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTGTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCTCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTcc                                                                                          <  1:4609/251‑1 (MQ=255)
                                                                                                                     aTGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCTAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATtgctgc                                                                              >  1:147934/1‑251 (MQ=255)
                                                                                                                                      ccTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGACTATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTGAACCGCTGGCGGTGGCTGCCGCGCACGCGGCGCTCCTCATTATCGACAAAGGGTCCCTCTGCGAACGCGCGGATCAACTGGGCCACCGTCTCAAAAACACGTTGATCGATGCCAAAGAAAGCGTTACGGCCCTTGCTGCGGTACGCGGCCTGGGGt                                                             >  1:187075/1‑251 (MQ=255)
                                                                                                                                         cgcgGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGTACGCGGCCTGGGGTCg                                                           <  2:98111/250‑1 (MQ=255)
                                                                                                                                                      ccGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGTACGCGGCCTGGGGTCGATGATTGCGGTaga                                             <  1:200241/251‑1 (MQ=255)
                                                                                                                                                                      tCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGGCACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGTGAACGCGCGAATCAACTGGGCCGGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGTACGCGGCCTGGGGTCGATGATTGCGGTAGAGTTTAACGATCCGCa                              >  1:165725/1‑250 (MQ=255)
                                                                                                                                                                      tCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGGCACCTACGCCGGTAACCCGCTGGCGGCGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACAAGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGTACGCGGCCTGGGGTCGATGATTGCGGGAGAGTTTAACGATCCGCa                              >  1:137629/1‑250 (MQ=255)
                                                                                                                                                                         gTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGTACGCGGCCTGGGGTCGGTGATTGGGGTAGAGGTTAACGATCCGCAAACg                          >  2:131663/1‑251 (MQ=255)
                                                                                                                                                                              gcgAATCTTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACCTTATCGACAAAGAATCACTCTGCGCACGCGCGAATCCACCGGGCCAGCGTCTCAAAAAACAGGTGATTAGTGCCAAAAAAAACGTTCCGGCCCTTGCTGCGGGAAGCGGCCCGGGGGCGAGGATTGCGGTAGAGTTTAAAGATCCGccacaggggg                      >  2:385780/1‑241 (MQ=255)
                                                                                                                                                                              gcgAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGTACGCGGCCTGGGGACGAGGATTGCGGTAGAGTTTAACGATCCGCAAACGGGCg                      >  1:292059/1‑250 (MQ=255)
                                                                                                                                                                                          gACGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGTACGCGGCCTGGGGTCGATGATTGCGGTAGAGTTTAACGATCCGCAAACGGGCGAGCCGTCAgcggc         <  1:220271/251‑1 (MQ=255)
                                                                                                                                                                                           aCGCACCCGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGTACGCGGCCTGGGGTCGATGATTGCGGTAGAGTTTAACGATCCGCAAACGGGCGAGCCGTCAgcggcg        <  1:111470/251‑1 (MQ=255)
                                                                                                                                                                                                 ccGCGCCGGGCGGGCTTGGCGACACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGTACGCGGCCTGGGGTCGATGATTGCGGTAGAGTTTAACGATCCGCAAACGGGCGAGCCGTCAGCGGCGATTGcc  >  1:95277/1‑250 (MQ=255)
                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      
ATTCGCCGCCTGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATGGATCATTACGCCGATAAGCCGGATTTAATGACGATGGCGAAAAGCCTCGCGGGCGGGATGCCGCTTTCGGGCGTGGTCGGTAACGCGAATATTATGGACGCACCCGCGCCGGGCGGGCTTGGCGGCACCTACGCCGGTAACCCGCTGGCGGTGGCTGCCGCGCACGCGGTGCTCAACATTATCGACAAAGAATCACTCTGCGAACGCGCGAATCAACTGGGCCAGCGTCTCAAAAACACGTTGATTGATGCCAAAGAAAGCGTTCCGGCCATTGCTGCGGTACGCGGCCTGGGGTCGATGATTGCGGTAGAGTTTAACGATCCGCAAACGGGCGAGCCGTCAGCGGCGATTGCA  >  NC_000913/1366219‑1366662

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