Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,926,543 2 bp→TT coding (236‑237/1365 nt) ppnN → nucleotide 5'‑monophosphate nucleosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,926,5430CT98.1% 174.6 / NA 53P79L (CCC→CTT) ppnNnucleotide 5'‑monophosphate nucleosidase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (29/22);  minor base G (0/1);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,926,5440CT98.1% 170.5 / NA 53P79L (CCC→CTTppnNnucleotide 5'‑monophosphate nucleosidase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (29/23);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

ATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATT  >  NC_000913/2926320‑2926784
                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                                
atTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCAGTACTGAACTACGGTAGTTTGACCGATAACAGCAAAGAATTACTGTCTCGTTTTGAACATTTCGAAATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGCGGCTTTTGTCGATGGGCg                                                                                                                                                                                                                        >  1:585180/1‑251 (MQ=255)
     cccGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAatt                                                                                                                                                                                                                    >  2:456653/1‑250 (MQ=255)
         cTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTcg                                                                                                                                                                                                               >  2:105630/1‑251 (MQ=255)
                ccATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCtt                                                                                                                                                                                                          <  1:203843/249‑1 (MQ=255)
                    ggATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGCAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAgg                                                                                                                                                                                                    >  1:527432/1‑251 (MQ=255)
                          gttgcCGCAGCTGGAATTGGATCTGCTTAACCGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAAc                                                                                                                                                                                               <  2:149792/250‑1 (MQ=255)
                                  aGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTcgcg                                                                                                                                                                                      <  2:454513/251‑1 (MQ=255)
                                      gtaaGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTc                                                                                                                                                                                   <  1:461745/248‑1 (MQ=255)
                                       gAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGGCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGTGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCt                                                                                                                                                                                 <  2:362188/251‑1 (MQ=255)
                                         aGTGGATAGGCTTAACCGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGCTCCTGc                                                                                                                                                                               <  2:101066/251‑1 (MQ=255)
                                               tatGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGAc                                                                                                                                                                          >  2:538660/1‑250 (MQ=255)
                                                     tAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATtctc                                                                                                                                                                    <  1:255737/250‑1 (MQ=255)
                                                       aaCGCACCGCCAGCCGCGAGCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGAAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTc                                                                                                                                                                 >  1:444216/1‑251 (MQ=255)
                                                        aCGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGGAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTGGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCg                                                                                                                                                                <  2:245638/251‑1 (MQ=255)
                                                             ccGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTAc                                                                                                                                                           <  2:355383/251‑1 (MQ=255)
                                                               gCCAGCATCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAATTTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGAATAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACgg                                                                                                                                                         <  1:209812/251‑1 (MQ=255)
                                                                ccAGCCGCGACCTCTATCATCTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTGGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACggg                                                                                                                                                        <  2:417816/251‑1 (MQ=255)
                                                                  agcagcGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCa                                                                                                                                                      <  1:456653/251‑1 (MQ=255)
                                                                   gcagcGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCa                                                                                                                                                      >  1:70513/1‑250 (MQ=255)
                                                                       cGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATcc                                                                                                                                                 >  2:83654/1‑251 (MQ=255)
                                                                         aCCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCAt                                                                                                                                               >  1:569782/1‑251 (MQ=255)
                                                                                tCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCCTCGTTTACGGGCAAATCCATAACACcg                                                                                                                                        >  1:540523/1‑251 (MQ=255)
                                                                                           gCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTccc                                                                                                                             <  1:83654/251‑1 (MQ=255)
                                                                                           gCAACTGTTCACTTGCCGAACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGAGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGAGTGACATTCTCTTCGTTTACGGGCAAATCCCTAACACCGTTCGttttacc                                                                                                                             >  2:234430/1‑248 (MQ=255)
                                                                                              aCTGTTCCCTTGCCGTACTTAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTGGTCTATGGGCGAATTATGCGCTCTGTGCAGGCCAACTTGTTCGCGGTCCTGCGTGTCATGCTCTTCGTTTACGGGCAAATCCATAACACCGGTCGTTTGCCCatc                                                                                                                          <  2:353053/251‑3 (MQ=255)
                                                                                                    cACTTGTCGTACTGAACGCCTGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTGTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAAATTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTAGTTTTCCCAACCTGAAt                                                                                                                    <  2:416943/251‑1 (MQ=255)
                                                                                                        tGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCGGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAAtct                                                                                                                  <  2:303108/249‑1 (MQ=255)
                                                                                                               cTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTc                                                                                                         >  1:518906/1‑251 (MQ=255)
                                                                                                                 gAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGTCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTc                                                                                                         <  1:346189/249‑1 (MQ=255)
                                                                                                                 gAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGATGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTTCGGGCAAATCCATAAAACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCg                                                                                                       >  2:124349/1‑251 (MQ=255)
                                                                                                                  aacccgggtggTTTGACCGATAACAGCAGAGCATTGCTGTATGGTTTTGAAAATTTAGAAATTAACGTCTTGCGCCGTAAACGCGGCGTAAGGCTGGACTTGATTACTCTTCCGGAAGAGTCTTTTGACGATGGGCGAATTATTCTCGATTTGCAGGCCGACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACCGGCACATCCATATCAATGTTCGTTGTCCCCACCTGCATCTCGACAACTCCGt                                                                                                      >  2:51884/11‑251 (MQ=255)
                                                                                                                  aaCTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGTCAAATCCATAACACCGTTCGTTTTccccagatcagaacccgtctgaactc                                                                                                      >  1:426573/1‑228 (MQ=255)
                                                                                                                        ggTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACAt                                                                                                <  1:538660/251‑1 (MQ=255)
                                                                                                                         gTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATc                                                                                               >  1:375445/1‑251 (MQ=255)
                                                                                                                                gACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAAcc                                                                                        >  1:560906/1‑251 (MQ=255)
                                                                                                                                       aaCAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCtt                                                                                 <  2:222849/251‑1 (MQ=255)
                                                                                                                                              aaGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATc                                                                          >  1:269659/1‑251 (MQ=255)
                                                                                                                                                gAATTGCTGTCTCGTTTTGAAACTTTCGCTATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACGCCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTCACCTGGTCTTTTCCATCtt                                                                        >  2:166246/1‑251 (MQ=255)
                                                                                                                                                     gctgccccgtgttGAAAAGTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGAGGGTCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTTGTCTTTTCCATCTTGCGt                                                                    <  2:396514/239‑1 (MQ=255)
                                                                                                                                                       tGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCAACATCACTAACCTGGTCTTTTCCATCTTGCGTAAc                                                                 >  2:415685/1‑251 (MQ=255)
                                                                                                                                                             gTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTcgc                                                           >  2:423675/1‑251 (MQ=255)
                                                                                                                                                                  acgctcttccgatctAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCt                                                      >  1:38623/15‑251 (MQ=255)
                                                                                                                                                                     caattCGATATTAACGTCTTGAGCCGTGAACGCGGCGTAACGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTACGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGc                                                    >  2:183880/4‑250 (MQ=255)
                                                                                                                                                                               ttAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGa                                          >  1:362230/1‑250 (MQ=255)
                                                                                                                                                                                           gcCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCACATAtg                             >  2:447072/1‑251 (MQ=255)
                                                                                                                                                                                                  aCGCGGCGTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTct                      >  2:393805/1‑251 (MQ=255)
                                                                                                                                                                                                   cgcgGCGTAAAGCTGGAACTGATTAATCGTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCGCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCGTGCGTAACTCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTctg                     <  1:189764/251‑1 (MQ=255)
                                                                                                                                                                                                      ggcgTAAAGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTctgctg                  >  1:315448/1‑251 (MQ=255)
                                                                                                                                                                                                         gTAAAGCTGGAACTGATTAATCTTCCGGAATAGGCTTTTGTCGATGAGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACATGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTgggg               >  1:174498/1‑251 (MQ=255)
                                                                                                                                                                                                             aGCTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGt           <  2:121302/251‑1 (MQ=255)
                                                                                                                                                                                                               cTGGAACTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTc          >  1:429924/1‑250 (MQ=255)
                                                                                                                                                                                                                   aaCTGATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTc      >  2:20663/1‑250 (MQ=255)
                                                                                                                                                                                                                       gATTAATCTTCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAAtt  <  1:20663/250‑1 (MQ=255)
                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                                
ATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGATTAATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATT  >  NC_000913/2926320‑2926784

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: