Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,741,228 C→T D119N (GAT→AAT)  ysaA ← putative electron transport protein YsaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,741,2280CT98.4% 206.7 / NA 61D119N (GAT→AAT) ysaAputative electron transport protein YsaA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (37/23);  minor base G (0/1);  total (37/24)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GCCTACATGATTTCTGCAATATATTGAATTTGGAAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCTGGCATAAACAAAGCGCACTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATCGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCT  >  NC_000913/3740987‑3741472
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
gCCTACATGATTTCTGCAATATATTGAATTTGGAAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCTGGCATAAACAAAGCGCACTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAAt                                                                                                                                                                                                                                             >  1:359413/1‑251 (MQ=255)
        gATTTCTGCAATATATTGAATTTGGAAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCTGGCATAAACAAAGCGCACTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCg                                                                                                                                                                                                                                     >  1:117898/1‑251 (MQ=255)
        gATTTCTGCAATATATTGAATTTGGAAGAAGTTGTAGGCCGGATAAGGCGTTTACGCCGCATCTGGCATAACCAAAGCGCACTGTGCCGGCGGTGAGAAACGCCGGAAAACGTTGCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACGGCAACGCCTTTGGCGGGCAGGCTTCAACACAGGCCTGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTc                                                                                                                                                                                                                                      <  2:407550/250‑1 (MQ=255)
           ttCTGCAATATATTAAATTTGGAAGACTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCTGGCATAAACAAAGCGCACTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCCACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAAGCGCCCTCGCCt                                                                                                                                                                                                                                  >  2:79365/1‑251 (MQ=255)
               gCAATATATTGAATTTGGAAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCTGGCATAAACAAAGCGCACTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCttt                                                                                                                                                                                                                                <  1:200103/249‑1 (MQ=255)
               gCAATATATTGAATTTGGAAGAATTTGTAGGCCGGATAAGGCGATTACGCCGCATCTGGCATAAACAAAGCGCACTGTTCCGGCGTTGCGAAACGCCGGAAAACGTTTAAAACAGGCTGCTGCCGTAGCCGATGCAGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCAGTCTCCCGATGCCAGCACAGCTTGCACTTAATCGCCCTCGCCttttt                                                                                                                                                                                                                              >  2:211247/1‑251 (MQ=255)
                 aaTATATTGAATTTGGAAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCTGGCATAAACAAAGCGCACTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTgc                                                                                                                                                                                                                            <  1:109971/251‑1 (MQ=255)
                    atatTGAATTTGGAAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCTGGCATAAACAAAGCGCCCTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGcgcg                                                                                                                                                                                                                         >  1:189214/1‑251 (MQ=255)
                     tatTGAATTTGGAAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCTGGCATAAACAAAGCGCACTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGcgcg                                                                                                                                                                                                                         >  2:218507/1‑250 (MQ=255)
                       ttGAATTTGGATGAATTTGAACGCCGGATACGGCGTTTACGCCGCATCTGACATAAACACAGCGCACTGTTCCGGCGTTGAGCAACGCCGGCAAACGTTTCAAACAGGCTGCTGCCGTCGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACCCCGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCga                                                                                                                                                                                                                        >  2:57298/1‑249 (MQ=255)
                       ttGAATTTGGAAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCTGGCATAAACAAAGCGCACTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCga                                                                                                                                                                                                                        >  2:302389/1‑249 (MQ=255)
                              tGGAAGAATTTGTAGGCCGGATACGGCGTTTACGCCGCATCTGGCATAAACAAAGCGCACTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTAAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGACGAAACGAc                                                                                                                                                                                                               >  2:54837/1‑251 (MQ=255)
                               ggAAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCTGGCATAAACAAAGCGCACTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGAc                                                                                                                                                                                                               <  1:250141/250‑1 (MQ=255)
                                   gaaTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCTGGCATAACCAAAGCGCACTGTTCCGGCGTTGAGAAACGCCGGACAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTcc                                                                                                                                                                                                           >  2:104697/1‑250 (MQ=255)
                                     aTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCTGGCATAAACAAAGCGCACTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATc                                                                                                                                                                                                        <  1:331600/251‑1 (MQ=255)
                                        tGTAGGCCGGATACGGCGTTTACGCCGCATCTGGCATAAACAAAGCGCACTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTCAACCAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGc                                                                                                                                                                                                      >  2:290069/1‑250 (MQ=255)
                                           aGGCCGGATAAGGCGTTTACGCCGCATCTGGCATAAACAAAGCGCACTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCAcc                                                                                                                                                                                                   >  1:132916/1‑250 (MQ=255)
                                                                   cATCTGGCATAAACAAAGCGCACTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGTCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCtt                                                                                                                                                                          <  1:54837/251‑1 (MQ=255)
                                                                   cATCTGGCATAAAAAAAGCGCACTGTTCCGGCGTTGAGAAACGACGGAACACGTTTAAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGCTGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACAAAACGGTCAAACCAGCATACCGCtt                                                                                                                                                                          >  2:350186/1‑251 (MQ=255)
                                                                        ggCATAAACAAAGCGCACTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTAc                                                                                                                                                                      >  2:290905/1‑250 (MQ=255)
                                                                         gCATAAACCAAGAGCCCTGTTCCGGCGTTGAGCAACGACGCAAAACGTTTCACACAGGCTGCTGCCGTAGCCGCTGAAGCTGCACTTTCTCGCCATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGACAGCACAGATTGCAATTAATCGCCCTCTCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCCAGCCAGCATACAGCTTTTAc                                                                                                                                                                      >  2:205563/1‑249 (MQ=255)
                                                                            aaaaCAAAGCGCACTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCGGCCGTAGCCTATGCCGCTGCACTTTCCCGACATCCCAGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAAGCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCa                                                                                                                                                                 <  2:119921/250‑1 (MQ=255)
                                                                             aaaCAAAGCGCACTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCa                                                                                                                                                                 >  1:296769/1‑250 (MQ=255)
                                                                                      gcacggttccgaggtggtgaAACGCCGTAAAACGTTTCAAACAGGCTGCTGCCGTAGCAGAAGCCGCTGCAATTTCGCGACATCCATGCACTGCACCGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCGCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTGACAGCCAATGCAACgtg                                                                                                                                                       <  2:175544/233‑1 (MQ=255)
                                                                                       cACTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGAACTGCAACGCCTTTGTCTGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACgtgt                                                                                                                                                      >  1:269896/1‑251 (MQ=255)
                                                                                                ggCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGCCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAACAGGGCAAGCCAGAATAACGCTTTTACAGCCCCTGCACCGtgttttatcaac                                                                                                                                             >  2:253219/1‑244 (MQ=255)
                                                                                                 gCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGGCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGCAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACg                                                                                                                                            >  1:278277/1‑251 (MQ=255)
                                                                                                  cGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACg                                                                                                                                            >  1:51418/1‑250 (MQ=255)
                                                                                                       agaAACGCCGGAAACCGTTTCAACCAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCTAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGACAat                                                                                                                                       >  1:108542/1‑250 (MQ=255)
                                                                                                        gaAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGATCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAatat                                                                                                                                     >  1:55922/1‑251 (MQ=255)
                                                                                                               cGGAAAACTTTTCAAATAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCGTTGTCGGGCAGGCTTCAAAACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATg                                                                                                                              <  2:45579/251‑1 (MQ=255)
                                                                                                                       gTTTCAAACAGACTGCTGCCGTAGCCGATGCCGCTGCACTTTCTAGAAATCCATGCACTGCAACGCCATTGTCGGCCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATAGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCCACAAACGGGCAAGACAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAACATGCCACTGTTCGCGGc                                                                                                                      >  1:213148/1‑251 (MQ=255)
                                                                                                                          tCAAACAGGCTGCTGCCGTCGCCGCTGCCGCTGCACTTTCTCGACCTACATGCCCTGGACCGCATTTGCCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGCTGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACAGCCTCCATCGCACCACACGGGCAAGCCAGCATACAGCTTTTCCAGCCAATGCAACGAGTTTGTTCAACGAAAATATGCCCATGTTcgggatcgg                                                                                                                   >  2:21205/1‑244 (MQ=255)
                                                                                                                          tCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACAGCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGAGGCCAGCACAGATTGCACTTAATCGCCCGCGCCTTTGTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTcg                                                                                                                          <  1:21205/244‑1 (MQ=255)
                                                                                                                              aCAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATcgc                                                                                                               <  2:153635/251‑1 (MQ=255)
                                                                                                                                     gctgcCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCCTACAGCTTTTACAGCCAATGCAACGAGTTTGTTCAACGAAAATCTGCCCATGTTCGCGGCTTCTCGCGTCAAc                                                                                                         >  2:281020/1‑250 (MQ=255)
                                                                                                                                             aGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGa                                                                                                 <  1:38475/250‑1 (MQ=255)
                                                                                                                                                 gATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCATCAGGGCAGACATTc                                                                                            >  1:283452/1‑251 (MQ=255)
                                                                                                                                                    gCCGCTGCACTTTCCCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGTCAGACATTcgcg                                                                                         <  2:359413/251‑1 (MQ=255)
                                                                                                                                                            aCTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGc                                                                                 >  1:109685/1‑251 (MQ=255)
                                                                                                                                                                       aTCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACt                                                                      >  1:351965/1‑251 (MQ=255)
                                                                                                                                                                           aTGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACtgatg                                                                  <  2:160054/251‑1 (MQ=255)
                                                                                                                                                                               aCTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGAc                                                                >  1:412056/1‑249 (MQ=255)
                                                                                                                                                                                cTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCAGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGAAATTCGCGCACGGTGCCTCTTCACACTGATGACAgg                                                             >  2:256268/1‑251 (MQ=255)
                                                                                                                                                                                  gccaCGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAgg                                                             <  1:28485/246‑1 (MQ=255)
                                                                                                                                                                                       gCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCCCAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACCCTGATGACAGGCTACTGc                                                      >  1:408068/1‑251 (MQ=255)
                                                                                                                                                                                                   cAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCcagc                                           >  2:117966/1‑250 (MQ=255)
                                                                                                                                                                                                    aGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCcagc                                           >  1:15690/1‑249 (MQ=255)
                                                                                                                                                                                                     ggCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCCCAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGt                                        <  2:109685/251‑1 (MQ=255)
                                                                                                                                                                                                       cTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTgg                                      >  1:255695/1‑251 (MQ=255)
                                                                                                                                                                                                          cAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCCGTGGTCt                                   >  2:11863/1‑251 (MQ=255)
                                                                                                                                                                                                            acacAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGACGGGGTCCAGCAGTGGTCttt                                 >  2:25889/1‑251 (MQ=255)
                                                                                                                                                                                                              acaGGCCGGTCCCGTCTCCCGATGCCAGCACAGATGGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAACCGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGGTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCTCACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTaa                               <  1:276789/251‑1 (MQ=255)
                                                                                                                                                                                                                        cccGTCTACCGATGCCAGCACAGATTGCACTTAAACGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCAACGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGaa                      <  2:94717/250‑1 (MQ=255)
                                                                                                                                                                                                                           gTCTCCCGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGGGGTCTTTAATGACACGAATACg                  <  1:324774/251‑1 (MQ=255)
                                                                                                                                                                                                                                ccGATGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCAGCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGaaa             <  2:15690/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aTGCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCGGTGGTCTTTAATGACACGAATACGGGAAATaa          <  2:117898/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gCCAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAAc        <  1:7347/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAACCGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTc      <  2:353346/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cAGCACAGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTc      <  1:273545/251‑1 (MQ=255)
                                                                                                                                                                                                                                           acaGATTGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCt  >  1:85857/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
GCCTACATGATTTCTGCAATATATTGAATTTGGAAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCTGGCATAAACAAAGCGCACTGTTCCGGCGTTGAGAAACGCCGGAAAACGTTTCAAACAGGCTGCTGCCGTAGCCGATGCCGCTGCACTTTCTCGACATCCATGCACTGCAACGCCTTTGTCGGGCAGGCTTCAACACAGGCCGGTCCCGTCTCCCGATGCCAGCACAGATCGCACTTAATCGCCCTCGCCTTTTTGCGCGAAGAAACGACCTCCATCGCACCAAACGGGCAAGCCAGCATACAGCTTTTACAGCCAATGCAACGTGTTTGTTCAACGAAAATATGCCCATGTTCGCGGCTTATCGCGTCAACAGGGCAGACATTCGCGCACGGTGCATCTTCACACTGATGACAGGCTACTGCCGTGGTCCAGCAGTGGTCTTTAATGACACGAATACGGGAAATAAACTCGTCT  >  NC_000913/3740987‑3741472

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