Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,963,299 Δ1 bp 9.5% coding (931/1485 nt) gpp ← guanosine pentaphosphatase/exopolyphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,963,2990G. 9.5% 307.4 / 27.7 95coding (931/1485 nt)gppguanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base G (50/36);  new base . (3/6);  total (53/42)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

AGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCA  >  NC_000913/3963155‑3963447
                                                                                                                                                 |                                                                                                                                                    
aGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTg                                                                                                                                                  >  4:670570/1‑149 (MQ=255)
 ggCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGc                                                                                                                                                 <  4:550194/149‑1 (MQ=255)
  gcgcTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTTGTTTTCCACCTGATCGAATAAGTTAGCCGCAACTCTGGCTACGCGCTGTGCCTGATCAATATCAATCATAAAGTGGCGTTGAATATTACGCAG‑CGTACGGCTGCg                                                                                                                                                >  3:773157/1‑149 (MQ=255)
     cTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCGCCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCCGTAC‑GCTGCgtaa                                                                                                                                             <  4:209816/149‑4 (MQ=255)
      tgatgaGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAata                                                                                                                                            >  4:818068/1‑149 (MQ=255)
      tgatgaGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAata                                                                                                                                            >  2:662320/1‑149 (MQ=255)
       gatgaGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCCGTAC‑GCTGCgtaata                                                                                                                                           <  4:20004/149‑6 (MQ=255)
       gatgaGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCCGTAC‑GCTGCgtaata                                                                                                                                           <  1:499690/149‑6 (MQ=255)
        atgaGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTTTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCCGTAC‑GCTGCgtaata                                                                                                                                           <  4:283026/148‑6 (MQ=255)
        atgaGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATc                                                                                                                                          <  1:545042/149‑1 (MQ=255)
         tgaGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATcc                                                                                                                                         <  2:933920/149‑1 (MQ=255)
          gaGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGAAG‑CGTACGGCTGCGAATATCCt                                                                                                                                        >  4:761812/1‑149 (MQ=255)
           aGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTg                                                                                                                                       >  1:150126/1‑149 (MQ=255)
            gCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGc                                                                                                                                      <  1:608485/149‑1 (MQ=255)
             cAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCt                                                                                                                                     <  1:69057/149‑1 (MQ=255)
              aaaTCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTc                                                                                                                                    >  3:894428/1‑149 (MQ=255)
                      gCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCa                                                                                                                            >  4:26333/1‑149 (MQ=255)
                          aTTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAGTATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATg                                                                                                                        >  4:664669/1‑149 (MQ=255)
                           ttGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGc                                                                                                                       <  1:422466/149‑1 (MQ=255)
                             gCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCaa                                                                                                                     <  2:343532/149‑1 (MQ=255)
                              cTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAAc                                                                                                                    >  1:356805/1‑149 (MQ=255)
                               ttCAAGATGCCGTTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACa                                                                                                                   >  3:541322/1‑149 (MQ=255)
                                  aaGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGc                                                                                                                >  3:735974/1‑149 (MQ=255)
                                      tGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTa                                                                                                            <  1:877527/149‑1 (MQ=255)
                                       gCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAg                                                                                                           <  1:454060/149‑1 (MQ=255)
                                        ccATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGa                                                                                                          >  4:615922/1‑149 (MQ=255)
                                          aTTCGTTTTCCACCTGATCGACGACGTTAATCGCAGCTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGCCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCGCGACGGCCAGATGCAACATGCCGTAGAcc                                                                                                        >  4:775211/1‑149 (MQ=255)
                                          aTTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGAcc                                                                                                        >  2:816298/1‑149 (MQ=255)
                                          aTTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGAcc                                                                                                        >  1:422142/1‑149 (MQ=255)
                                               ttttCCACCTGATCGTAGAAGTTAGCCGCGCCTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATTAAGCGACGTTGAATATGACGGAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCGCGATGCAACATGCCGTAGACCAGGcc                                                                                                   >  4:775093/1‑149 (MQ=255)
                                               ttttCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGTcc                                                                                                   >  2:649944/1‑149 (MQ=255)
                                                      ccTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGc                                                                                            >  4:653069/1‑149 (MQ=255)
                                                      ccTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATAGCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGc                                                                                            <  1:873470/149‑1 (MQ=255)
                                                       cTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCa                                                                                           >  4:474944/1‑149 (MQ=255)
                                                         gATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGAGCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAgc                                                                                         <  1:662320/149‑1 (MQ=255)
                                                             gaagaaGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCAc                                                                                     >  3:432319/1‑149 (MQ=255)
                                                                   gTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCAccgcccg                                                                               >  2:880991/1‑149 (MQ=255)
                                                                    ttAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGc                                                                                                    >  3:593639/1‑127 (MQ=255)
                                                                    ttAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGc                                                                                                    <  4:593639/127‑1 (MQ=255)
                                                                    ttAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCACAccgcccgc                                                                              >  4:791130/1‑149 (MQ=255)
                                                                        ccGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAgg                                                                          <  2:356805/149‑1 (MQ=255)
                                                                         cGCAACCTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAAACCAGGCCTTCACGCAGCGCACCGCCCGCCAggg                                                                         <  3:615922/149‑1 (MQ=255)
                                                                          gCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGt                                                                        >  3:290655/1‑149 (MQ=255)
                                                                          gCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGt                                                                        >  3:35476/1‑149 (MQ=255)
                                                                            aaCTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCa                                                                      <  4:735974/149‑1 (MQ=255)
                                                                              cTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATa                                                                    <  4:894428/149‑1 (MQ=255)
                                                                               tttGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGa                                                                                                          >  4:129058/1‑110 (MQ=255)
                                                                               tttGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGa                                                                                                          <  3:129058/110‑1 (MQ=255)
                                                                               tttGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATac                                                                   >  2:220249/1‑149 (MQ=255)
                                                                                  ggCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACt                                                                >  4:787911/1‑149 (MQ=255)
                                                                                     tACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGaa                                                             <  1:576469/149‑1 (MQ=255)
                                                                                      aCGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAat                                                            >  3:409778/1‑149 (MQ=255)
                                                                                         cgcTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATAtt                                                         <  1:816298/149‑1 (MQ=255)
                                                                                          gcTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTc                                                        <  1:220249/149‑1 (MQ=255)
                                                                                           cTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACTCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCa                                                       <  3:791130/149‑1 (MQ=255)
                                                                                           cTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAata                                                                                                                                            >  1:407387/1‑64 (MQ=255)
                                                                                           cTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAata                                                                                                                                            <  2:407387/64‑1 (MQ=255)
                                                                                            tgtgCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAg                                                      <  2:150126/149‑1 (MQ=255)
                                                                                              tgCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGtt                                                    <  4:35476/149‑1 (MQ=255)
                                                                                                ccTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCg                                                  >  1:725726/1‑149 (MQ=255)
                                                                                                   gATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTa                                               >  2:420575/1‑149 (MQ=255)
                                                                                                    aTCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTaa                                              >  3:696471/1‑149 (MQ=255)
                                                                                                    aTCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTaa                                              >  1:632749/1‑149 (MQ=255)
                                                                                                          atCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATgg                                        >  1:387597/1‑149 (MQ=255)
                                                                                                             gATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGa                                     >  4:848648/1‑149 (MQ=255)
                                                                                                              aTCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGAt                                    >  2:685887/1‑149 (MQ=255)
                                                                                                               tCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCa                                                       >  4:70220/1‑129 (MQ=255)
                                                                                                               tCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATc                                   >  2:46947/1‑149 (MQ=255)
                                                                                                               tCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCACACACTGAATATTCa                                                       <  3:70220/129‑1 (MQ=255)
                                                                                                                 aTAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGc                                       <  2:613990/143‑1 (MQ=255)
                                                                                                                 aTAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGc                                       >  1:613990/1‑143 (MQ=255)
                                                                                                                 aTAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAg                                 >  2:580939/1‑149 (MQ=255)
                                                                                                                 aTAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAg                                 >  3:511448/1‑149 (MQ=255)
                                                                                                                  tAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAgg                                >  2:213134/1‑149 (MQ=255)
                                                                                                                     aGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATc                             <  1:873745/149‑1 (MQ=255)
                                                                                                                     aGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATc                             >  3:555442/1‑149 (MQ=255)
                                                                                                                        ggcgTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGcc                          <  4:541322/149‑1 (MQ=255)
                                                                                                                         gcgTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCa                         <  1:567537/149‑1 (MQ=255)
                                                                                                                          cgTTGAATCTTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCATCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCTGTAAAAATGGCGATCCTGATCGCCat                        <  3:775211/149‑2 (MQ=255)
                                                                                                                           gTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCatt                       <  1:46947/149‑3 (MQ=255)
                                                                                                                            ttGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCatt                       <  1:851535/148‑3 (MQ=255)
                                                                                                                            ttGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCatt                       >  2:851535/1‑146 (MQ=255)
                                                                                                                            ttGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGAc                      >  4:499824/1‑149 (MQ=255)
                                                                                                                            ttGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGAc                      <  4:511448/149‑1 (MQ=255)
                                                                                                                                  tattcGCAGCCGTAC‑GCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATg                                         >  3:267516/5‑124 (MQ=255)
                                                                                                                                  tattcGCAGCCGTAC‑GCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATg                                         <  4:267516/120‑1 (MQ=255)
                                                                                                                                   ttACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCgg               >  3:273405/1‑149 (MQ=255)
                                                                                                                                   ttACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCgg               >  4:47480/1‑149 (MQ=255)
                                                                                                                                    tACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCggg              <  1:580939/149‑1 (MQ=255)
                                                                                                                                     tcGCAGCCGTAC‑GCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGc                                       >  1:27688/2‑123 (MQ=255)
                                                                                                                                     tcGCAGCCGTAC‑GCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGc                                       <  2:27688/122‑1 (MQ=255)
                                                                                                                                     tcGCAGCCGTAC‑GCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGa             >  3:263205/2‑149 (MQ=255)
                                                                                                                                      cGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAATTGGCGATCAGGATCGCCAGACCACTCGGGaa            <  2:4300/149‑1 (MQ=255)
                                                                                                                                      cGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGaa            >  1:887152/1‑149 (MQ=255)
                                                                                                                                        cAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAaca          <  1:420575/149‑1 (MQ=255)
                                                                                                                                                 ggCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCa  >  4:633532/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCA  >  NC_000913/3963155‑3963447

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑