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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,560,3450CA 5.4% 429.1 / 10.4 129R294L (CGT→CTT) purAadenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base C (62/60);  new base A (3/4);  total (65/64)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 129 total aligned reads displayed.

TTGATGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAG  >  NC_002947/5560198‑5560492
                                                                                                                                                   |                                                                                                                                                   
ttGATGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGc                                                                                                                                                   <  1:4914319/150‑1 (MQ=255)
  gATGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCtt                                                                                                                                                 <  2:4847222/150‑1 (MQ=255)
  gATGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCtt                                                                                                                                                 >  2:4247127/1‑150 (MQ=255)
   aTGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTg                                                                                                                                                >  1:4291303/1‑150 (MQ=255)
   aTGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTg                                                                                                                                                >  2:2572404/1‑150 (MQ=255)
    tGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTgg                                                                                                                                               >  2:4279995/1‑150 (MQ=255)
    tGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTgg                                                                                                                                               >  1:2136478/1‑150 (MQ=255)
    tGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTgg                                                                                                                                               >  1:2136415/1‑150 (MQ=255)
     ggTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGc                                                                                                                                              <  1:3240102/150‑1 (MQ=255)
     ggTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGc                                                                                                                                              <  1:3240069/150‑1 (MQ=255)
     ggTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGc                                                                                                                                              <  1:2261774/150‑1 (MQ=255)
      gTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGcc                                                                                                                                             >  1:1034270/1‑150 (MQ=255)
       tttCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCa                                                                                                                                            >  1:915908/1‑150 (MQ=255)
         tCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAg                                                                                                                                          <  1:1033969/150‑1 (MQ=255)
         tCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAg                                                                                                                                          <  2:1789647/150‑1 (MQ=255)
          ccAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCTTGGCCAAGCGTGGCCacga                                                                                                                                         <  2:1209402/150‑4 (MQ=255)
          ccAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCTTGGCCAAGCGTGGCCacga                                                                                                                                         <  2:1778322/150‑4 (MQ=255)
          ccAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGCCCACGCTTGGCCAAgg                                                                                                                                         >  2:5345052/1‑150 (MQ=255)
           cAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGt                                                                                                                                        >  1:2088770/1‑150 (MQ=255)
           cAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGt                                                                                                                                        <  1:908860/150‑1 (MQ=255)
            aGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATTGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTa                                                                                                                                       >  1:4039079/1‑150 (MQ=255)
            aGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTa                                                                                                                                       >  1:280720/1‑150 (MQ=255)
              gCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCAAGCGTGGCCacgagttc                                                                                                                                     <  2:3129710/150‑8 (MQ=255)
                    cAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGt                                                                                                                               <  1:2151436/150‑1 (MQ=255)
                         cGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCtcgt                                                                                                                          <  1:5088335/150‑1 (MQ=255)
                             cAGCTTGGTCAGGCAGAGGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGaa                                                                                                                      <  2:490985/150‑1 (MQ=255)
                              aGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCCTCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAAc                                                                                                                     >  2:3558949/1‑150 (MQ=255)
                                 ttGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGt                                                                                                                  >  1:1865738/1‑150 (MQ=255)
                                   ggTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTtc                                                                                                                <  1:2847575/150‑1 (MQ=255)
                                   ggTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTtc                                                                                                                <  1:2724180/150‑1 (MQ=255)
                                   ggTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTtc                                                                                                                <  2:3713793/150‑1 (MQ=255)
                                    gTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTtcg                                                                                                               <  1:3588907/150‑1 (MQ=255)
                                     tCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTtcgg                                                                                                              >  1:3410429/1‑150 (MQ=255)
                                       aGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTtcggtc                                                                                                            >  1:5218298/1‑150 (MQ=255)
                                         gCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCTTGGCCAAGCGTGGCCacga                                                                                                                                         <  2:1935035/119‑4 (MQ=255)
                                         gCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTtcggtcgg                                                                                                          >  1:1292331/1‑150 (MQ=255)
                                            gATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGaa                                                                                                       <  1:3844441/150‑1 (MQ=255)
                                              tGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACg                                                                                                     <  2:3566067/150‑1 (MQ=255)
                                                 ccGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGAc                                                                                                  >  2:928224/1‑150 (MQ=255)
                                                   gagaTGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCg                                                                                                <  2:208943/150‑1 (MQ=255)
                                                      aTGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGaa                                                                                             <  2:2780523/150‑1 (MQ=255)
                                                          tGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCa                                                                                         >  1:4142741/1‑150 (MQ=255)
                                                            ttGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCAcg                                                                                       >  1:4457002/1‑150 (MQ=255)
                                                               aCGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGa                                                                                    <  2:2713258/150‑1 (MQ=255)
                                                                 gTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgt                                                                                  >  1:2258188/1‑150 (MQ=255)
                                                                 gTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgt                                                                                  >  1:807170/1‑150 (MQ=255)
                                                                    gATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtggt                                                                               >  1:4896174/1‑150 (MQ=255)
                                                                    gATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtggt                                                                               <  2:3841997/150‑1 (MQ=255)
                                                                      tGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGgtgt                                                                             >  1:4697603/1‑150 (MQ=255)
                                                                      tGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGgtgt                                                                             <  2:915908/150‑1 (MQ=255)
                                                                        gcgcGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAg                                                                           >  2:1111575/1‑150 (MQ=255)
                                                                        gcgcGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAg                                                                           >  2:586634/1‑150 (MQ=255)
                                                                          gcgACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGc                                                                         >  2:2969869/1‑150 (MQ=255)
                                                                              cGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTg                                                                     <  1:479574/150‑1 (MQ=255)
                                                                               gCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTgg                                                                    <  1:3403625/150‑1 (MQ=255)
                                                                                 aGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTg                                                                  >  1:2859774/1‑150 (MQ=255)
                                                                                 aGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTg                                                                  >  1:134083/1‑150 (MQ=255)
                                                                                     tCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATAc                                                              >  2:3690644/1‑150 (MQ=255)
                                                                                     tCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATAc                                                              >  1:4694887/1‑150 (MQ=255)
                                                                                     tCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATAc                                                              >  2:5048029/1‑150 (MQ=255)
                                                                                     tCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATAc                                                              >  2:5010436/1‑150 (MQ=255)
                                                                                       aCGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCg                                                            >  2:4293635/1‑150 (MQ=255)
                                                                                        cGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCGCGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGa                                                           >  1:2860806/1‑150 (MQ=255)
                                                                                         ggCATCGAACCAGCCGCATCGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGATCAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAg                                                          <  2:552430/150‑1 (MQ=255)
                                                                                         ggCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGa                                                                                              >  2:549305/1‑114 (MQ=255)
                                                                                         ggCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGa                                                                                              <  1:549305/114‑1 (MQ=255)
                                                                                          gCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGAACTGAACCCACGCGAGTGGTGTAGGCCTTGTTGATACCGAgg                                                         <  1:3109893/150‑1 (MQ=255)
                                                                                               gAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTa                                                    >  2:4757579/1‑150 (MQ=255)
                                                                                                  ccAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCCCGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTc                                                 >  1:1808656/1‑150 (MQ=255)
                                                                                                    aGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCGCGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGa                                               <  2:2860806/150‑1 (MQ=255)
                                                                                                    aGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGa                                               <  2:744944/150‑1 (MQ=255)
                                                                                                    aGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGa                                               <  1:4283844/150‑1 (MQ=255)
                                                                                                     gCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAg                                              <  1:4115869/150‑1 (MQ=255)
                                                                                                        gCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCCCGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTa                                           >  2:13277/1‑150 (MQ=255)
                                                                                                             gacgGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCg                                      <  1:4311014/150‑1 (MQ=255)
                                                                                                              acgGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCgg                                     >  1:209248/1‑150 (MQ=255)
                                                                                                              acgGGCACGGCCGGCGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCgg                                     >  2:4980687/1‑150 (MQ=255)
                                                                                                                gggCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCggg                                                                                                         <  2:2102373/80‑1 (MQ=255)
                                                                                                                gggCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGc                                   >  1:2812788/1‑150 (MQ=255)
                                                                                                                gggCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGc                                   >  1:683175/1‑150 (MQ=255)
                                                                                                                 ggCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGcc                                  <  2:5218298/150‑1 (MQ=255)
                                                                                                                  gCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtggtg                                                                              >  1:2087106/1‑105 (MQ=255)
                                                                                                                  gCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtggtg                                                                              <  2:2087106/105‑1 (MQ=255)
                                                                                                                       gCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGcc                            <  2:1034270/150‑1 (MQ=255)
                                                                                                                              gtggAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCg                     <  1:2220890/150‑1 (MQ=255)
                                                                                                                                ggAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGt                   >  1:5213976/1‑150 (MQ=255)
                                                                                                                                  aaCCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTgg                 <  1:2572404/150‑1 (MQ=255)
                                                                                                                                  aaCCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTgg                 <  2:1768722/150‑1 (MQ=255)
                                                                                                                                   aCCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGa                >  2:5282891/1‑150 (MQ=255)
                                                                                                                                      gAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAAt             <  1:2875442/150‑1 (MQ=255)
                                                                                                                                         accttGGCCAAGCGTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGACCCGGTGGAAATacc          >  1:1890992/5‑150 (MQ=255)
                                                                                                                                          tCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATacca         <  1:634682/150‑1 (MQ=255)
                                                                                                                                          tCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATacca         >  2:2886102/1‑150 (MQ=255)
                                                                                                                                           cttGGCCAAGCGTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTCCATCGGGCCGACGCCGGAACCGGTGGAAATaccac        >  2:3292621/3‑150 (MQ=255)
                                                                                                                                              ggCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCgg     >  1:2083479/1‑150 (MQ=255)
                                                                                                                                               gCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGc    <  1:4279995/150‑1 (MQ=255)
                                                                                                                                                 cACGCTTTGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAg  >  2:2293582/1‑150 (MQ=255)
                                                                                                                                                 cACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGa                                               >  2:5243871/1‑105 (MQ=255)
                                                                                                                                                 cACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGa                                               <  1:5243871/105‑1 (MQ=255)
                                                                                                                                                 cACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGt                                            <  2:2155397/108‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTGATGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAG  >  NC_002947/5560198‑5560492

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: