Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,781,049 C→T intergenic (+32/‑166) lldD → / → trmL L‑lactate dehydrogenase/tRNA (cytidine/uridine‑2'‑O)‑ribose methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,0490CT94.4% 105.6 / ‑3.9 37intergenic (+32/‑166)lldD/trmLL‑lactate dehydrogenase/tRNA (cytidine/uridine‑2'‑O)‑ribose methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (18/16);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACC  >  NC_000913/3780808‑3781289
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
ctgGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTCCTAAATCTGCTAGAAAAAGAGATGAAAGTGGCAATGACGCTGACTGGCGCGAAATCGATTAGCGAAATTTCGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGAATAGGCGTTTGCGCCCCCTCACCCTAAACCTCTCCTTCAggg                                                                                                                                                                                                                                             >  1:255056/1‑248 (MQ=255)
  gggTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGAGTAGCTAACCTGGTAAAACTGATCGAAAAAGAGGTGAAAGAGGGGATGACACTGAATGGCGGGGAATTGAGCAGGGGAATTACGCCAGACTACGTGGGGGAGGGGGTGGGTGAAGAGGTGTCCGCGGCCCAGGGTCCCATGGCGGGAGGGGAGGCGGACTGGTCGTTTGCCGCCCCTCACACCAACCCTCTCCTTCAgggggcgg                                                                                                                                                                                                                                        >  1:267188/1‑246 (MQ=255)
       gTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGAcc                                                                                                                                                                                                                                   <  1:293088/251‑1 (MQ=255)
             tCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCg                                                                                                                                                                                                                              <  1:90219/250‑1 (MQ=255)
                   aTGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCtg                                                                                                                                                                                                                        <  1:279801/250‑1 (MQ=255)
                                  cGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCAct                                                                                                                                                                                                        <  1:84601/251‑1 (MQ=255)
                                       cAGGCGGGTGTAGCTAACCTGCTTAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGGGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGa                                                                                                                                                                                                    >  2:242355/1‑250 (MQ=255)
                                                  aGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCa                                                                                                                                                                                        <  1:130487/251‑1 (MQ=255)
                                                      aaCCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCAct                                                                                                                                                                                      <  1:165884/249‑1 (MQ=255)
                                                        ccTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttt                                                                                                                                                                                  >  1:275065/1‑251 (MQ=255)
                                                              aaaTCTGATCGACAAAGAGATGAAAGAGGCGATGACGCTGACAGGAGCGCAATCAATAAGCGAAATTACGAAAGATTCAATGGTGCAGGGGCTAGGTAAAGAGATGCCTGAGGCACTGGCTCCCATGGCGCAAGGGAGTGCGGCATAGTCGTTTGCGCCCCCTCAACCTAACCCTCTCCTTAAGGAAGAGGGGACAGTTCGGAGCTGTATGTACTCCCTCACTCTGAAACGACACCGCCATCTTTTTTtct                                                                                                                                                                            >  1:39347/1‑251 (MQ=255)
                                                              aaaTCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtct                                                                                                                                                                            >  1:233596/1‑251 (MQ=255)
                                                               aaTCTGAGCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAAGCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCTCCTCACCCTAACCCTCGCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTAGGTACTCCCCCACTCTGAAACGACACCGCCCTCTTTTTTtctc                                                                                                                                                                           <  2:170715/251‑1 (MQ=255)
                                                                   tGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTAGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTc                                                                                                                                                                      >  2:149537/1‑251 (MQ=255)
                                                                       cGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGccc                                                                                                                                                                  <  1:111582/251‑1 (MQ=255)
                                                                       cGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGccc                                                                                                                                                                  >  2:285912/1‑251 (MQ=255)
                                                                       cGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGccc                                                                                                                                                                  <  1:34896/251‑1 (MQ=255)
                                                                        gAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGcccc                                                                                                                                                                 >  1:171899/1‑251 (MQ=255)
                                                                              gagaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                                                                                                                           <  2:62449/251‑1 (MQ=255)
                                                                                   gAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGTGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTCCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                                                                                                                                       <  1:133973/250‑1 (MQ=255)
                                                                                             aTGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTg                                                                                                                                            >  2:65613/1‑251 (MQ=255)
                                                                                               gACGCTGCCTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCtgtatgt                                                                                                                                                                                                                   >  2:301374/1‑179 (MQ=255)
                                                                                               gACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCtgtatgt                                                                                                                                                                                                                   <  1:301374/179‑1 (MQ=255)
                                                                                                  gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGATAGGGCCggg                                                                                                                                               <  2:107468/243‑1 (MQ=255)
                                                                                                  gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGATAGGGCCggg                                                                                                                                               >  1:107468/1‑243 (MQ=255)
                                                                                                                     aTCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCCCTGCTCAc                                                                                                                    >  1:262556/1‑251 (MQ=255)
                                                                                                                                   gCAAGATACGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGCGGGGACCGTTCGGCGCTGTCTGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAACAGGCCGCACTGCTCACCAATTTTCTGCTAt                                                                                                      >  2:106533/1‑251 (MQ=255)
                                                                                                                                             gCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc                                                                                             >  1:253706/1‑251 (MQ=255)
                                                                                                                                                          gggggTAAAGGGTTCCCTGGGGCACGGCCTCCCATGGCAAAAGGGAATGCGGCATAGTCTTTTGCGCCCCCCCCCCCTAACCCTCCCCTTCAGGGGGGGGGGCCCGTCGGGCGCTGTTTGTACTCCCTCACTCTGAAACAACCCCGCACTTTTTTTTTCTCCCCCGCCCCCCCGGGGAGAGGCCCGGGGTGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                                                                <  2:233596/249‑1 (MQ=255)
                                                                                                                                                            gggaaaaGAGTGGCCGGCGGCACTGGCCCCCGTGGGAAAAGGAATGGCGGCAGAGTCGTTGGCCCCCCCCCCCCTTACCCCTCCCCTGCAGGGAGGGGGGCCCGTCGGGCGGTGTATGTACCCCCTCCCTCGGAAACGACCCCGCCTTTTTTTTTCCCCCCTCCCCCCCCCGGGGAGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                                                                              <  2:32083/247‑1 (MQ=255)
                                                                                                                                                                            gCGGCACTGGCCCCCAGGGCGAAAGGAAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCCCCTCCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACCCCCCCACTCTGAAACGACCCCGCCCTTTTTTTTTCCCCCCCCCCCCCCCGGGGAGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTa                                                              <  2:253706/251‑1 (MQ=255)
                                                                                                                                                                                         ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc                                                   >  2:144089/1‑249 (MQ=255)
                                                                                                                                                                                                gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTcc                                          >  2:174196/1‑251 (MQ=255)
                                                                                                                                                                                                aaaaggaagtGCGGCATAGTGTTTGCCCCCCCTCCCCCCAACCCACCCCGTCCGGGGAGGGGGGCCTGTTGGGCGCTGTATGTACCCCCCCGCTCTGAAGCGGCCCCGCCTTCTTTTGTTCCCCCTCGCCCCTCCGGGGAGGGGGCGGGGTTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTTCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTGCTTTACGAACCAGAAATTcc                                          <  2:14502/242‑1 (MQ=255)
                                                                                                                                                                                                                 gagccgttgccgcCCCCTCCCCGTAACTCCCTCCTCCAGGGAGAGGGGACCGTTGGGCGCGTTATGGCCCCCCCCACTCGGAAACGACGCGGCCCTTTTTTTTTCCCCGTCGCCCCCCGGGGGAGAGGGCGGGGGTGAGGGAAAAAGGCCGCACTGCTCTCCAATTTTCTGCGAGCCTGCCCCCGCACTAAGGGGGCAGGATGCTAGACAGCTTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAAc                          <  1:27586/241‑1 (MQ=255)
                                                                                                                                                                                                                          gcgcCCCCTCACCCTAACCCTCTCCTTCAGGGAGGGGGGCCCGTTCGGCGCTGTATGTCCCCCCTCCCTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCg                  <  2:75369/250‑1 (MQ=255)
                                                                                                                                                                                                                                         aaCCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCACTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCAGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATAcc  >  1:35974/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
CTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACC  >  NC_000913/3780808‑3781289

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