Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 6,050,487 C→G 6.0% P450R (CCG→CGG)  spoT → bifunctional (p)ppGpp synthase/hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154516,050,4870CG 6.0% 234.9 / 7.4 67P450R (CCG→CGG) spoTbifunctional (p)ppGpp synthase/hydrolase
Reads supporting (aligned to +/- strand):  ref base C (25/38);  new base G (2/2);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-02

AGCGTGAAGATCGACCTGTTCCCGGACGAGGTCTACGTGTTCACGCCCAAAGGCCGGATCATGGAGTTGCCCAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTCGACAGCAGCCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAG  >  AE015451/6050273‑6050720
                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
agCGTGAAGATCGACCTGTTCCCGGACGAGGTCTACGTGTTCACGCCCAAAGGCCGGATCATGGAGTTGCCCAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAAcc                                                                                                                                                                                                                          <  1:214639/232‑1 (MQ=255)
   gTGAAGATCGACCTGTTCCCGGACGAGGTCTACGTGTTCACGCCCAAAGGCCGGATCATGGAGTTGCCCAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATg                                                                                                                                                                                                                   >  1:1042024/1‑236 (MQ=255)
             aCCTGTTCCCGGACTAGGTCTACGTGTTCACGCCCAAAGGCCGGATCATGGAGTTGCCCAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACt                                                                                                                                                                                                           <  1:559440/234‑1 (MQ=255)
                     ccGGACGAGGTCTACGTGTTCACGCCCAAAGGCCGGATCATGGAGTTGCCCAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGtct                                                                                                                                                                                                  >  2:865147/1‑235 (MQ=255)
                                 tACGTGTTCACGCCCAAAGGCCGGATCATGGAGTTGCCCAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCACACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGgcacg                                                                                                                                                                                     >  1:89142/1‑236 (MQ=255)
                                            gCCCAAAGGCCGGATCATGGAGTTGCCCAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAgg                                                                                                                                                                                         >  1:20381/1‑221 (MQ=255)
                                            gCCCAAAGGCCGGATCATGGAGTTGCCCAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATc                                                                                                                                                                           >  1:717280/1‑235 (MQ=255)
                                                         aTCATGGAGTTGCCCAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCa                                                                                                                                                              >  1:915295/1‑235 (MQ=255)
                                                            aTGGAGTTGCCCAAAGTCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCACCCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCg                                                                                                                                                                                                                        <  2:348387/174‑1 (MQ=255)
                                                             tGGAGTTGCCCAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTgga                                                                                                                                                                                                    >  2:126498/1‑192 (MQ=255)
                                                              ggAGTTGCCCAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAg                                                                                                                                                                                          <  1:491650/202‑1 (MQ=255)
                                                                       cAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGCAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGgcacgc                                                                                                                                                                                    >  1:351011/1‑199 (MQ=255)
                                                                        aaaGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTc                                                                                                                                              >  1:22549/1‑236 (MQ=255)
                                                                         aaGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTc                                                                                                                                              <  1:625878/235‑1 (MQ=255)
                                                                          aGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCt                                                                                                                                               <  1:90613/233‑1 (MQ=255)
                                                                                 gcGGCCGTCGCCTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACa                                                                                                                                                       <  1:315814/217‑1 (MQ=255)
                                                                                  cGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAgcgc                                                                                                                                                   >  1:585450/1‑221 (MQ=255)
                                                                                  cGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCAt                                                                                                                                     >  1:1064771/1‑235 (MQ=255)
                                                                                             ttCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCg                                                                                                                          <  1:523121/235‑1 (MQ=255)
                                                                                              tCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCGCGGGCGCTCGGCCAAACCCGGCATg                                                                                                                                                                                                                   <  2:665770/145‑1 (MQ=255)
                                                                                                gCCCACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCg                                                                                                                          <  1:189844/232‑1 (MQ=255)
                                                                                                 ccTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCACCCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgtg                                                                                                                                                                                                      <  2:749087/155‑1 (MQ=255)
                                                                                                 ccTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAgcgc                                                                                                                     <  1:244856/236‑1 (MQ=255)
                                                                                                  cTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCaa                                                                                                                                                                                           >  1:1077035/1‑165 (MQ=255)
                                                                                                   caCGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCg                                                                                                                          <  1:102550/228‑1 (MQ=255)
                                                                                                        ggTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGGGAAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTc                                                                                                                                              <  2:1034330/204‑1 (MQ=255)
                                                                                                         gTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCctgctg                                                                                                               <  2:931355/234‑1 (MQ=255)
                                                                                                          tCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAg                                                                                                                                          <  2:505071/206‑1 (MQ=255)
                                                                                                            cacaCCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAg                                                                                                                        <  2:429725/222‑1 (MQ=255)
                                                                                                                       cGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTg                                                                                               >  1:1060817/1‑236 (MQ=255)
                                                                                                                       cGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCCCAACTTTgggg                                                                                                                                                                                                     >  2:417723/1‑134 (MQ=255)
                                                                                                                       cGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCCCGACAGTGGAAATCGTCAGCGCCCCg                                                                                                                                                                                                                                          >  2:338453/1‑97 (MQ=255)
                                                                                                                            aCAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCg                                                                                                                          <  2:667684/204‑1 (MQ=255)
                                                                                                                                ttGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGcacgcacg                                                                                                                                                                                 >  2:418454/1‑145 (MQ=255)
                                                                                                                                   cATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGaa                                                                                                                                                                               >  2:648538/1‑144 (MQ=255)
                                                                                                                                    aTCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCg                                                                                                                          <  1:78410/196‑1 (MQ=255)
                                                                                                                                       gCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACa                                                                                                           <  1:706166/208‑1 (MQ=255)
                                                                                                                                           gCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTCGa                                                                                        <  2:903620/223‑1 (MQ=255)
                                                                                                                                                   aaCCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTCGACAGCAGCCTGGAGAAAATcc                                                                    <  1:865150/235‑1 (MQ=255)
                                                                                                                                                    aCCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCg                                                                                                                          <  2:820913/180‑1 (MQ=255)
                                                                                                                                                     ccgccgCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTCGa                                                                                        <  2:982474/213‑1 (MQ=255)
                                                                                                                                                      cgccgcCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTCGACAGCAGCCTGGAGAAAATCCCCCa                                                                <  1:746316/236‑1 (MQ=255)
                                                                                                                                                            cTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACCCGCTCAAGc                                                                                                                                                           >  2:797421/1‑139 (MQ=255)
                                                                                                                                                                        tCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTCGACAGCAGCCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTa                                               <  1:1021960/235‑1 (MQ=255)
                                                                                                                                                                                gCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTCGACAGCAGCCTGGAGAAAATCCCCCAGGAACGc                                                         <  2:1064768/217‑1 (MQ=255)
                                                                                                                                                                                 cTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTCGACAGCAGCCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCgc                                     >  1:264448/1‑235 (MQ=255)
                                                                                                                                                                                     aaaGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACtgtgtg                                                                                                                                                                                                      <  1:260596/71‑1 (MQ=255)
                                                                                                                                                                                     aaaGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACt                                                                                                                                                                                                           >  2:1063267/1‑66 (MQ=255)
                                                                                                                                                                                     aaaGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACt                                                                                                                                                                                                           <  1:1063270/66‑1 (MQ=255)
                                                                                                                                                                                     aaaGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgtg                                                                                                                                                                                                      >  2:260593/1‑71 (MQ=255)
                                                                                                                                                                                     aaaGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAg                                                                                                                                                      >  2:397081/1‑119 (MQ=255)
                                                                                                                                                                                      aaGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTCGACAGCAGCCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTAcc                                <  1:1072454/236‑1 (MQ=255)
                                                                                                                                                                                       aGCGGCTGGACAGGGGAAATCGTCAGCGCCGCGGGCGCGCGGCCAAACCCGGCATGGCTCAACTGTGTGGTCTCGGGCGAGGCACGCACGGATATCCGCCa                                                                                                                                                                      <  2:18300/101‑1 (MQ=255)
                                                                                                                                                                                       aGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTCGAcagc                                                                                    <  2:717277/183‑1 (MQ=255)
                                                                                                                                                                                            cTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATc                                                                                                                                    >  1:865715/1‑130 (MQ=255)
                                                                                                                                                                                             tCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTCGACAGCAGCCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGa                         <  1:567621/236‑1 (MQ=255)
                                                                                                                                                                                               gACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTCGACAGCAGCCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAg                        >  1:753040/1‑235 (MQ=255)
                                                                                                                                                                                                aCAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTc                                                                                                                                              <  2:790559/116‑1 (MQ=255)
                                                                                                                                                                                                          tCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTCGACAGCAGCCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGa                         <  1:8197/223‑1 (MQ=255)
                                                                                                                                                                                                            gTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTCGACAGCAGCCTGGAGaaaa                                                                       <  2:521271/175‑1 (MQ=255)
                                                                                                                                                                                                             cgaGCGCCCGGGGCGCTCGGCCAAACCCGGCATGGCTCAACt                                                                                                                                                                                                           <  2:943462/40‑1 (MQ=255)
                                                                                                                                                                                                             cgaGCGCCCGGGGCGCTCGGCCAAACCCGGCATGGCTCAACt                                                                                                                                                                                                           >  1:943465/3‑42 (MQ=255)
                                                                                                                                                                                                             cgaGCGCCCGGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCAcgc                                                                                                                                                                   >  2:686726/3‑82 (MQ=255)
                                                                                                                                                                                                             cgaGCGCCCGGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCAcgc                                                                                                                                                                   <  1:686729/80‑1 (MQ=255)
                                                                                                                                                                                                               aGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTCGACAGCAGCCTGGAGAAAATCCCCCAGGAACg                                                          <  2:1060814/185‑1 (MQ=255)
                                                                                                                                                                                                                 cgcCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGGGAGgca                                                                                                                                                                                       >  2:490452/1‑58 (MQ=255)
                                                                                                                                                                                                                 cgcCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTCGACAGCAGCCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCt                            <  1:345949/213‑1 (MQ=255)
                                                                                                                                                                                                                    cccGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTCGACAGCAGCCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAg  <  1:14019/236‑1 (MQ=255)
                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
AGCGTGAAGATCGACCTGTTCCCGGACGAGGTCTACGTGTTCACGCCCAAAGGCCGGATCATGGAGTTGCCCAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACTGGCTTCGACAGCAGCCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAG  >  AE015451/6050273‑6050720

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