Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,851,242 G→T E96* (GAG→TAG)  relA → ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,851,2420GT98.3% 1043.4 / 3.2 297E96* (GAG→TAG) relAATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base T (136/156);  total (139/158)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

Only 100 of 297 total aligned reads displayed.

GTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTGGGCTCGCAGGCGCAGGTA  >  NC_002947/1851095‑1851385
                                                                                                                                                   |                                                                                                                                               
gTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAg                                                                                                                                               >  2:4609491/1‑150 (MQ=255)
gTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAg                                                                                                                                               <  1:6342719/150‑1 (MQ=255)
gTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAg                                                                                                                                               >  1:4830756/1‑150 (MQ=255)
 tAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAgg                                                                                                                                              <  2:5465892/150‑1 (MQ=255)
   gagaAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGc                                                                                                                                            <  2:5203036/150‑1 (MQ=255)
     gaAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCaa                                                                                                                                          >  2:1956091/1‑150 (MQ=255)
      aaaaaGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAg                                                                                                                                         >  2:2874746/1‑150 (MQ=255)
      aaaaaGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAg                                                                                                                                         <  2:3754063/150‑1 (MQ=255)
      aaaaaGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAg                                                                                                                                         >  2:2126212/1‑150 (MQ=255)
       aaaaGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAgg                                                                                                                                        >  1:930447/1‑150 (MQ=255)
       aaaaGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAgg                                                                                                                                        >  2:5281431/1‑150 (MQ=255)
        aaaGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGt                                                                                                                                       >  2:1555342/1‑150 (MQ=255)
             gCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCAccc                                                                                                                                  >  1:1961482/1‑150 (MQ=255)
                     gCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAg                                                                                                                          <  2:1748779/150‑1 (MQ=255)
                     gCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAg                                                                                                                          <  1:4275392/150‑1 (MQ=255)
                     gCCAAGCAGTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAg                                                                                                                          <  2:4199044/150‑1 (MQ=255)
                      ccAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAgg                                                                                                                         <  2:2903428/150‑1 (MQ=255)
                        aaGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTc                                                                                                                       >  1:4797396/1‑150 (MQ=255)
                          gCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTcag                                                                                                                     >  2:2224584/1‑150 (MQ=255)
                          gCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTcag                                                                                                                     >  2:2259201/1‑150 (MQ=255)
                             ttCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTcagcca                                                                                                                  <  1:961830/150‑1 (MQ=255)
                                  gggCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGt                                                                                                             <  2:5847906/150‑1 (MQ=255)
                                     cGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTg                                                                                                          <  1:1555381/150‑1 (MQ=255)
                                     cGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTg                                                                                                          <  2:5120395/150‑1 (MQ=255)
                                      ggATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTgg                                                                                                         >  1:3983632/1‑150 (MQ=255)
                                      ggATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTgg                                                                                                         <  2:2069253/150‑1 (MQ=255)
                                          ggTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCg                                                                                                     <  1:6418717/150‑1 (MQ=255)
                                          ggTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCg                                                                                                     <  2:2960503/150‑1 (MQ=255)
                                           gTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCgg                                                                                                    >  2:2048099/1‑150 (MQ=255)
                                           gTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCgg                                                                                                    >  2:3221034/1‑150 (MQ=255)
                                           gTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCgg                                                                                                    >  2:5910919/1‑150 (MQ=255)
                                             aCGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCggtg                                                                                                  <  2:4526711/150‑1 (MQ=255)
                                               gTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCggtggt                                                                                                >  1:3853033/1‑150 (MQ=255)
                                               gTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCggtggt                                                                                                >  1:4897357/1‑150 (MQ=255)
                                               gTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTAGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCggtggt                                                                                                <  1:6541209/150‑1 (MQ=255)
                                                  cAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTc                                                                                             >  1:819451/1‑150 (MQ=255)
                                                  cAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTc                                                                                             <  1:5322609/150‑1 (MQ=255)
                                                   aGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTcc                                                                                            <  2:133138/150‑1 (MQ=255)
                                                    gCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCa                                                                                           >  1:3395011/1‑150 (MQ=255)
                                                    gCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCa                                                                                           >  1:3791727/1‑150 (MQ=255)
                                                      ttCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGttt                                                                                                           <  1:6131103/132‑1 (MQ=255)
                                                       tCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGc                                                                                        <  2:4282948/150‑1 (MQ=255)
                                                       tCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGc                                                                                        <  1:4139130/150‑1 (MQ=255)
                                                          aggcaggcCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGa                                                                                     <  1:2107123/150‑1 (MQ=255)
                                                            gcctgcCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATc                                                                                   >  2:3051025/5‑150 (MQ=255)
                                                            gcaggcCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGa                                                                                                                           >  2:3497922/1‑110 (MQ=255)
                                                            gcaggcCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATc                                                                                   >  1:306853/1‑150 (MQ=255)
                                                              ctgcCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCCGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGa                                                                                 >  1:2198504/3‑150 (MQ=255)
                                                              ctgcCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGAGGTTTGGCCCGGTGGTGTCCAAGCTGATCGa                                                                                 >  1:5492481/3‑150 (MQ=255)
                                                               ggcCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGAc                                                                                <  2:1124420/150‑1 (MQ=255)
                                                                gcCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACg                                                                               >  1:4303016/1‑150 (MQ=255)
                                                                  cTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGgt                                                                             <  2:2375580/150‑1 (MQ=255)
                                                                  cTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGgt                                                                             >  1:5195699/1‑150 (MQ=255)
                                                                   tGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGgtg                                                                            >  2:6327097/1‑150 (MQ=255)
                                                                         tCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGtgctgc                                                                      <  1:1512291/150‑1 (MQ=255)
                                                                           gCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTgcgc                                                                    >  1:458168/1‑150 (MQ=255)
                                                                             cGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCAt                                                                  >  1:1386872/1‑150 (MQ=255)
                                                                              gAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTTCTGCGCATg                                                                 >  2:325477/1‑150 (MQ=255)
                                                                               aaaTTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATgg                                                                <  2:1018360/150‑1 (MQ=255)
                                                                               aaaTTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATgg                                                                <  1:2354874/150‑1 (MQ=255)
                                                                                      ggCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCAt                                                         <  1:2870515/150‑1 (MQ=255)
                                                                                       gCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATc                                                        >  1:5307056/1‑150 (MQ=255)
                                                                                       gCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATc                                                        <  2:2646994/150‑1 (MQ=255)
                                                                                          gACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGt                                                     <  2:2753107/150‑1 (MQ=255)
                                                                                             cTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGcc                                                  >  1:1010501/1‑150 (MQ=255)
                                                                                                aaGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGc                                               <  1:6458816/150‑1 (MQ=255)
                                                                                                aaGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGc                                               <  1:5396892/150‑1 (MQ=255)
                                                                                                aaGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGc                                               <  1:146894/150‑1 (MQ=255)
                                                                                                    tCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCa                                           <  1:3484057/150‑1 (MQ=255)
                                                                                                    tCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCa                                           >  1:786489/1‑150 (MQ=255)
                                                                                                     cGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAg                                          <  1:3646177/150‑1 (MQ=255)
                                                                                                      gACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGc                                         >  2:1370065/1‑150 (MQ=255)
                                                                                                        ccAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGATGCGCATGGCCGCCATCAGTGCCAGCCTCAGccc                                       >  1:1035710/1‑150 (MQ=255)
                                                                                                          aGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCAc                                     >  2:1653040/1‑150 (MQ=255)
                                                                                                           ggACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACg                                    <  2:6467542/150‑1 (MQ=255)
                                                                                                            gACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGa                                   >  1:3567693/1‑150 (MQ=255)
                                                                                                             aCTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGAc                                  >  1:5779576/1‑150 (MQ=255)
                                                                                                               tCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGc                                                                                                        >  1:1785543/1‑78 (MQ=255)
                                                                                                               tCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAg                                >  1:1080474/1‑150 (MQ=255)
                                                                                                                  cTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCg                             <  1:5671290/150‑1 (MQ=255)
                                                                                                                   tggtggCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGc                            >  1:2742644/1‑150 (MQ=255)
                                                                                                                     gtggCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTg                          <  1:4423619/150‑1 (MQ=255)
                                                                                                                           gCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTg                    <  2:4830756/150‑1 (MQ=255)
                                                                                                                           gCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTg                    <  1:1653040/150‑1 (MQ=255)
                                                                                                                            cGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTgg                   <  1:3534566/150‑1 (MQ=255)
                                                                                                                            cGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTgg                   <  1:5172203/150‑1 (MQ=255)
                                                                                                                            cGGTCATCCACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTgg                   <  2:1511096/150‑1 (MQ=255)
                                                                                                                             ggTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTggg                  >  2:2306951/1‑150 (MQ=255)
                                                                                                                                cATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTGGGCTc               <  2:3613247/150‑1 (MQ=255)
                                                                                                                                 aTCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTGGGCTcg              >  2:627305/1‑150 (MQ=255)
                                                                                                                                 aTCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTGGGCTcg              >  1:6004022/1‑150 (MQ=255)
                                                                                                                                  tCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCggt                                                                                                   >  2:862053/1‑64 (MQ=255)
                                                                                                                                  tCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTGGGCTcgc             >  2:3922410/1‑150 (MQ=255)
                                                                                                                                     aCCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTGGGCTcgcagg          <  2:4794408/150‑1 (MQ=255)
                                                                                                                                     aCCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTGGGCTcgcagg          >  1:2940772/1‑150 (MQ=255)
                                                                                                                                       cGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTGGGCTcgcaggcg        >  2:684245/1‑150 (MQ=255)
                                                                                                                                        gCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTGGGCTcgcaggcgc       >  1:1211512/1‑150 (MQ=255)
                                                                                                                                         cTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTGGGCTcgcaggcgca      <  2:961069/150‑1 (MQ=255)
                                                                                                                                           cGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTGGGCTcgcaggcgcagg    <  2:6124134/150‑1 (MQ=255)
                                                                                                                                            ggTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTGGGCTCGCAGGCGCAGGt   <  1:2537079/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTGGGCTCGCAGGCGCAGGTA  >  NC_002947/1851095‑1851385

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