Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,558,426 T→G 5.1% K293T (AAG→ACG)  purA ← adenylosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,558,4260TG 5.1% 225.1 / 3.7 80K293T (AAG→ACG) purAadenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base T (34/41);  new base G (2/2);  total (37/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

ATGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTG  >  NC_002947/5558279‑5558575
                                                                                                                                                   |                                                                                                                                                     
aTGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTg                                                                                                                                                     <  2:1406309/150‑1 (MQ=255)
aTGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTg                                                                                                                                                     <  2:1956391/150‑1 (MQ=255)
 tGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTgg                                                                                                                                                    >  1:731400/1‑150 (MQ=255)
   gTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGCTGCTGTTGACGTCGATGGCGCGACGCAGGATGACGGCATCGCACCAGCCGCAACGACGGGCACGGCCGGGGGTGGAACCGAACTCGTGGCCACGCTTGGcc                                                                                                                                                  >  1:197985/1‑150 (MQ=255)
     ttCCAGGCCGTCCAGTACGTCGAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCaa                                                                                                                                                <  2:141093/150‑1 (MQ=255)
       ccAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGCACTCTTGGCCAAGCGTGGCCacga                                                                                                                                              <  1:1471924/150‑4 (MQ=255)
         aGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTa                                                                                                                                            >  2:712554/1‑150 (MQ=255)
         aGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTa                                                                                                                                            <  1:805058/150‑1 (MQ=255)
           gCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGc                                                                                                                                                   <  2:197150/141‑1 (MQ=255)
           gCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGc                                                                                                                                                   >  1:197150/1‑141 (MQ=255)
           gCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGc                                                                                                                                          >  2:1062663/1‑150 (MQ=255)
             cGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCAc                                                                                                                                        >  2:1917260/1‑150 (MQ=255)
                 cAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCTTGGCCAAGCGTGGCCac                                                                                                                                                <  1:2209618/138‑2 (MQ=255)
                 cAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCTTGGCCAAGCGTGGCCac                                                                                                                                                >  2:2209618/1‑137 (MQ=255)
                 cAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGt                                                                                                                                    >  1:2094241/1‑150 (MQ=255)
                  aGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGtc                                                                                                                                   >  1:793633/1‑150 (MQ=255)
                  aGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGtc                                                                                                                                   <  1:924865/150‑1 (MQ=255)
                   gTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGtct                                                                                                                                  >  1:2105230/1‑150 (MQ=255)
                    tACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGtctc                                                                                                                                 >  1:1882713/1‑150 (MQ=255)
                       gTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCtcgtc                                                                                                                              <  1:1062663/150‑1 (MQ=255)
                        tCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCtcgtcg                                                                                                                             >  1:601785/1‑150 (MQ=255)
                          cAGTTTGTTCATGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGaa                                                                                                                           <  2:1882713/150‑1 (MQ=255)
                           aGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAAc                                                                                                                          >  1:333242/1‑150 (MQ=255)
                            gCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCtt                                                                                                                                                      <  2:1869158/121‑1 (MQ=255)
                            gCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCtt                                                                                                                                                      >  1:1869158/1‑121 (MQ=255)
                                 gTCAGGCAGATGCCGGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTtcg                                                                                                                    <  2:1070645/150‑1 (MQ=255)
                                  tCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTtcgg                                                                                                                   >  2:370939/1‑150 (MQ=255)
                                  tCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCAGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTtcgg                                                                                                                   >  1:747951/1‑150 (MQ=255)
                                   cAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTtcggt                                                                                                                  <  2:2105230/150‑1 (MQ=255)
                                          aTGACCGAGATCCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGATGGTCGAACCGAACTAGTGGACCCGCTCGGACCAGGTAGCCCCGGTCTCGTCGAACAGTTCGGTCGGGAAc                                                                                                           <  1:1917260/150‑1 (MQ=255)
                                             cccGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGa                                                                                                        >  2:2051076/1‑150 (MQ=255)
                                                gagaTGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCg                                                                                                     <  1:798495/150‑1 (MQ=255)
                                                gagaTGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCg                                                                                                     >  1:979242/1‑150 (MQ=255)
                                                       tGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAg                                                                                                                                               <  2:1642974/101‑1 (MQ=255)
                                                       tGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAg                                                                                                                                               >  1:1642974/1‑101 (MQ=255)
                                                       tGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCa                                                                                              <  2:1111039/150‑1 (MQ=255)
                                                       tGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGATTCGGCCAGGCTCGCCAACACCGCACGGTTCTAACAGTTCGGTCGGGAACGGACCGGAACCCa                                                                                              >  1:1686818/1‑150 (MQ=255)
                                                            aCGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCaa                                                                                                                                                >  2:303665/1‑95 (MQ=255)
                                                            aCGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCaa                                                                                                                                                <  1:303665/95‑1 (MQ=255)
                                                            aCGTCGATGGCGCGACGCAGGATCACGGAATCGCCCCCGCGGGCCTGGAGCTGACGGGCGGTGGTGGAACCGCACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGa                                                                                         <  1:2271543/150‑1 (MQ=255)
                                                              gTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgt                                                                                       <  2:2022122/150‑1 (MQ=255)
                                                              gTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgt                                                                                       >  2:1316122/1‑150 (MQ=255)
                                                               tagATGGCGCGACGCAGTTTCACTGTATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtg                                                                                      <  2:1298441/148‑1 (MQ=255)
                                                               tCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtg                                                                                      <  2:333242/150‑1 (MQ=255)
                                                               tCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtg                                                                                      >  2:1559303/1‑150 (MQ=255)
                                                                 gATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtggt                                                                                    >  2:1185569/1‑150 (MQ=255)
                                                                  aTGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAAATACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAgtggtg                                                                                   <  1:712554/150‑1 (MQ=255)
                                                                       gcgACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGc                                                                              <  1:1786796/150‑1 (MQ=255)
                                                                       gcgACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGc                                                                              <  1:770764/150‑1 (MQ=255)
                                                                       gcgACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGc                                                                              >  2:770764/1‑150 (MQ=255)
                                                                           cGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTg                                                                          <  2:2094241/150‑1 (MQ=255)
                                                                            gCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTgg                                                                         >  1:1149493/1‑150 (MQ=255)
                                                                            gCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTgg                                                                         >  1:139795/1‑150 (MQ=255)
                                                                            gCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTgg                                                                         <  1:1585686/150‑1 (MQ=255)
                                                                            gCAGAATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTgg                                                                         <  2:139795/150‑1 (MQ=255)
                                                                                gATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGAt                                                                     >  2:1832216/1‑150 (MQ=255)
                                                                                  tCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACATTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGTCCTTGGCGATAc                                                                   >  2:911578/1‑150 (MQ=255)
                                                                                  tCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGa                                                                                         <  2:655895/128‑1 (MQ=255)
                                                                                  tCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGa                                                                                         >  1:655895/1‑128 (MQ=255)
                                                                                    aCGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGt                                                                                                                                             <  2:644330/74‑1 (MQ=255)
                                                                                    aCGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGt                                                                                                                                             >  1:644330/1‑74 (MQ=255)
                                                                                     cGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGa                                                                <  2:979242/150‑1 (MQ=255)
                                                                                     cGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGa                                                                <  2:731400/150‑1 (MQ=255)
                                                                                       gCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAgg                                                              <  2:793633/150‑1 (MQ=255)
                                                                                        cATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGa                                                             <  2:1588244/150‑1 (MQ=255)
                                                                                         aTCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGAt                                                            <  2:747951/150‑1 (MQ=255)
                                                                                          tCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATg                                                           <  1:370939/150‑1 (MQ=255)
                                                                                             aaCCACCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAATGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAg                                                        <  2:1225565/150‑1 (MQ=255)
                                                                                                 aGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGa                                                    >  1:1679886/1‑150 (MQ=255)
                                                                                                    cGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGt                                                 <  1:896533/150‑1 (MQ=255)
                                                                                                      cAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTAc                                               <  1:837533/150‑1 (MQ=255)
                                                                                                                         tggtggAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAAc                            <  2:2158363/150‑1 (MQ=255)
                                                                                                                               aaCCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTgg                      <  2:601785/150‑1 (MQ=255)
                                                                                                                                aCCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAgg                                                  <  1:565498/121‑1 (MQ=255)
                                                                                                                                aCCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAgg                                                  >  2:565498/1‑121 (MQ=255)
                                                                                                                                aCCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGa                     >  2:1795165/1‑150 (MQ=255)
                                                                                                                                 ccGAACTCGTGGCCACGCTTGGCCAGGGTAGCACCGGTCTCGTCGAACAGTtc                                                                                                                     <  2:1663068/53‑1 (MQ=255)
                                                                                                                                 ccGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTtc                                                                                                                     >  1:1663068/1‑53 (MQ=255)
                                                                                                                                   gAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCg                                                     <  2:2179937/115‑1 (MQ=255)
                                                                                                                                   gAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCg                                                     >  1:2179937/1‑115 (MQ=255)
                                                                                                                                      accttGGCCAAGCGTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATacc               >  2:2253968/5‑150 (MQ=255)
                                                                                                                                              cACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAg       >  1:1305340/1‑150 (MQ=255)
                                                                                                                                                   ttGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAgtggtg  >  1:1955953/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATGGTTTCCAGGCCGTCCAGTACGTCCAGCTTGGTCAGGCAGATGCCCGAGATGCTGTTGACGTCGATGGCGCGACGCAGGATCACGGCATCGAACCAGCCGCAACGACGGGCACGGCCGGTGGTGGAACCGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTG  >  NC_002947/5558279‑5558575

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