Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,089,816 A→G 5.2% S14S (AGT→AGC pgaB ← poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine N‑deacetylase and beta‑1,6 glycoside hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,089,8160AG 5.2% 363.1 / 6.8 117S14S (AGT→AGCpgaBpoly‑beta‑1,6‑N‑acetyl‑D‑glucosamine N‑deacetylase and beta‑1,6 glycoside hydrolase
Reads supporting (aligned to +/- strand):  ref base A (52/58);  new base G (2/4);  total (55/62)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

Only 100 of 121 total aligned reads displayed.

TGGTCGGCAGCTTCGTCTTCAACGTTATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGTTTTTCCCAGCGTAGCGTTGCGCCAGCATCAATC  >  NC_000913/1089667‑1089965
                                                                                                                                                     |                                                                                                                                                     
tGGTCGGCAGCTTCGTCTTCAACGTTATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATtata                                                                                                                                                       >  1:437667/1‑150 (MQ=255)
  gTCGGCAGCTTCGTCTTAAACGTTATGCCATGAAATCGCTAAAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTAGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATAcg                                                                                                                                                     >  2:662802/1‑149 (MQ=255)
     ggCAGCTTCGTCTTCAACGTTATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTAGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCAc                                                                                                                                                  >  2:395982/1‑150 (MQ=255)
        aGCTTCGTCTTCAACGTTATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGTATAATTATGCTCAccgc                                                                                                                                               <  1:228272/150‑3 (MQ=255)
        aGCTTCGTCTTCAACGTTATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGTATAATTATGCTCAccgc                                                                                                                                               <  2:1398169/150‑3 (MQ=255)
         gCTTCGTCTTCAACGTTATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGTATAATTATGCTCAccgcg                                                                                                                                              <  1:1170895/150‑4 (MQ=255)
         gCTTCGTCTTCAACGTTATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGc                                                                                                                                              >  2:1490018/1‑150 (MQ=255)
          cTTCGTCTTCAACGTTATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCa                                                                                                                                             <  1:948901/150‑1 (MQ=255)
              gTCTTCAACGTTATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAg                                                                                                                                         <  2:1197235/150‑1 (MQ=255)
              gTCTTCAACGTTATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAg                                                                                                                                         <  2:1081044/150‑1 (MQ=255)
              gTCTTCAACGTTATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAg                                                                                                                                         <  2:501772/150‑1 (MQ=255)
              gTCTTCAACGTTATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAg                                                                                                                                         >  2:1029551/1‑150 (MQ=255)
                 ttCAACGTTATGCCAAGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGgag                                                                                                                                      >  2:357413/1‑150 (MQ=255)
                 ttCAACGTTAAGCCATGAAATCGCTACAAAACAATTATGCGGCCACCGAACCTCGGCGACAAAAGATTCGCGATCCTGTGGCGCTATAAATGATGGTCTTGACTGGCTAATGCACGCGGTGAGCATATTTATCCTGACCAGCATCAGgag                                                                                                                                      >  2:499466/1‑150 (MQ=255)
                   cAACGTTATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATAc                                                                                                                                                      <  2:243559/132‑1 (MQ=255)
                   cAACGTTATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATAc                                                                                                                                                      >  1:243559/1‑132 (MQ=255)
                   cAACGTTATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGga                                                                                                                                       >  1:809217/1‑147 (MQ=255)
                   cAACGTTATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGga                                                                                                                                       <  2:809217/147‑1 (MQ=255)
                   cAACGTTATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGat                                                                                                                                    >  1:161191/1‑150 (MQ=255)
                    aaCGTTATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGata                                                                                                                                   <  1:1080077/150‑1 (MQ=255)
                        ttATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGATCGGCGAGTAAAGATTCGCGATCCTGTGTCTGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGAtattta                                                                                                                               <  1:542980/150‑1 (MQ=255)
                             ccATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTcc                                                                                                                          >  1:433071/1‑150 (MQ=255)
                                 gAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTa                                                                                                                      <  2:855540/150‑1 (MQ=255)
                                 gAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTa                                                                                                                      >  2:5952/1‑150 (MQ=255)
                                  aaaTCGCTACAAAACGATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTAc                                                                                                                     <  2:1316087/150‑1 (MQ=255)
                                  aaaTCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTAc                                                                                                                     <  1:1453005/150‑1 (MQ=255)
                                        cTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACa                                                                                                               >  1:1418050/1‑150 (MQ=255)
                                        cTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACa                                                                                                               <  1:732692/150‑1 (MQ=255)
                                             aaaCCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCCCCAGCATCAGGAGATATTTATTTCCATTACGTAACATAttt                                                                                                          >  2:1032681/1‑150 (MQ=255)
                                             aaaCCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATAttt                                                                                                          >  1:758495/1‑150 (MQ=255)
                                                ccATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGAACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATc                                                                                                       >  1:783252/1‑150 (MQ=255)
                                                  aTTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCt                                                                                                     >  1:1144806/1‑150 (MQ=255)
                                                      tGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTaaa                                                                                                 >  1:1141972/1‑150 (MQ=255)
                                                       gCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTaaaa                                                                                                >  2:1175352/1‑150 (MQ=255)
                                                        cGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTTAGCATAATTATACTCACCAGCATCAGGACATATTTATTTCCATTACGTAACATATTTATCCTTAAAAt                                                                                               <  2:1477525/150‑1 (MQ=255)
                                                        cGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAAt                                                                                               >  1:1499856/1‑150 (MQ=255)
                                                         ggCCACGGTTGCTCCGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATc                                                                                              >  2:177366/1‑150 (MQ=255)
                                                            cACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGa                                                                                           <  2:1144806/150‑1 (MQ=255)
                                                              cGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAAt                                                                                         <  2:739813/150‑1 (MQ=255)
                                                               ggTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCGCCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATg                                                                                        <  1:1181307/150‑1 (MQ=255)
                                                                gTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGt                                                                                       >  2:1220511/1‑150 (MQ=255)
                                                                  tGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCa                                                                                     <  2:433071/150‑1 (MQ=255)
                                                                     tCGCCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCatat                                                                                  >  2:404008/1‑150 (MQ=255)
                                                                      cGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGCCATATc                                                                                 >  1:80438/1‑150 (MQ=255)
                                                                       ggCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCg                                                                                >  2:822686/1‑150 (MQ=255)
                                                                       ggCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCg                                                                                <  1:1038973/150‑1 (MQ=255)
                                                                       ggCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCg                                                                                <  1:1490018/150‑1 (MQ=255)
                                                                        gCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGa                                                                               >  2:817731/1‑150 (MQ=255)
                                                                         cGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGaa                                                                              <  1:5952/150‑1 (MQ=255)
                                                                              aaaGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCaa                                                                         >  2:476599/1‑150 (MQ=255)
                                                                               aaGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAAc                                                                        <  2:692970/150‑1 (MQ=255)
                                                                                aGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACg                                                                       <  2:437667/150‑1 (MQ=255)
                                                                                aGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACg                                                                       <  1:478408/150‑1 (MQ=255)
                                                                                     cgcgATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGta                                                                     >  1:1221250/1‑147 (MQ=255)
                                                                                     cgcgATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGta                                                                     <  2:1221250/147‑1 (MQ=255)
                                                                                      gcgATGCTGTGGCTGTATAAATGATGTTCTTGACTGGCTAATGCACGCTGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAg                                                                 <  1:1382687/150‑1 (MQ=255)
                                                                                         aTCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATTCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTtg                                                              <  1:138627/150‑1 (MQ=255)
                                                                                         aTCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGATCATAATTATACTCACCAGCATCTGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATTAGTtg                                                              <  2:664237/150‑1 (MQ=255)
                                                                                         aTCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTtg                                                              <  2:281586/150‑1 (MQ=255)
                                                                                            cTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTgtgt                                                           <  2:873663/150‑1 (MQ=255)
                                                                                             tgtgGAGTTATAAATGATGTTCTTGACTTGCTAATGCACGCGTTGAGCATAATTATACTCACCAGCATCAGGATATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGtt                                                          <  1:600486/150‑1 (MQ=255)
                                                                                                ggCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTtctc                                                       <  2:559821/150‑1 (MQ=255)
                                                                                                 gCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTctct                                                      <  1:1115092/150‑1 (MQ=255)
                                                                                                 gCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTctct                                                      <  2:87112/150‑1 (MQ=255)
                                                                                                 gCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGCCATATCGAATTCAACGTATAAGTTGTGTTctct                                                      <  2:80438/150‑1 (MQ=255)
                                                                                                  cGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTTTTCTCTg                                                     <  2:161191/150‑1 (MQ=255)
                                                                                                  cGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTg                                                     <  2:1499856/150‑1 (MQ=255)
                                                                                                  cGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTg                                                     <  2:90874/150‑1 (MQ=255)
                                                                                                  cGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTg                                                     <  2:933358/150‑1 (MQ=255)
                                                                                                     tataAATTATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCa                                                  <  1:668524/150‑1 (MQ=255)
                                                                                                     tataAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCa                                                  >  1:1358304/1‑150 (MQ=255)
                                                                                                        aaaTGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCg                                               <  1:42452/150‑1 (MQ=255)
                                                                                                        aaaTGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCg                                               <  1:1124508/150‑1 (MQ=255)
                                                                                                           tgatgTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCa                                            >  1:1340176/1‑150 (MQ=255)
                                                                                                             atgTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATa                                          <  2:847515/150‑1 (MQ=255)
                                                                                                               gTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAg                                        >  1:272161/1‑150 (MQ=255)
                                                                                                               gTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAg                                        <  2:1358304/150‑1 (MQ=255)
                                                                                                                ttCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAgg                                       >  1:263902/1‑150 (MQ=255)
                                                                                                                 tCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGt                                      >  2:469861/1‑150 (MQ=255)
                                                                                                                    tGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGt                                   >  1:919844/1‑150 (MQ=255)
                                                                                                                       cTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGtttt                                >  2:191047/1‑150 (MQ=255)
                                                                                                                       cTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGtttt                                <  1:395982/150‑1 (MQ=255)
                                                                                                                         gtcTAATGCCAGAGGTGAGCATAATTATACTCCCAAGCTTCAGGAGATTTTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTAATATCGAATTAAACGTATAAGTTGTGGTATCTGTCACCGTCATAATGTCGTTTTTc                              >  2:1439126/3‑150 (MQ=255)
                                                                                                                             aaTGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGTTTTTCCCag                          <  1:585254/150‑1 (MQ=255)
                                                                                                                              aTGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGTTTTTCCCagc                         >  1:1042054/1‑150 (MQ=255)
                                                                                                                               tGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGTTTTTCCCagcg                        >  2:1516585/1‑150 (MQ=255)
                                                                                                                                gCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGTTTTTCCCagcgt                       >  2:1399379/1‑150 (MQ=255)
                                                                                                                                 cACGCGGTGATCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGTTTTTCCCagcgta                      <  1:662802/150‑1 (MQ=255)
                                                                                                                                      ggTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGTTTTTCCCAGCGTAGCGtt                 <  2:1042054/150‑1 (MQ=255)
                                                                                                                                        tGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGTTTTTCCCAGCGTAGCGTTgc               >  2:583731/1‑150 (MQ=255)
                                                                                                                                         taGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGTTTTTCCCAGCGTAGCGTTgcg              <  2:333645/149‑1 (MQ=255)
                                                                                                                                          aGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGTTTTTCCCa                           >  1:805516/1‑136 (MQ=255)
                                                                                                                                          aGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGTTTTTCCCAGCGTAGCGTTgcgc             >  1:232883/1‑150 (MQ=255)
                                                                                                                                           gCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGTTTTTCCCAGCGTAGCGTTGCGcc            >  1:873334/1‑150 (MQ=255)
                                                                                                                                             aTAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAGTGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGTTTTTCCCAGCGTAGCGTTGCGCCAg          <  2:153087/150‑1 (MQ=255)
                                                                                                                                               aaTTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGTTTTTCCCAGCGTAGCGTTGCGCCAGCa        <  1:546353/150‑1 (MQ=255)
                                                                                                                                               aaTTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGTTTTTCCCAGCGTAGCGTTGCGCCAGCa        <  1:476599/150‑1 (MQ=255)
                                                                                                                                                   ataCTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGTTTTTCCCAGCGTAGCGTTGCGCCAGCatcaa    >  1:1014655/1‑150 (MQ=255)
                                                                                                                                                    taCTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGTTTTTCCCAGCGTAGCGTTGCGCCAGCatcaat   >  1:776291/1‑150 (MQ=255)
                                                                                                                                                     aCTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGTTTTTCCCAGCGTAGCGTTGCGCCAGCatcaatc  <  1:1527208/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGTCGGCAGCTTCGTCTTCAACGTTATGCCATGAAATCGCTACAAAACCATTATGCGGCCACGGTTGCTCGGCGAGTAAAGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCATATCGAATTCAACGTATAAGTTGTGTTCTCTGTCACCGTCATAAGGTCGTTTTTCCCAGCGTAGCGTTGCGCCAGCATCAATC  >  NC_000913/1089667‑1089965

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