Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,435,750 C→T G272G (GGG→GGA cpdB ← 2',3'‑cyclic‑nucleotide 2'‑phosphodiesterase/3'‑nucleotidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,435,7500CT98.3% 184.3 / NA 59G272G (GGG→GGAcpdB2',3'‑cyclic‑nucleotide 2'‑phosphodiesterase/3'‑nucleotidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/30);  minor base G (1/0);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCCCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTT  >  NC_000913/4435504‑4436000
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          
cTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCtt                                                                                                                                                                                                                                                        >  2:100769/1‑251 (MQ=255)
     tgacgcccTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGat                                                                                                                                                                                                                                                    >  1:82223/8‑250 (MQ=255)
     tGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTACGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCCTTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGat                                                                                                                                                                                                                                                    >  1:302294/1‑250 (MQ=255)
       aCGAACTGGAGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGTTATCAGCTCCTTCGATATc                                                                                                                                                                                                                                                 >  2:385603/1‑251 (MQ=255)
                  tgtgGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCtt                                                                                                                                                                                                                                      >  2:81074/1‑251 (MQ=255)
                          cGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGAGGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAAGCTTTACCTgg                                                                                                                                                                                                                               <  1:63484/250‑1 (MQ=255)
                                cGGCTTTGATTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTGTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAAc                                                                                                                                                                                                                         <  2:10328/250‑1 (MQ=255)
                                              tCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGCTGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAACACGGCGTGAGCATGGcc                                                                                                                                                                                                          >  1:189918/1‑251 (MQ=255)
                                                  ccAGCTGGCTGTCTGCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCGCACCAAGATGATCGCCCCTCATGCCTGGTATTGCCGCCGGAACATCATTCAGTGTGCCTTGGGCGATATCAGCTCCTTCGATAGCATCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAAc                                                                                                                                                                                                      <  1:32048/251‑1 (MQ=255)
                                                    aGCTTGCTGTCTTCCGCCGCGCGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGACATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGCTATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACAt                                                                                                                                                                                                    >  1:135991/1‑251 (MQ=255)
                                                    aGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACAt                                                                                                                                                                                                    >  2:155085/1‑251 (MQ=255)
                                                    aGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACAt                                                                                                                                                                                                    <  1:422405/251‑1 (MQ=255)
                                                      cTTGCTATCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAACTCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTCCTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGATCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATa                                                                                                                                                                                                   >  2:50861/1‑250 (MQ=255)
                                                        tgctgcattCCGCAGAGAGGGATTTATTATTAGCGATGTCGTAAAGCGGACCTGCTGCCGCTTTCGCCTGCGTCACCAGACATTTACCGCTGTCATAAATGAGTTGTAAGTCGGCCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACAGCATTCAGCGTGCCTTTGGCGATATCAGCGCCTTCGCTATCAGCAAAATCTTTACCTGGGAAAAAGGCGTGAGCATGGCCAAACATAATg                                                                                                                                                                                                >  2:22130/8‑251 (MQ=255)
                                                                    ccgcGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGccc                                                                                                                                                                                    >  1:279721/1‑251 (MQ=255)
                                                                         aGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAAt                                                                                                                                                                               <  1:416658/251‑1 (MQ=255)
                                                                           ggATTTTTTATTAGCGATGTCGTAAATCGGTCGTGTTTCCGCTTTCGCCTGCGTCACCGGCCATTTACCGCGGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCGTTGGCGATATCAGCTCCGTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAAttt                                                                                                                                                                             <  1:290050/251‑1 (MQ=255)
                                                                            gATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTACGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTc                                                                                                                                                                            <  2:302294/251‑1 (MQ=255)
                                                                             aTTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCa                                                                                                                                                                           <  1:145684/251‑1 (MQ=255)
                                                                                 tttATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCTGGATTTCACTg                                                                                                                                                                        <  1:392853/250‑1 (MQ=255)
                                                                                   tattaGCGATGTCGGAAATCGGTCGTGCTTCCGCTTTCGCCTGAGTCACCTGCCATGTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACAGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCCGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATCATGGCGTTAACGCCCGGAATTTCACTGAgg                                                                                                                                                                     >  2:92291/1‑251 (MQ=255)
                                                                                         cGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCTCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAg                                                                                                                                                                                                                                               <  2:392901/171‑1 (MQ=255)
                                                                                         cGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAg                                                                                                                                                                                                                                               >  1:392901/1‑171 (MQ=255)
                                                                                              ccGTAAATCGGTCGTGCGTCCGCTTTCGCCTGCGACCCCTGCCATTTCCCGCTGTCATTACGTAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGTAAAACATAATGGCGTTAACGCCCGGATTTTCATTGAGGTAATAAACTGa                                                                                                                                                          <  2:272870/250‑1 (MQ=255)
                                                                                               cGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAg                                                                                                                                                         <  1:155085/251‑1 (MQ=255)
                                                                                               cGTAAATCGGGCGTGCTTCCGCGTTCGCCTGCGTCACCTGCCATTTCCCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAg                                                                                                                                                         <  2:50122/251‑1 (MQ=255)
                                                                                                gTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGt                                                                                                                                                        <  2:176399/251‑1 (MQ=255)
                                                                                                  aaaTCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCAATTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGtt                                                                                                                                                       <  1:69967/250‑1 (MQ=255)
                                                                                                   aagcGGTCGTGCTTCCGCTTTCGCCTGCGTCTCCTGCCATTTACCGCTGTCATTGCTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGtttt                                                                                                                                                     <  2:31949/248‑1 (MQ=255)
                                                                                                   aaTCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGtttt                                                                                                                                                     >  1:113431/1‑251 (MQ=255)
                                                                                                   aaTCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTGAACGCCCGGAATTTCACTGATGTAATAAACGTAGtttt                                                                                                                                                     <  2:187051/251‑1 (MQ=255)
                                                                                                        gTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGcc                                                                                                                                                >  1:324867/1‑251 (MQ=255)
                                                                                                            ggCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCa                                                                                                                                            <  2:366345/250‑1 (MQ=255)
                                                                                                              cTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCACGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACt                                                                                                                                          >  1:243881/1‑251 (MQ=255)
                                                                                                               ttCCGCTTTCGCCCGCGTCACCTGCCATTTACAGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCAAATGCCTGGCATTACCGCCGGAACACCATTCCGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACtt                                                                                                                                         >  1:160901/1‑251 (MQ=255)
                                                                                                                cccGCTTTCGCCTGCGTCAACTGCCATGTACCGCTGTCATTACTGAGTTGTAAGTCGACCCCGCCAAGATGATCGCCCCACATGCCTTGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAAGTTCACTGAGGTAATAAACTGAGTTTTCCGCGATCACtt                                                                                                                                         <  2:203254/249‑1 (MQ=255)
                                                                                                                 ccGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACttt                                                                                                                                        <  2:113431/250‑1 (MQ=255)
                                                                                                                           cTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCgg                                                                                                                             >  1:11/1‑251 (MQ=255)
                                                                                                                              cGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGc                                                                                                                            >  1:351819/1‑249 (MQ=255)
                                                                                                                                           ttACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAACCCGGAAtg                                                                                                             >  2:357607/1‑251 (MQ=255)
                                                                                                                                                gCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGcc                                                                                                         >  1:452778/1‑250 (MQ=255)
                                                                                                                                                         aCTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGaca                                                                                                <  1:100769/250‑1 (MQ=255)
                                                                                                                                                           tGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGACCGacaac                                                                                              >  1:303660/1‑250 (MQ=255)
                                                                                                                                                                       cGACCACACCAAGACGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCGTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCAcc                                                                                  <  1:466590/250‑1 (MQ=255)
                                                                                                                                                                              acCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCCGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTcgc                                                                          >  1:379548/1‑251 (MQ=255)
                                                                                                                                                                                cAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCgcgc                                                                        >  2:396638/1‑251 (MQ=255)
                                                                                                                                                                                      gatCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCa                                                                  <  2:134052/251‑1 (MQ=255)
                                                                                                                                                                                      gatCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTGTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCa                                                                  <  2:350051/251‑1 (MQ=255)
                                                                                                                                                                                         cGCCCCTCATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAgg                                                                <  1:255955/250‑1 (MQ=255)
                                                                                                                                                                                             ccACATGCCTGGCATTACCTCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACga                                                           <  2:135721/251‑2 (MQ=255)
                                                                                                                                                                                             ccACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGt                                                           <  1:81074/251‑1 (MQ=255)
                                                                                                                                                                                                      tGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACCACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTacgc                                                   >  2:267702/1‑250 (MQ=255)
                                                                                                                                                                                                         cATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACgg                                               <  1:397892/251‑1 (MQ=255)
                                                                                                                                                                                                                        acCATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCAt                                >  1:40717/1‑251 (MQ=255)
                                                                                                                                                                                                                          cATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAATTCTTTACCTGGGAAAACGGCGTGAGCATGGCAAAACATAATGGCGTTAACGACCGGCATTTCACTGAGGTAATACCCTGCGTTTCCCGCCACCACTTTATCCGGCTCGGCAGCTAGCCCGGAATGTGACAGAACGACAACCACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGGACGGTTTCGGTAATATCATTc                              >  1:63298/1‑251 (MQ=255)
                                                                                                                                                                                                                          cATTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCAtgc                              <  2:128021/251‑3 (MQ=255)
                                                                                                                                                                                                                           aTTCAGCGTGCCTTTGGCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAACCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGTTTTCGGTAATATCATTcc                             >  2:409465/1‑250 (MQ=255)
                                                                                                                                                                                                                                         gggCGATATCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCg               <  2:279721/250‑1 (MQ=255)
                                                                                                                                                                                                                                              atatCAGCTCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATa           >  2:424647/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      cccTTCGATATCATCAAAATCATTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACTCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGTATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCTCTTTCCCGGATAAATTAGCtt  <  1:40990/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          
CTTGCTGACGAACTGGCGTGTGGCATCGTGATCGGCTTTGAGTGTTTCTACCAGCTTGCTGTCTTCCGCCGCGAGGGATTTTTTATTAGCGATGTCGTAAATCGGTCGTGCTTCCGCTTTCGCCTGCGTCACCTGCCATTTACCGCTGTCATTACTGAGTTGTAAGTCGACCACACCAAGATGATCGCCCCACATGCCTGGCATTACCGCCGGAACACCATTCAGCGTGCCTTTGGCGATATCAGCCCCTTCGATATCAGCAAAATCTTTACCTGGGAAAACGGCGTGAGCATGGCCAAACATAATGGCGTTAACGCCCGGAATTTCACTGAGGTAATAAACTGAGTTTTCCGCCATCACTTTATACGGATCGGCAGATAGCCCGGAATGTGCCAGAACGACAACAACATCGGCACCTTTCTCGCGCATTTCAGGCACGTATTTGCGCACGGTTTCGGTAATATCATTCACCGTCACTTTCCCGGATAAATTAGCTT  >  NC_000913/4435504‑4436000

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