Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,851,157 A→C 6.3% A67A (GCA→GCC relA → ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,851,1570AC 6.3% 228.3 / 6.7 80A67A (GCA→GCCrelAATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base A (40/35);  new base C (2/3);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

TCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGA  >  NC_002947/1851009‑1851302
                                                                                                                                                    |                                                                                                                                                  
tcGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCcaggcag                                                                                                                                                   <  2:105520/150‑1 (MQ=255)
tcGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCcaggcag                                                                                                                                                   <  1:701768/150‑1 (MQ=255)
tcGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCcaggcag                                                                                                                                                   >  2:1254964/1‑150 (MQ=255)
  gAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCaggcaggc                                                                                                                                                 <  1:597040/150‑1 (MQ=255)
  gAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAACCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGCGCGGATGGTACGTCCAGCTTCCaggcaggc                                                                                                                                                 <  1:1279408/150‑1 (MQ=255)
   aaGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGcc                                                                                                                                                >  1:393189/1‑150 (MQ=255)
     gCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTg                                                                                                                                              >  2:542892/1‑150 (MQ=255)
       aTGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGa                                                                                                                                            >  1:1062381/1‑150 (MQ=255)
          gTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAAt                                                                                                                                         <  1:1921715/150‑1 (MQ=255)
           ttGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCCTGCCTGGaagct                                                                                                                                        <  2:1970688/150‑4 (MQ=255)
           ttGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCCTGCCTGGaagct                                                                                                                                        <  2:1338238/150‑4 (MQ=255)
           ttGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATc                                                                                                                                        <  2:2094099/150‑1 (MQ=255)
            tGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCCTGCCTGGaagctg                                                                                                                                       <  2:316055/150‑5 (MQ=255)
              gACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGcc                                                                                                                                     >  2:568423/1‑150 (MQ=255)
                  aTGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGaaa                                                                                                                                 >  2:1828689/1‑150 (MQ=255)
                  aTGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGaaa                                                                                                                                 <  2:358687/150‑1 (MQ=255)
                       gTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTg                                                                                                                            <  1:568423/150‑1 (MQ=255)
                         aaGCGTGGATTGGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGc                                                                                                                          <  2:965794/150‑1 (MQ=255)
                          aGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAACATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCt                                                                                                                         <  1:2127384/150‑1 (MQ=255)
                           gCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTg                                                                                                                        >  2:960396/1‑150 (MQ=255)
                             gTCTATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGAc                                                                                                                      <  2:870413/150‑1 (MQ=255)
                              tCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGAcc                                                                                                                     >  1:1470532/1‑150 (MQ=255)
                                  ttCGGCACTGGATCGCGCAGCGCTGAAAAAAGCCTGCGAGGTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCAAGGAAATCGCCGAAATTCTGGCTGACCTCaa                                                                                                                 >  1:1576388/1‑150 (MQ=255)
                                   tCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAg                                                                                                                >  1:372675/1‑150 (MQ=255)
                                     ggCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCt                                                                                                              >  2:463960/1‑150 (MQ=255)
                                      gCACTCGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGCGGTAGAGAAAAAGGGCAACCCCGCCACGTTTCGCTGGGCTGATGGTACGTTCATCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTc                                                                                                             <  1:542892/150‑1 (MQ=255)
                                         cTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGAc                                                                                                          <  1:876282/150‑1 (MQ=255)
                                           ggATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCa                                                                                                        <  1:1591739/150‑1 (MQ=255)
                                           ggATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCa                                                                                                        <  1:463960/150‑1 (MQ=255)
                                           ggATCGCGCAGCGCTGAAAGAAGACTGGGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCa                                                                                                        <  2:1387377/150‑1 (MQ=255)
                                            gATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAg                                                                                                       <  1:925458/150‑1 (MQ=255)
                                                       gcTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGCCTCCCtggt                                                                                            <  1:2301255/150‑1 (MQ=255)
                                                         tGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCtggtgg                                                                                          >  2:275804/1‑150 (MQ=255)
                                                         tGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCtggtgg                                                                                          <  1:2227101/150‑1 (MQ=255)
                                                          gAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGc                                                                                         <  1:275804/150‑1 (MQ=255)
                                                               aagCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCgg                                                                                    <  1:1828689/150‑1 (MQ=255)
                                                                 gCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCTCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTc                                                                                  <  2:1069666/150‑1 (MQ=255)
                                                                 gCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTc                                                                                  >  2:1619129/1‑150 (MQ=255)
                                                                    tGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGAcc                                                                                                                     <  2:2332735/112‑1 (MQ=255)
                                                                    tGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATc                                                                               >  2:566931/1‑150 (MQ=255)
                                                                    tGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAAGCATGCCTGGAAATCGCCGAAATTCTGGCTGAcc                                                                                                                     >  1:2332735/1‑112 (MQ=255)
                                                                               tctTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCa                                                                                                        <  2:1258483/114‑1 (MQ=255)
                                                                               tctTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCa                                                                                                        >  1:1258483/1‑114 (MQ=255)
                                                                               tctTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGt                                                                   >  1:2152662/1‑150 (MQ=255)
                                                                               tctTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGt                                                                   >  2:977867/1‑150 (MQ=255)
                                                                                 ttGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTgc                                                                 >  1:1746108/1‑150 (MQ=255)
                                                                                    aGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGcc                                                                                                                                     <  2:1239539/80‑1 (MQ=255)
                                                                                    aGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGcc                                                                                                                                     >  1:1239539/1‑80 (MQ=255)
                                                                                      gTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAg                                                            >  1:330079/1‑150 (MQ=255)
                                                                                      gTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAg                                                            <  1:1266556/150‑1 (MQ=255)
                                                                                      gTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCATGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAg                                                            >  2:116152/1‑150 (MQ=255)
                                                                                         gagaAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGAc                                                                                                          <  1:194284/102‑1 (MQ=255)
                                                                                         gagaAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGAc                                                                                                          >  2:194284/1‑102 (MQ=255)
                                                                                           gaAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCaa                                                       <  2:1746108/150‑1 (MQ=255)
                                                                                            aaaaaGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGa                                                             >  1:430476/1‑143 (MQ=255)
                                                                                            aaaaaGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGa                                                             <  2:430476/143‑1 (MQ=255)
                                                                                            aaaaaGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAg                                                      <  1:703106/150‑1 (MQ=255)
                                                                                            aaaaaGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAg                                                      >  2:1307746/1‑150 (MQ=255)
                                                                                            aaaaaGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAg                                                      >  2:647622/1‑150 (MQ=255)
                                                                                              aaaGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGt                                                    <  1:647622/150‑1 (MQ=255)
                                                                                                 gggCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCAc                                                 >  1:122090/1‑150 (MQ=255)
                                                                                                 gggCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCAc                                                 >  1:113446/1‑150 (MQ=255)
                                                                                                   gCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCAccc                                               >  1:957508/1‑150 (MQ=255)
                                                                                                   gCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCAccc                                               <  2:1062381/150‑1 (MQ=255)
                                                                                                      aCCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCg                                            >  1:2018762/1‑150 (MQ=255)
                                                                                                        ccGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCCAAACGATGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGcc                                          >  2:686526/1‑150 (MQ=255)
                                                                                                            ccAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAgg                                      <  1:977867/150‑1 (MQ=255)
                                                                                                              aaGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTc                                    >  2:1640450/1‑150 (MQ=255)
                                                                                                              aaGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTc                                    >  2:282838/1‑150 (MQ=255)
                                                                                                                  aTTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTcagcc                                >  1:1717155/1‑150 (MQ=255)
                                                                                                                  aTTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTcagcc                                >  1:1718288/1‑150 (MQ=255)
                                                                                                                  aTTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTcagcc                                >  2:1417479/1‑150 (MQ=255)
                                                                                                                           cGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTg                       <  2:1470532/150‑1 (MQ=255)
                                                                                                                            ggATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTAACGCGCGACCAGGCCTCCCTGGTGGCTGCGGTCAGCTACCGCTCGGTGGGCGAGGGCAAGGTCACCCTCGCAGAGGTCAGCCAGCGGTTTgg                      <  2:1348962/150‑1 (MQ=255)
                                                                                                                                  tACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCggt                >  1:1158485/1‑150 (MQ=255)
                                                                                                                                   aCGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCggtg               <  1:1417479/150‑1 (MQ=255)
                                                                                                                                   aCGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCggtg               <  2:2152662/150‑1 (MQ=255)
                                                                                                                                     gTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCggt                <  2:1547286/147‑1 (MQ=255)
                                                                                                                                     gTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCggt                >  1:1547286/1‑147 (MQ=255)
                                                                                                                                      tCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGgtggtg            >  1:1102229/1‑150 (MQ=255)
                                                                                                                                        cAGCTTACAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTc          <  1:1307746/150‑1 (MQ=255)
                                                                                                                                           tttCCAGGACTGCCTGGAAAACGACGAAACTCCCGCTAACCCCAAGCTCGACAAGGACTCCCAGGTGGCGGAGGTCATCTACCGCTCGGTGCGCTAGGTCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCaa       >  2:883028/2‑150 (MQ=255)
                                                                                                                                            ttCCAGGCCTGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAg      >  2:1814049/1‑150 (MQ=255)
                                                                                                                                            ttCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTGTGCTGGGGTGCGCGAGGGCAAGGTCACCCTCGCCGACGTCAGCACGCGGTTTGGGCCGGTGTTGTCCAAg      >  1:43461/1‑150 (MQ=255)
                                                                                                                                              ccAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGCCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCt    >  1:1110698/1‑150 (MQ=255)
                                                                                                                                                aggcaggcCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAg                                                            >  1:2217244/1‑92 (MQ=255)
                                                                                                                                                aggcaggcCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAg                                                            <  2:2217244/92‑1 (MQ=255)
                                                                                                                                                aggcaggcCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGACCGCAGGTGTCCAATCTGa  >  2:1945555/1‑150 (MQ=255)
                                                                                                                                                  gcaggcCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATc                                                                               <  1:1751008/72‑1 (MQ=255)
                                                                                                                                                  gcaggcCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATc                                                                               >  2:1751008/1‑72 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGA  >  NC_002947/1851009‑1851302

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