Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,164,123 C→G A162G (GCA→GGA)  btuB → cobalamin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,164,1230CG97.6% 136.5 / ‑3.0 42A162G (GCA→GGA) btuBcobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (21/20);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

AGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTGCGCGGCTATGGCT  >  NC_000913/4163883‑4164365
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
aGCTTTCATCTAGTTTTAGTCGCGGTACAAATGCCAGTCAGGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTGGTCCAGCGTGTTGAATAGATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGGGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGa                                                                                                                                                                                                                                          <  2:59728/251‑1 (MQ=255)
 ttgctcATCTATTTTTATGCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGa                                                                                                                                                                                                                                          <  2:75644/247‑1 (MQ=255)
   tctcATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAgggggggggagg                                                                                                                                                                                                                                       >  1:75644/1‑247 (MQ=255)
    ctcATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGGGGATATAATCACGACGCGCGATGAACCCGGAACGGAAATTTAAGGAGGGTGGGGAAGc                                                                                                                                                                                                                                       >  1:216532/1‑250 (MQ=255)
           aTTTTTATTCGCGGTACAACTGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGGGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGt                                                                                                                                                                                                                                 <  1:322856/249‑1 (MQ=255)
                   tCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGAGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGTGAATATCATCACGACGGGGGATGAACCAGGAACGGAGATTTCAGGGGGGTGGGGACGAAATAGTTATCAGAACt                                                                                                                                                                                                                       >  1:166089/1‑251 (MQ=255)
                               tGCCAGTCATGTGTTGGTGTTAATTGAGGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTACGATGCAATAGGCGGGGAGGTGAATATCCTCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGGGGGTGGGGAAGCAATAGATATCCAAACTATGATGTctct                                                                                                                                                                                                            >  1:303512/1‑250 (MQ=255)
                                         gtgtTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaac                                                                                                                                                                                                  >  1:20174/1‑250 (MQ=255)
                                               gtgtGAATTGATGGAGTACGCCTGAATCTGGCGGGGGTGAGGGGTTCTGCCGACCTTCGCCAGGTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGGGGGCCGAGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaacatcaggg                                                                                                                                                                                           <  2:331816/251‑7 (MQ=255)
                                                  ttAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTggggg                                                                                                                                                                                         >  1:172312/1‑250 (MQ=255)
                                                      ttGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGACTATATCCGTGGGCCGAGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGGGAATATCATCACGACACGCGATGAACCCGGAACGGAAAATTCCGGAGGGTGGGGAAGCAATCGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATa                                                                                                                                                                                      >  1:417413/1‑249 (MQ=255)
                                                      ttGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATa                                                                                                                                                                                      >  1:114833/1‑249 (MQ=255)
                                                                 cGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGGGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAAc                                                                                                                                                                         >  1:119618/1‑251 (MQ=255)
                                                                           tGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGGGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGTGAATATCATCACGAAGCGCGATGAAACCGGAACGGAAATTTCAGGGGGGGGGGGAAGCCATAGTTATCAGAACTTTGATGTCTCTACGCAGCAACCACTAGGGGGTAAGACACGGGTAACGCTGTTgggg                                                                                                                                                               >  1:2654/1‑250 (MQ=255)
                                                                                ggggAGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGAAGCGCGATGAACCCGGAACGGGAATTTCAGGAGGGTGGGGGAGCAATAGTTATAAGAACTATGATTTCTCTACGCAGCAACAACTGGGGGATAAGAAACGGGTAACGCTGTTgggagaat                                                                                                                                                           >  2:219476/1‑245 (MQ=255)
                                                                                  ggTGAGTTGTTCTGCCGACATTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGTCGCGCGACGCAACCGGAACGGGAAATTCAGGCGGGTGGGGAAGGATTAGTTATCTGAACTATGATTTCTCTACGAGGCACCAACTGGGGGATAAGAAACGGGTaagcattatggtgggttaat                                                                                                                                                        >  2:83419/1‑233 (MQ=255)
                                                                                               gCCGTCCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATgg                                                                                                                                           <  2:126228/251‑1 (MQ=255)
                                                                                               gCCGACCTTAGCCAGGTCCCTATTTCGCTTGTCCAGCGTGTTTAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATGTCAGGAGGGTGGGGAAGCAATAGGTATCAGAGCTATGATGTCTCTACGCATCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCAt                                                                                                                                             <  1:196910/249‑1 (MQ=255)
                                                                                                           cAGTTACCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGCGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAAATTCAGGAGGGTGGGGAAGCAATAGATATCAGAACTATGATGTCTCTACGCCGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATTCCCATCCTCATGGTTATGAtgttgt                                                                                                                               >  1:270197/1‑251 (MQ=255)
                                                                                                             gTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGttgttg                                                                                                                             <  1:38372/251‑1 (MQ=255)
                                                                                                                     ttGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGGGGGTGGGGGAGCAATAGGTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGGTTATGCCCATACTCATGGTTATGAGTTTGTTGCCTTTggg                                                                                                                     >  2:270391/1‑250 (MQ=255)
                                                                                                                          cTTGTCCAGCGTGTTGAATATATCAGTGGGCCGCGCACCGCTGGTTATGGTTCCGATGCAATAGGCGGAGGGGTGAATATCATCGCGCCGCGCGATGAACCCGGAACGGGAAGTTCAGGTGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGAAACCACTGGGGGAGAAGACCCGGGTAACGCAGTTGGGAGGATATGCCCCAACTCCTGGTTATGGTGTTGTTGCATATGGTAATAc                                                                                                                >  2:244485/1‑251 (MQ=255)
                                                                                                                            tGTCCAGCGGTTTGAATAGTTCCGTGGGCCGCGCTCCTCTGTTTATGTTTCCGATGCAATAGGCGGGGTGGTTAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCATCATCTGGGGGATAATACACGGGTAACGCTTTTGGTCGATTATGCCCATACGCATGGTTATGATGTTGTTGCCTATCGTAATAcc                                                                                                               <  2:390912/250‑1 (MQ=255)
                                                                                                                               ccAGAGTGTTGAATATATTCGTGGGCCGCGCTCCGCGGTTTAGGGTGCCGATGCAATAGGCGGGGTGGTGAATATCATCACGCCGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGTCACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGaa                                                                                                           <  2:336612/251‑1 (MQ=255)
                                                                                                                                 aGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGcc                                                                                                                           <  1:253845/233‑1 (MQ=255)
                                                                                                                                 aGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGGGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCCACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCCTACTCATGGTTTTGATGttgttggcagatgggaagagggtcgg                                                                                                         >  2:253845/1‑231 (MQ=255)
                                                                                                                                       ttGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAgcg                                                                                                   >  1:60548/1‑251 (MQ=255)
                                                                                                                                        tGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCg                                                                                                              >  1:315418/1‑239 (MQ=255)
                                                                                                                                        tGAATATAGCCGTGGGCCGCGCTCCTCTGTTTATGGTTCCGATGCAAGAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCg                                                                                                              <  2:315418/239‑1 (MQ=255)
                                                                                                                                           atatatCCGTGGGCCGCGCTCCGCTGTTTATGGTGCCGATGCAATAGGCGGGGTGGGGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGcaga                                                                                               <  2:211025/251‑1 (MQ=255)
                                                                                                                                            tatatCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGACTATCATCACGACGCGAGATGAACCCGGAACGGAACTTTCAGGTGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGcagac                                                                                              >  1:49061/1‑251 (MQ=255)
                                                                                                                                               atCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTGGCCTATGGTAATACCGGAACGCAAGCGcagacag                                                                                            <  1:394827/250‑1 (MQ=255)
                                                                                                                                                 ccGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGAt                                                                                          <  1:289789/250‑1 (MQ=255)
                                                                                                                                                      ggCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAATCGCAGACAGATAACGAt                                                                                    <  2:114833/251‑1 (MQ=255)
                                                                                                                                                                                gATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATgg                                                          <  1:150243/251‑1 (MQ=255)
                                                                                                                                                                                                   tgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACttt                                       >  1:53951/1‑251 (MQ=255)
                                                                                                                                                                                                     aaTATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTAGGGCGCGCTGGAGCATAACTTTAc                                     <  2:60548/251‑1 (MQ=255)
                                                                                                                                                                                                        atcatcACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTg                                   >  1:318962/1‑250 (MQ=255)
                                                                                                                                                                                                          catcaCGACGCGTGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATg                                <  2:303512/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aaTTTCAGGCGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATaa                                           >  2:295870/1‑210 (MQ=255)
                                                                                                                                                                                                                                        aaTTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATaa                                           <  1:295870/210‑1 (MQ=255)
                                                                                                                                                                                                                                        aaTTGCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTGCGCggctatggct  <  1:270391/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
AGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTGCGCGGCTATGGCT  >  NC_000913/4163883‑4164365

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