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

AGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGC‑GGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATGAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGCGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGG  >  NC_000913/4616806‑4617289
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             
agcCATTCGGGATCGTTCGCCAGAAAATGGGCTAAATCGTGTAAGACATTGCGCTGTTTTGT‑CGGCGCTTGACCTTTAAACAGATCGCAGATGACTTCTTCATCCAGGGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTTTGCTTCCAGGGCCAGAAAATTGCGCGTCGGCTCCGGGCTAAGCACTGGTCTGGAGACAATTTGCTTGCATCGTTGTTGCAGTGCATTGTCATcaga                                                                                                                                                                                                                                            <  2:65698/251‑1 (MQ=255)
       cGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGC‑GGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATGTGCTGGCAACGTTGTGGCAGTGCATTTTCATCAGAAGTGGGc                                                                                                                                                                                                                                     <  1:210658/251‑1 (MQ=255)
        ggAACCGTTCGCCAGAAAACGGGCGTAAGCGGGTAAGACATAGCGC‑GGTTTGTACGGCGCACGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGGGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGTGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCa                                                                                                                                                                                                                                    <  2:367994/251‑1 (MQ=255)
         gAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGC‑GGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCAt                                                                                                                                                                                                                                   <  1:262755/251‑1 (MQ=255)
              gTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGC‑GGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCAtcttct                                                                                                                                                                                                                              <  1:138898/251‑1 (MQ=255)
                      gAAAACGGGCGTAATCGGGTAAGACATAGCGC‑GGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCttt                                                                                                                                                                                                                      <  1:417893/251‑1 (MQ=255)
                                     cGGGTAAGACATAGCGC‑GGTTTGTACGGCGCATGACCTTCAAACAGATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGTAAGGAAGCTGTGGGTAAGGCCGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTAGGGACAATTTGCTGGCCACGTTGTTGCAGTGCCTTTTCGTCAGAAGTGGGCATCTTCTTTTCCTTTTGTGCCGCAGGTGAt                                                                                                                                                                                                       >  2:378102/1‑251 (MQ=255)
                                         tAAGACATAGCGC‑GGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGCAGGTGATgcgc                                                                                                                                                                                                   >  1:632588/1‑251 (MQ=255)
                                                           tGTACAGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTAGGCAGGAAGCTGTTGGAAAGGCAGATTGTTTTCTGCTTCCGGTGCCCGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAGGTGGGCATATTCTTTTCCTTTTATGCCGAAGGGGCTGCTCAATTTTAGGAAGTTTCg                                                                                                                                                                                  >  2:450632/1‑251 (MQ=255)
                                                            gTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCTCGGCGGGCTTCGGCAGGAAGCTGTGGGGAAGGCAGATTGTTTTCTTCTTCCAGTGCCAGAAAATGGCGCTTCTGTTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTTCGTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑t                                                                                                                                                                                <  2:29389/251‑1 (MQ=255)
                                                                       gACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCTATTGTAAGAAGTTTCG‑TGATGTTCACtt                                                                                                                                                                     <  1:558489/251‑1 (MQ=255)
                                                                            tCAAACCTATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCGTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGa                                                                                                                                                                  <  2:120131/249‑1 (MQ=255)
                                                                             cAAACATATCGCAGGTTACACCTTCAACCAGCTTGCGGCGGGCTTTGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGTCGCTTCTGCTCCGGGCTAAGCATGGGGCTGGTGACAATTTGCTTGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTGATGCTGAAGGTGATGCGCCATTGTAATAAGTGTCG‑TGATGTTCACTTTGATcc                                                                                                                                                               <  2:640453/251‑1 (MQ=255)
                                                                                    atCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCCGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCACCGTTGTTGCCTTGCATTTTCATCAGAAGGGGGCATCTTCTTTTCCTTTTATGCCGAAGGAGATGTGCCATTGTAAGACGTTTCGGTG‑TTTTCACTTTGATCTTGATg                                                                                                                                                          >  2:224296/1‑249 (MQ=255)
                                                                                             tACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTTCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTGGCTGGCAACGTTGTTGCAGTGCAGTTTCATCAGAAGTGGGCATCTTCGTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGccac                                                                                                                                                <  1:436829/250‑1 (MQ=255)
                                                                                               cacCTTCATACAGCGCGCGGCGGGCTTCGGCAGGAAGCTGGGGGTATGGCAGATTGTTTTCTGCTTCCATTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGAGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGGGGGCATCTGCTTTTTCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCGGATgaggtggccgccac                                                                                                                                             <  2:632588/251‑14 (MQ=255)
                                                                                                 ccTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTg                                                                                                                                           >  1:586572/1‑251 (MQ=255)
                                                                                                   ttCATACAGAGCGCGGCGGGCTTAGGCAGGACGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCCTCCGAAGTGGGCATCTTCTTTTCATTTTATGACGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACttac                                                                                                                                         >  2:621642/1‑248 (MQ=255)
                                                                                                        ccAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAATATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATGTGCTTGCAACGTTGTTTCAGTTCATTTTCATCAGAAGTGGTCATCTTCTTTTCCTTTTATTCCGAAGGTTATGCGCCATTGTAAGAAGTTTCA‑TGATGTTCACTTTGATCCTGATGAGTTTTACTCCACTGACGcac                                                                                                                                     <  2:247554/250‑2 (MQ=255)
                                                                                                              cgcgGCGGGCTTCGGCAGGAAGCTGTGGGGAAGGCCGATTGTTTTCTGCTTCCCGTGCCAGAAACTGGCGCTTCTGCTCCGGGCTAAACACTGGGCTGGTGACAATTTGCTGGCAACGGTGTTGCAGGGCCTTTTCATCAGAAGTGGGCATCTTCTTTTACTTTTATGCCGAAGGGGATGCGCCATTGTAAGAAGATTCG‑TGTTGTTCCATTTGATCCTGATGAGTTTGCCACCAATGACGCATTCAtttt                                                                                                                              >  2:585157/1‑250 (MQ=255)
                                                                                                                    gggCGTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTg                                                                                                                        >  2:84906/1‑251 (MQ=255)
                                                                                                                      gctccGGCAGGAAGCTGTGGGTAAGGCCGATTGGTTGCTGCGTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGACATCTTCTTTTCCTTTTATGCCGAAGGGGATGCGCCATTGTAAGAAGGTTCG‑TTATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGcattcagatcggaagagc                                                                                                                      <  2:7152/247‑13 (MQ=255)
                                                                                                                       cTTCGGCAGGAAGCTGTGTGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATGAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGCGTTTGCCACCACTGACGCATTCATTTGAAAGTGAat                                                                                                                     >  1:264486/1‑251 (MQ=255)
                                                                                                                          cGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAattat                                                                                                                  >  2:602927/1‑251 (MQ=255)
                                                                                                                            gCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTAttt                                                                                                                <  1:21040/251‑1 (MQ=255)
                                                                                                                              aGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCCGTGCCAGAACATGGCGCTTCTGCTCCGGGCTAAGCACTGCGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTACGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGCAAGTGACTTATTTg                                                                                                               >  1:286090/1‑250 (MQ=255)
                                                                                                                              aGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGACAGAAAATGGCGATTCTGCTACGGACTAAGCACTGGGCTGGTGACAATATGCTGGCAACGTTGTTGCAGAGAATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGAAGAAGGTGATGCGCCATTGAACAAAGTTTCG‑TGATGTTCTATTTGATCCTGCTGAGTGTGCCCCCACTGACGCATCCATTGGAAAGTGAATTATTTg                                                                                                               >  2:299227/1‑250 (MQ=255)
                                                                                                                                  aGCTGTGAGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGCAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAAcc                                                                                                           >  2:578898/1‑250 (MQ=255)
                                                                                                                                       tgGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAAcc                                                                                                           >  1:358858/1‑245 (MQ=255)
                                                                                                                                       tgGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAAcc                                                                                                           >  2:355232/1‑245 (MQ=255)
                                                                                                                                                        tttCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCaa                                                                                     <  1:581960/250‑1 (MQ=255)
                                                                                                                                                           cggCTTCCAGTTCCATAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACTTTGTTTCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTGTATGGCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTTATCCTGATGAGGTTGCCAGCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTATAGTGATGCAAACt                                                                                  <  2:643834/248‑1 (MQ=255)
                                                                                                                                                               ttCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACt                                                                                                                                            <  1:416377/188‑1 (MQ=255)
                                                                                                                                                               ttCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACt                                                                                                                                            >  2:416377/1‑188 (MQ=255)
                                                                                                                                                                   aGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAg                                                                         <  1:578898/251‑1 (MQ=255)
                                                                                                                                                                                      tCTGCTCCGGGCTAAGCACTGGGCGGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGAtg                                                       <  2:20898/250‑1 (MQ=255)
                                                                                                                                                                                         gCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCttt                                                                                                                                                                                                                            >  1:298734/1‑83 (MQ=255)
                                                                                                                                                                                         gCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCttt                                                                                                                                                                                                                            <  2:298734/83‑1 (MQ=255)
                                                                                                                                                                                         gCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGAtgtata                                                   <  1:258327/251‑4 (MQ=255)
                                                                                                                                                                                         gCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGAtgtat                                                    <  1:98734/250‑3 (MQ=255)
                                                                                                                                                                                          cTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATtg                                                            <  1:415650/241‑1 (MQ=255)
                                                                                                                                                                                          cTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATtg                                                            >  2:415650/1‑241 (MQ=255)
                                                                                                                                                                                                gCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATGAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGCGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAgt                                            >  1:214486/1‑251 (MQ=255)
                                                                                                                                                                                                           ggCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCgg                                   >  1:597761/1‑249 (MQ=255)
                                                                                                                                                                                                           ggCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAg                                 <  2:351250/251‑1 (MQ=255)
                                                                                                                                                                                                            gCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGt                                <  2:412322/251‑1 (MQ=255)
                                                                                                                                                                                                                gTGACAATGTGCTGGGAACGTTGTTGCAGTGCTGTTTTATGAGAAGTTGGCATCTTCTCTTCCTTTTATGCCGGAGGTGATGCGCCAGTGTAAGAAGATTCG‑TGATGTTCACGTTGATCCTGATGCTTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCATATCGCATCACAGTGATGCAAACTTGTCAGTAGAGTTCCTTGATTGTGAAGTGTATCGAAGTGTGTTGCGGATTAGAt                            <  2:284513/251‑1 (MQ=255)
                                                                                                                                                                                                                                aaCGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTTAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACaa            <  1:673577/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gttgCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCGTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCg       <  1:585157/250‑1 (MQ=255)
                                                                                                                                                                                                                                          cAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAgg  >  1:371801/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             
AGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGC‑GGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGCGCGCGGCGGGCTTCGGCAGGAAGCTGTGGGTAAGGCAGATTGTTTTCTGCTTCCAGTGCCAGAAAATGGCGCTTCTGCTCCGGGCTAAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATGAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCG‑TGATGTTCACTTTGATCCTGATGCGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGG  >  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: ‑