Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,129,487 T→C 5.9% N55S (AAT→AGT)  wcaE ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,129,4870TC 5.9% 209.0 / 2.7 68N55S (AAT→AGT) wcaEputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (29/35);  new base C (2/2);  total (31/37)
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.14e-01

CGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGA  >  NC_000913/2129340‑2129628
                                                                                                                                                   |                                                                                                                                             
cGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCAt                                                                                                                                               <  2:374578/148‑1 (MQ=255)
 ggACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCTTGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCAtt                                                                                                                                              >  8:358938/1‑148 (MQ=255)
  gACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTga                                                                                                                                            <  7:82559/149‑1 (MQ=255)
   aCAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCAtt                                                                                                                                              >  8:114479/1‑146 (MQ=255)
       aaTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATTGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGAttt                                                                                                                                       <  3:309128/149‑1 (MQ=255)
       aaTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCTTGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGAttt                                                                                                                                       <  7:358938/149‑1 (MQ=255)
            gCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCa                                                                                                                                   <  1:201718/148‑1 (MQ=255)
            gCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCa                                                                                                                                   <  1:363783/148‑1 (MQ=255)
                cATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCACTGAGATTTTCCAGATAc                                                                                                                              >  3:220796/1‑149 (MQ=255)
                cATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCACTGAGATTTTCCAGATAc                                                                                                                              >  3:226003/1‑149 (MQ=255)
                   tCTGATGAAAAATATCGCCCAAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCa                                                                                                                           >  3:351840/1‑149 (MQ=255)
                    cTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATGTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCa                                                                                                                           <  4:21180/148‑1 (MQ=255)
                       atgaAAAATATCGCCCGAATTGAGAAACAACGCGAACTTTCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACgg                                                                                                                        <  5:65076/148‑1 (MQ=255)
                                atCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTcc                                                                                                                                    >  4:251621/1‑127 (MQ=255)
                                atCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTcc                                                                                                                                    <  3:251621/127‑1 (MQ=255)
                                 tCGCCCGAATTGAGAAACAACGCGAACTTTCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGtt                                                                                                             <  6:342202/149‑1 (MQ=255)
                                 tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATAGCGTTATCTGGCGCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGtt                                                                                                             <  3:199611/149‑1 (MQ=255)
                                 tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGtt                                                                                                             <  8:337153/149‑1 (MQ=255)
                                   gCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGTCGTCGTTg                                                                                                            <  4:191946/148‑1 (MQ=255)
                                    cccGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTgg                                                                                                           >  7:131901/1‑148 (MQ=255)
                                     ccGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCa                                                                                                                                   >  7:126067/1‑123 (MQ=255)
                                     ccGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCa                                                                                                                                   <  8:126067/123‑1 (MQ=255)
                                     ccGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGaa                                                                                                         >  2:258227/1‑149 (MQ=255)
                                       gAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAAcc                                                                                                       <  1:29840/149‑1 (MQ=255)
                                       gAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGCAAcc                                                                                                       <  7:145301/149‑1 (MQ=255)
                                          ttGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGc                                                                                                     >  4:97405/1‑148 (MQ=255)
                                             agaAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATc                                                                                                 >  5:326364/1‑149 (MQ=255)
                                              gaAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCCTGGCGTCGTAGCTACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTTATGATGCCATTGAGATTTTCCAGATCCTCACGGGTGCCGTCGTTGGCACCGCCATCg                                                                                                >  8:255981/1‑149 (MQ=255)
                                                    acGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGacaac                                                                                           >  3:339040/1‑148 (MQ=255)
                                                    acGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGacaac                                                                                           >  5:186628/1‑148 (MQ=255)
                                                           cTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCAtt                                                                                   >  3:210689/1‑149 (MQ=255)
                                                                cTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAg                                                                              <  1:337262/149‑1 (MQ=255)
                                                                 ttGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGc                                                                             >  7:230171/1‑149 (MQ=255)
                                                                    tgCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGa                                                                          <  7:114479/149‑1 (MQ=255)
                                                                    tgCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCACTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGa                                                                          <  4:220796/149‑1 (MQ=255)
                                                                    tgCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCACTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGa                                                                          <  4:226003/149‑1 (MQ=255)
                                                                      ccATCGCAATGCCTTTGTTCATGGCTTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGata                                                                        <  3:138801/149‑1 (MQ=255)
                                                                       cATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGatat                                                                       >  4:123808/1‑149 (MQ=255)
                                                                       cATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGatat                                                                       <  4:168801/149‑1 (MQ=255)
                                                                           gCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTc                                                                   <  6:186628/149‑1 (MQ=255)
                                                                            cAATGCCTGTGTTCATGGCGTCGTAGATACCGTTGTCTGGCTCGCTGACAAAGCGTAGGTTAGAGATACGGTTGAGATTTTCCTGGTGCTCACGGGTGCCGTCGTTGGAGCGGCCCTCGAGAGCAAGCCGTTCGAAGCTGATGTCttgc                                                                  >  8:381096/1‑147 (MQ=255)
                                                                                 ccTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCt                                                              <  7:381096/148‑1 (MQ=255)
                                                                                  cTTTTTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTgc                                                            <  3:97405/149‑1 (MQ=255)
                                                                                   tttGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTgcg                                                           >  2:284950/1‑149 (MQ=255)
                                                                                         cATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCa                                                                                                                                   >  1:336520/1‑71 (MQ=255)
                                                                                         cATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCa                                                                                                                                   <  2:336520/71‑1 (MQ=255)
                                                                                         cATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGAt                                                     >  2:33708/1‑149 (MQ=255)
                                                                                             gCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGa                                                                                               <  4:248586/103‑1 (MQ=255)
                                                                                             gCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGa                                                                                               >  3:248586/1‑103 (MQ=255)
                                                                                               gtcgtAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCa                                               <  4:339040/149‑1 (MQ=255)
                                                                                               gtcgtAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCa                                               <  4:351840/149‑1 (MQ=255)
                                                                                                       tACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGGGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCa                                       <  6:32901/149‑1 (MQ=255)
                                                                                                                   gCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGa                             <  3:123808/147‑1 (MQ=255)
                                                                                                                     tCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCTCCAGATGCGCCAGCGAGGCATGTGTTTTGACTAt                         >  4:255346/1‑149 (MQ=255)
                                                                                                                     tCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTAt                         >  4:382566/1‑149 (MQ=255)
                                                                                                                        cTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTTCCGTCGTTg                                                                                                            <  1:241803/63‑1 (MQ=255)
                                                                                                                           aCAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTc                   <  1:258227/149‑1 (MQ=255)
                                                                                                                             aaaGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGa                 <  1:33708/149‑1 (MQ=255)
                                                                                                                               aGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAgg               >  7:402763/1‑149 (MQ=255)
                                                                                                                                    aGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCCTTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTAc           >  2:234975/1‑148 (MQ=255)
                                                                                                                                       ttAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAAGCGCCATCGACATCTATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGTGGCTTGTGTTTTGACTTTCACTTCGAGGTTACgtta       >  3:292539/1‑146 (MQ=255)
                                                                                                                                       ttAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGa                 >  2:177918/1‑139 (MQ=255)
                                                                                                                                       ttAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGa                 <  1:177918/139‑1 (MQ=255)
                                                                                                                                       ttAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCAATCTACAACAATCCATTCGAAGCTGATATCTTCCAACTTCGCCAGATTAGCCAGCGAGGCATGTGTTTTGACTATCCCTTAGAGGTTACGaaa       >  3:292536/1‑149 (MQ=255)
                                                                                                                                         aaaGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAAc      <  8:402763/148‑1 (MQ=255)
                                                                                                                                             aTACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGa                                                                          <  1:38776/76‑1 (MQ=255)
                                                                                                                                             aTACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGa                                                                          >  2:38776/1‑76 (MQ=255)
                                                                                                                                             aTACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGa  >  8:1978/1‑148 (MQ=255)
                                                                                                                                             aTACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGa  >  3:301771/1‑148 (MQ=255)
                                                                                                                                                ccATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCa                                                                                     >  5:10795/1‑62 (MQ=255)
                                                                                                                                                ccATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCa                                                                                     <  6:10795/62‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGA  >  NC_000913/2129340‑2129628

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