Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 475,884 C→T T169I (ACC→ATC)  ybaY → PF09619 family lipoprotein YbaY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913475,8840CT97.9% 154.2 / NA 48T169I (ACC→ATC) ybaYPF09619 family lipoprotein YbaY
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (24/22);  minor base G (0/1);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTGCGTACTGAAGGTAAACAGTCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGT  >  NC_000913/475644‑476121
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
cGGTGCGTAATGAAGGTAAACAGTCACAATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAAAGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGATCAGGCCAGCGGTGGCGCAACGACTATCGTACATTCg                                                                                                                                                                                                                                     >  2:225222/1‑251 (MQ=255)
    gCGTACTGAAGGTAAACAGTCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCCCCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACtt                                                                                                                                                                                                                                 <  1:215825/251‑1 (MQ=255)
         cTGAAGGTAAACAGTCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCAccc                                                                                                                                                                                                                            >  2:261440/1‑250 (MQ=255)
                       ccACCATTCAGCTTTGTTCTGTCATTTACCCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGGGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATcc                                                                                                                                                                                                              <  2:361651/250‑1 (MQ=255)
                        ccccATCCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCg                                                                                                                                                                                                             <  2:121262/249‑1 (MQ=255)
                          ccATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATCTTCAGCCGAACGCCCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGTTAAACTGGTATTTATCACAGATACCGTTCAGCCGGTGATCAACCAGGGCGGACCTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTc                                                                                                                                                                                                           >  1:345743/1‑251 (MQ=255)
                                agctttgtgcgttCATTTAACCCGGCATATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGGTTACCGTGAATGACAACCTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGc                                                                                                                                                                                                     <  2:282199/239‑1 (MQ=255)
                                            tCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTcc                                                                                                                                                                                         >  1:227262/1‑251 (MQ=255)
                                                ttaccccGGCAGTTGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACGGGGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCGGACATTGGTGCCGGTTCAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTAc                                                                                                                                                                                      <  1:161077/246‑1 (MQ=255)
                                                ttAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACg                                                                                                                                                                                     >  2:174074/1‑251 (MQ=255)
                                                     ccGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAat                                                                                                                                                                                 >  1:63548/1‑250 (MQ=255)
                                                                       aaCGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAAcccc                                                                                                                                                              >  1:22003/1‑250 (MQ=255)
                                                                          gcgcgTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGGACAGCAATCCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCtctc                                                                                                                                                           <  1:55100/251‑1 (MQ=255)
                                                                          gcgcgTATTCGGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTAGCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACGAAAGCCGACCTGACATTGGTGCCGGGACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCctcgc                                                                                                                                                           <  1:55102/251‑3 (MQ=255)
                                                                                aTTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATAACGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGACAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCg                                                                                                                                                     >  1:43653/1‑251 (MQ=255)
                                                                                            tcgccgATTACCGTGAATAACAAACTGGTATTGATCCCCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGCCCGGACATTGGGGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGGGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGGTAACCCTCTCCGCCCGGAGAGGGTtag                                                                                                                                          <  2:98744/246‑1 (MQ=255)
                                                                                                       cGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCACTATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTCGTAGTTCCAGCg                                                                                                                               >  2:49178/1‑250 (MQ=255)
                                                                                                            aTGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACg                                                                                                                         >  2:283513/1‑251 (MQ=255)
                                                                                                               aCAAACCGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTg                                                                                                                      >  1:95134/1‑251 (MQ=255)
                                                                                                                 aaaCTGGTATTTATCACCGATACCGTTCAGCCGGCGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACACTCTCCGCCCGGAGAGGGGTAGGAGTTCCAGCGCTAACGCTGc                                                                                                                     >  2:250671/1‑250 (MQ=255)
                                                                                                                     gggTATTTATCCCCGATACCTTTCAGCCGGTGATCACCCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTc                                                                                                                <  1:356809/250‑1 (MQ=255)
                                                                                                                      ggTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTc                                                                                                                <  1:210114/250‑1 (MQ=255)
                                                                                                                        tATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGAGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACAGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGAt                                                                                                             >  1:29018/1‑251 (MQ=255)
                                                                                                                         aTTTACCACCGATACCGTTCAGCCGGTGACCAACCAGGGCGGAACTAAAGCCGACCTGACATGGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGAt                                                                                                             <  2:343158/250‑1 (MQ=255)
                                                                                                                           ttATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGAttt                                                                                                           <  1:301367/250‑1 (MQ=255)
                                                                                                                                cccGATACCGTTCAGCCGGTGACCACCCAGGGAGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTAGCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGc                                                                                                     <  2:158345/250‑1 (MQ=255)
                                                                                                                                 ccGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCa                                                                                                                                  >  1:251989/1‑221 (MQ=255)
                                                                                                                                 ccGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCa                                                                                                                                  <  2:251989/221‑1 (MQ=255)
                                                                                                                                   gATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGAACAGCAAACCGCCGTGCCGGTTCAGGCAAGCGGTGGCGAAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGATCCAGCTCCTACGCAATATTAAGCGACGTTAAACCTCTCCGCCCGGAGAGGGATAGTAGTTCCAGCGATAACGCTGCCAGTCGATTTACCCGCTTc                                                                                                  >  2:122837/1‑251 (MQ=255)
                                                                                                                                              gCCGGTGATCACCCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATc                                                                                       <  1:282634/251‑1 (MQ=255)
                                                                                                                                                           cAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCa                                                                          <  1:174074/251‑1 (MQ=255)
                                                                                                                                                             gggCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGa                                                                                                                                                                                                                                    <  2:83785/95‑1 (MQ=255)
                                                                                                                                                             gggCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGa                                                                                                                                                                                                                                    >  1:83785/1‑95 (MQ=255)
                                                                                                                                                                            cGACCTGACAGTGGTGCCGGTACAGCACACCGCCGTGCCGGTTCAGGCCATCGGTGGCGCAACGACAGTAGTACATTCGACTTCACCAACCCAGGTGACTCCGCATTCGGAAGTTCCAGCTCCTACGCAATATTAAGCGACGCTCACCCTCTCCGCCCGGGGAGGGCTAGTAGTTCCACCGATAACGCTGCAAGTCGATGTGTCCGCTTCCCGATACCCTCACTACTTCTGCCAGTAATGCCTGTCGCTg                                                          >  2:92194/1‑250 (MQ=255)
                                                                                                                                                                             gACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGAc                                                        >  2:331433/1‑251 (MQ=255)
                                                                                                                                                                                       tGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtc                                              >  1:349015/1‑251 (MQ=255)
                                                                                                                                                                                            ccGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccggtc                                         >  1:281427/1‑251 (MQ=255)
                                                                                                                                                                                             cggtccaccAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAgccggtcc                                        <  2:227262/243‑1 (MQ=255)
                                                                                                                                                                                             cGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccggtcc                                        <  1:113149/251‑1 (MQ=255)
                                                                                                                                                                                             cGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAgccggtcc                                        <  1:265869/251‑1 (MQ=255)
                                                                                                                                                                                               gTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATTCCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccggtccgg                                      <  1:42132/251‑1 (MQ=255)
                                                                                                                                                                                                           gccgTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATtg                          >  2:266285/1‑251 (MQ=255)
                                                                                                                                                                                                              gTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGc                       >  2:57378/1‑251 (MQ=255)
                                                                                                                                                                                                                   ggTTCAGGCCAGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTgg                  >  1:125440/1‑251 (MQ=255)
                                                                                                                                                                                                                     ttCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCAGCAGTTCCAGCTCCTACGCGATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCg                >  1:162490/1‑251 (MQ=255)
                                                                                                                                                                                                                             aGCGGTGGCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTGCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAacca        <  2:349015/251‑1 (MQ=255)
                                                                                                                                                                                                                                   ggCGTAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCgg   <  2:209449/250‑1 (MQ=255)
                                                                                                                                                                                                                                   ggCGCAACGACTATCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATCACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGt  >  1:210521/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
CGGTGCGTACTGAAGGTAAACAGTCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATCAACCAGGGCGGAACTAAAGCCGACCTGACATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACCCGGT  >  NC_000913/475644‑476121

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