Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,617,053 G→C H5D (CAT→GAT)  yjjI ← DUF3029 domain‑containing protein YjjI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,617,0530GC93.9% 151.7 / 2.3 49H5D (CAT→GAT) yjjIDUF3029 domain‑containing protein YjjI
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base C (19/27);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

AGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATGAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGCGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGG  >  NC_000913/4616806‑4617289
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             
agcCATTCGGGATCGTTCGCCAGAAAATGGGCTAAATCGTGTAAGACATTGCGCTGTTTTGTCGGCGCTTGACCTTTAAACAGATCGCAGATGACTTCTTCATCCAGGGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTTTGCTTCCAGGGCCAGAAAATTGCGCGTCGGCTCCGGGCTAAGCACTGGTCTGGAGACAATTTGCTTGCATCGTTGTTGCAGTGCATTGTCATcaga                                                                                                                                                                                                                                            <  2:65698/251‑1 (MQ=255)
       cGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATGTGCTGGCAACGTTGTGGCAGTGCATTTTCATCAGAAGTGGGc                                                                                                                                                                                                                                     <  1:210658/251‑1 (MQ=255)
        ggAACCGTTCGCCAGAAAACGGGCGTAAGCGGGTAAGACATAGCGCGGTTTGTACGGCGCACGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGGGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGTGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCa                                                                                                                                                                                                                                    <  2:367994/251‑1 (MQ=255)
         gAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCAt                                                                                                                                                                                                                                   <  1:262755/251‑1 (MQ=255)
              gTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCAtcttct                                                                                                                                                                                                                              <  1:138898/251‑1 (MQ=255)
                      gAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCttt                                                                                                                                                                                                                      <  1:417893/251‑1 (MQ=255)
                                     cGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACAGATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGTAAGGAAGCTGTGGGTAAGGCCGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTAGGGACAATTTGCTGGCCACGTTGTTGCAGTGCCTTTTCGTCAGAAGTGGGCATCTTCTTTTCCTTTTGTGCCGCAGGTGAt                                                                                                                                                                                                       >  2:378102/1‑251 (MQ=255)
                                         tAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGCAGGTGATgcgc                                                                                                                                                                                                   >  1:632588/1‑251 (MQ=255)
                                                           tGTACAGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTAGGCAGGAAGCTGTTGGAAAGGCAGATTGTTTTCTGCTTCCGGTGCCCGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAGGTGGGCATATTCTTTTCCTTTTATGCCGAAGGGGCTGCTCAATTTTAGGAAGTTTCg                                                                                                                                                                                  >  2:450632/1‑251 (MQ=255)
                                                            gTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCTCGGCGGGCTTCGGCAGGAAGCTGTGGGGAAGGCAGATTGTTTTCTTCTTCCAGTGCCAGAAAATGGCGCTTCTGTTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTTCGTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGt                                                                                                                                                                                <  2:29389/251‑1 (MQ=255)
                                                                       gACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCTATTGTAAGAAGTTTCGTGATGTTCACtt                                                                                                                                                                     <  1:558489/251‑1 (MQ=255)
                                                                            tCAAACCTATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCGTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGa                                                                                                                                                                  <  2:120131/249‑1 (MQ=255)
                                                                             cAAACATATCGCAGGTTACACCTTCAACCAGCTTGCGGCGGGCTTTGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGTCGCTTCTGCTCCGGGCTAAGCATGGGGCTGGTGACAATTTGCTTGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTGATGCTGAAGGTGATGCGCCATTGTAATAAGTGTCGTGATGTTCACTTTGATcc                                                                                                                                                               <  2:640453/251‑1 (MQ=255)
                                                                                    atCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCCGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCACCGTTGTTGCCTTGCATTTTCATCAGAAGGGGGCATCTTCTTTTCCTTTTATGCCGAAGGAGATGTGCCATTGTAAGACGTTTCGGTGTTTTCACTTTGATCTTGATg                                                                                                                                                          >  2:224296/1‑249 (MQ=255)
                                                                                             tACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTTCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTGGCTGGCAACGTTGTTGCAGTGCAGTTTCATCAGAAGTGGGCATCTTCGTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGccac                                                                                                                                                <  1:436829/250‑1 (MQ=255)
                                                                                               cacCTTCATACAGCGCGCGGCGGGCTTCGGCAGGAAGCTGGGGGTATGGCAGATTGTTTTCTGCTTCCATTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGAGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGGGGGCATCTGCTTTTTCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCGGATgaggtggccgccac                                                                                                                                             <  2:632588/251‑14 (MQ=255)
                                                                                                 ccTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTg                                                                                                                                           >  1:586572/1‑251 (MQ=255)
                                                                                                   ttCATACAGAGCGCGGCGGGCTTAGGCAGGACGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCCTCCGAAGTGGGCATCTTCTTTTCATTTTATGACGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACttac                                                                                                                                         >  2:621642/1‑248 (MQ=255)
                                                                                                        ccAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAATATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATGTGCTTGCAACGTTGTTTCAGTTCATTTTCATCAGAAGTGGTCATCTTCTTTTCCTTTTATTCCGAAGGTTATGCGCCATTGTAAGAAGTTTCATGATGTTCACTTTGATCCTGATGAGTTTTACTCCACTGACGcac                                                                                                                                     <  2:247554/250‑2 (MQ=255)
                                                                                                              cgcgGCGGGCTTCGGCAGGAAGCTGTGGGGAAGGCCGATTGTTTTCTGCTTCCCGTGCCAGAAACTGGCGCTTCTGCTCCGGGCTAAACACTGGGCTGGTGACAATTTGCTGGCAACGGTGTTGCAGGGCCTTTTCATCAGAAGTGGGCATCTTCTTTTACTTTTATGCCGAAGGGGATGCGCCATTGTAAGAAGATTCGTGTTGTTCCATTTGATCCTGATGAGTTTGCCACCAATGACGCATTCAtttt                                                                                                                              >  2:585157/1‑250 (MQ=255)
                                                                                                                    gggCGTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTg                                                                                                                        >  2:84906/1‑251 (MQ=255)
                                                                                                                      gctccGGCAGGAAGCTGTGGGTAAGGCCGATTGGTTGCTGCGTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGACATCTTCTTTTCCTTTTATGCCGAAGGGGATGCGCCATTGTAAGAAGGTTCGTTATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGcattcagatcggaagagc                                                                                                                      <  2:7152/247‑13 (MQ=255)
                                                                                                                       cTTCGGCAGGAAGCTGTGTGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATGAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGCGTTTGCCACCACTGACGCATTCATTTGAAAGTGAat                                                                                                                     >  1:264486/1‑251 (MQ=255)
                                                                                                                          cGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAattat                                                                                                                  >  2:602927/1‑251 (MQ=255)
                                                                                                                            gCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTAttt                                                                                                                <  1:21040/251‑1 (MQ=255)
                                                                                                                              aGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCCGTGCCAGAACATGGCGCTTCTGCTCCGGGCTAAGCACTGCGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTACGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGCAAGTGACTTATTTg                                                                                                               >  1:286090/1‑250 (MQ=255)
                                                                                                                              aGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGACAGAAAATGGCGATTCTGCTACGGACTAAGCACTGGGCTGGTGACAATATGCTGGCAACGTTGTTGCAGAGAATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGAAGAAGGTGATGCGCCATTGAACAAAGTTTCGTGATGTTCTATTTGATCCTGCTGAGTGTGCCCCCACTGACGCATCCATTGGAAAGTGAATTATTTg                                                                                                               >  2:299227/1‑250 (MQ=255)
                                                                                                                                  aGCTGTGAGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGCAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAAcc                                                                                                           >  2:578898/1‑250 (MQ=255)
                                                                                                                                       tgGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAAcc                                                                                                           >  1:358858/1‑245 (MQ=255)
                                                                                                                                       tgGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAAcc                                                                                                           >  2:355232/1‑245 (MQ=255)
                                                                                                                                                        tttCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCaa                                                                                     <  1:581960/250‑1 (MQ=255)
                                                                                                                                                           cggCTTCCAGTTCCATAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACTTTGTTTCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTGTATGGCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTTATCCTGATGAGGTTGCCAGCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTATAGTGATGCAAACt                                                                                  <  2:643834/248‑1 (MQ=255)
                                                                                                                                                               ttCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACt                                                                                                                                            <  1:416377/188‑1 (MQ=255)
                                                                                                                                                               ttCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACt                                                                                                                                            >  2:416377/1‑188 (MQ=255)
                                                                                                                                                                   aGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAg                                                                         <  1:578898/251‑1 (MQ=255)
                                                                                                                                                                                      tCTGCTCCGGGCTAAGCACTGGGCGGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGAtg                                                       <  2:20898/250‑1 (MQ=255)
                                                                                                                                                                                         gCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCttt                                                                                                                                                                                                                            >  1:298734/1‑83 (MQ=255)
                                                                                                                                                                                         gCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCttt                                                                                                                                                                                                                            <  2:298734/83‑1 (MQ=255)
                                                                                                                                                                                         gCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGAtgtata                                                   <  1:258327/251‑4 (MQ=255)
                                                                                                                                                                                         gCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGAtgtat                                                    <  1:98734/250‑3 (MQ=255)
                                                                                                                                                                                          cTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATtg                                                            <  1:415650/241‑1 (MQ=255)
                                                                                                                                                                                          cTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATtg                                                            >  2:415650/1‑241 (MQ=255)
                                                                                                                                                                                                gCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATGAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGCGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAgt                                            >  1:214486/1‑251 (MQ=255)
                                                                                                                                                                                                           ggCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCgg                                   >  1:597761/1‑249 (MQ=255)
                                                                                                                                                                                                           ggCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAg                                 <  2:351250/251‑1 (MQ=255)
                                                                                                                                                                                                            gCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGt                                <  2:412322/251‑1 (MQ=255)
                                                                                                                                                                                                                gTGACAATGTGCTGGGAACGTTGTTGCAGTGCTGTTTTATGAGAAGTTGGCATCTTCTCTTCCTTTTATGCCGGAGGTGATGCGCCAGTGTAAGAAGATTCGTGATGTTCACGTTGATCCTGATGCTTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCATATCGCATCACAGTGATGCAAACTTGTCAGTAGAGTTCCTTGATTGTGAAGTGTATCGAAGTGTGTTGCGGATTAGAt                            <  2:284513/251‑1 (MQ=255)
                                                                                                                                                                                                                                aaCGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTTAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACaa            <  1:673577/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gttgCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCGTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCg       <  1:585157/250‑1 (MQ=255)
                                                                                                                                                                                                                                          cAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAgg  >  1:371801/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             
AGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATGAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGCGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGG  >  NC_000913/4616806‑4617289

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