Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,563,372 C→A 5.6% R294L (CGT→CTT)  purA ← adenylosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,563,3720CA 5.6% 1032.4 / 32.8 304R294L (CGT→CTT) purAadenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base C (130/157);  new base A (7/10);  total (137/167)
Fisher's exact test for biased strand distribution p-value = 8.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 309 total aligned reads displayed.

TGTTGATGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTG  >  NC_002947/5563223‑5563521
                                                                                                                                                     |                                                                                                                                                     
 gTTGATGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCCCCACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACg                                                                                                                                                      >  2:1347725/1‑150 (MQ=255)
  ttGATGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGc                                                                                                                                                     >  1:3687843/1‑150 (MQ=255)
    gATGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCtt                                                                                                                                                   <  1:6180686/150‑1 (MQ=255)
    gATGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCtt                                                                                                                                                   >  2:4427695/1‑150 (MQ=255)
    gATGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCtt                                                                                                                                                   >  2:3135538/1‑150 (MQ=255)
    gATGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCtt                                                                                                                                                   <  2:6835496/150‑1 (MQ=255)
     aTGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTg                                                                                                                                                  >  1:6048722/1‑150 (MQ=255)
     aTGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTg                                                                                                                                                  <  1:2528970/150‑1 (MQ=255)
      tGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTgg                                                                                                                                                 <  1:5599872/150‑1 (MQ=255)
          ttCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCAAGCGTGGCCac                                                                                                                                             <  1:4907101/150‑2 (MQ=255)
          ttCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCAAGCGTGGCCac                                                                                                                                             <  1:4792159/150‑2 (MQ=255)
            ccAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAgg                                                                                                                                           <  1:1484922/150‑1 (MQ=255)
            ccAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCAAGCGTGGCCacga                                                                                                                                           <  2:20497/150‑4 (MQ=255)
            ccAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCAAGCGTGGCCacga                                                                                                                                           <  2:1313827/150‑4 (MQ=255)
             cAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGt                                                                                                                                          >  2:5526047/1‑150 (MQ=255)
             cAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGt                                                                                                                                          >  2:3893796/1‑150 (MQ=255)
               ggCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAg                                                                                                                                        >  2:6558852/1‑150 (MQ=255)
                 ccGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCa                                                                                                                                      >  2:2572553/1‑150 (MQ=255)
                   gTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCAcc                                                                                                                                    >  1:1208678/1‑150 (MQ=255)
                    tCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCg                                                                                                                                   >  2:5066927/1‑150 (MQ=255)
                    tCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACAGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCg                                                                                                                                   <  2:5021797/150‑1 (MQ=255)
                     ccAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCgg                                                                                                                                  >  2:1811619/1‑150 (MQ=255)
                       aGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGtc                                                                                                                                >  2:2786070/1‑150 (MQ=255)
                       aGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGtc                                                                                                                                <  1:103192/150‑1 (MQ=255)
                            gTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCtcgtc                                                                                                                           >  1:5697422/1‑150 (MQ=255)
                             tCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCtcgtcg                                                                                                                          >  1:3913459/1‑150 (MQ=255)
                                aGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAAc                                                                                                                       <  2:577320/150‑1 (MQ=255)
                                 gCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACa                                                                                                                      <  1:6690964/150‑1 (MQ=255)
                                   ttGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGt                                                                                                                    >  1:6778117/1‑150 (MQ=255)
                                   ttGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGt                                                                                                                    >  2:6080017/1‑150 (MQ=255)
                                    tGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGtt                                                                                                                   <  2:5387757/150‑1 (MQ=255)
                                       tCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTtcgg                                                                                                                >  2:689076/1‑150 (MQ=255)
                                         aGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTtcggtc                                                                                                              <  2:6751949/150‑1 (MQ=255)
                                         aGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTtcggtc                                                                                                              <  2:2948142/150‑1 (MQ=255)
                                          ggCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTtcggtcg                                                                                                             >  2:1941113/1‑150 (MQ=255)
                                               aTGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAAc                                                                                                        <  1:6364460/150‑1 (MQ=255)
                                                  cccGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGa                                                                                                     <  1:1778488/150‑1 (MQ=255)
                                                    cGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGAcc                                                                                                   >  1:4114434/1‑150 (MQ=255)
                                                    cGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGAcc                                                                                                   <  1:1657670/150‑1 (MQ=255)
                                                        aTGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGaa                                                                                               <  2:5577892/150‑1 (MQ=255)
                                                        aTGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGaa                                                                                               <  2:6373485/150‑1 (MQ=255)
                                                          gCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAAcc                                                                                             >  1:6873779/1‑150 (MQ=255)
                                                               tGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCAcgc                                                                                        <  1:779114/150‑1 (MQ=255)
                                                                 tcGTCGATGGCGCGACTCAGGATCACGGCCTCGACCCTGCCGCAACGACGGGCTCTGCCGGTGGTGGAACCGACCTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGATCAGTTCGGTCGGGAACGGACCGGAACCCACGCGa                                                                                      <  2:6626224/149‑1 (MQ=255)
                                                                  cGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAg                                                                                     <  2:6409824/150‑1 (MQ=255)
                                                                    tCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtg                                                                                   <  1:5966174/150‑1 (MQ=255)
                                                                    tCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtg                                                                                   >  1:5488574/1‑150 (MQ=255)
                                                                    tCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtg                                                                                   <  1:4109141/150‑1 (MQ=255)
                                                                    tCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtg                                                                                   >  1:3105900/1‑150 (MQ=255)
                                                                      gATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtggt                                                                                 <  1:515742/150‑1 (MQ=255)
                                                                      gATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtggt                                                                                 <  1:3298327/150‑1 (MQ=255)
                                                                       aTGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtggtg                                                                                <  1:5356454/150‑1 (MQ=255)
                                                                            gcgACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGc                                                                           >  1:125446/1‑150 (MQ=255)
                                                                            gcgACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGc                                                                           >  2:86820/1‑150 (MQ=255)
                                                                                cGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTg                                                                       <  2:4414510/150‑1 (MQ=255)
                                                                                 gCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTgg                                                                      >  1:907424/1‑150 (MQ=255)
                                                                                 gCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTgg                                                                      <  1:2019925/150‑1 (MQ=255)
                                                                                      aTCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATa                                                                 <  2:1487556/150‑1 (MQ=255)
                                                                                          cGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGa                                                             <  1:2845342/150‑1 (MQ=255)
                                                                                            gCATCGAACCAGCCGCCACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCTCCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAgg                                                           <  2:245396/150‑1 (MQ=255)
                                                                                            gCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCAAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAgg                                                           <  2:5619344/150‑1 (MQ=255)
                                                                                             cATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGa                                                          <  1:1085986/150‑1 (MQ=255)
                                                                                              aTCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGAt                                                         <  1:5538140/150‑1 (MQ=255)
                                                                                              aTCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGAt                                                         <  2:5968113/150‑1 (MQ=255)
                                                                                              aTCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGAt                                                         <  2:367334/150‑1 (MQ=255)
                                                                                               tCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATg                                                        <  1:3355109/150‑1 (MQ=255)
                                                                                                cGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGt                                                       >  1:2431416/1‑150 (MQ=255)
                                                                                                 gAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTa                                                      >  2:4401166/1‑150 (MQ=255)
                                                                                                     cAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCg                                                  >  2:2354375/1‑150 (MQ=255)
                                                                                                      aGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGa                                                 <  1:3203905/150‑1 (MQ=255)
                                                                                                      aGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGa                                                 <  2:2850462/150‑1 (MQ=255)
                                                                                                       gCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGa                                                             <  1:2725332/137‑1 (MQ=255)
                                                                                                        ccGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAgg                                               <  2:3622346/150‑1 (MQ=255)
                                                                                                        ccGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAgg                                               <  2:6694068/150‑1 (MQ=255)
                                                                                                                 cgGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCggg                                      <  1:2547381/150‑1 (MQ=255)
                                                                                                                  gggCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGc                                     <  1:6394209/150‑1 (MQ=255)
                                                                                                                  gggCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGGAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGc                                     <  1:2640945/150‑1 (MQ=255)
                                                                                                                   ggCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGcc                                    <  1:4463182/150‑1 (MQ=255)
                                                                                                                   ggCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGcc                                    >  1:4330389/1‑150 (MQ=255)
                                                                                                                    gCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCg                                   >  1:5926448/1‑150 (MQ=255)
                                                                                                                       cGGCCGGTGGTGGAACCGAACTCGTGGCCACGCGTGGCCAAGGTAGAACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGCCCGGAACCCACGCGAGGGGTGTAGGCCTTGGTGAGACCGAGGATGTAGTCGGGGTACATCGGGCCGACg                                >  1:3038679/1‑150 (MQ=255)
                                                                                                                             gtggtgGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGaa                          >  2:1109772/1‑150 (MQ=255)
                                                                                                                                  ggAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGt                     >  2:3069572/1‑150 (MQ=255)
                                                                                                                                  ggAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGt                     <  1:5615862/150‑1 (MQ=255)
                                                                                                                                   gAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTg                    <  2:4176029/150‑1 (MQ=255)
                                                                                                                                    aaCCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTgg                   <  2:4694515/150‑1 (MQ=255)
                                                                                                                                    aaCCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTgg                   <  1:5844016/150‑1 (MQ=255)
                                                                                                                                       cGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGaaa                <  1:869418/150‑1 (MQ=255)
                                                                                                                                         aaCTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATa              >  2:2726081/1‑150 (MQ=255)
                                                                                                                                            tCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCg                                                              >  2:1217483/1‑99 (MQ=255)
                                                                                                                                            tCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATacca           >  1:292007/1‑150 (MQ=255)
                                                                                                                                            ccttGGCCAAGCGTGGCCACGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATacca           >  1:6539473/4‑150 (MQ=255)
                                                                                                                                            ccttGGCCAAGCGTGGCCACGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATacca           >  2:3320765/4‑150 (MQ=255)
                                                                                                                                             cttGGCCAAGCGTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATaccac          >  2:4123784/3‑150 (MQ=255)
                                                                                                                                              gTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATaccacc         >  2:4928275/1‑150 (MQ=255)
                                                                                                                                              gTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATaccacc         >  2:4612968/1‑150 (MQ=255)
                                                                                                                                               tGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCg        >  1:2153685/1‑150 (MQ=255)
                                                                                                                                                  ccACGCTGGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGTATTTCCACCggtt     <  2:5224111/150‑3 (MQ=255)
                                                                                                                                                    aCGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGTATTTCCACCggttcc   <  1:1919159/150‑5 (MQ=255)
                                                                                                                                                     cGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAgtg  <  1:1229579/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTTGATGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTG  >  NC_002947/5563223‑5563521

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: