Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 772,181 T→C 6.0% intergenic (+39/‑160) PP_0663 → / → PP_0664 AsnC family transcriptional regulator/homoserine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947772,1810TC 6.0% 412.6 / 10.5 132intergenic (+39/‑160)PP_0663/PP_0664AsnC family transcriptional regulator/homoserine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (62/62);  new base C (5/3);  total (67/65)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

Only 100 of 132 total aligned reads displayed.

GAGGCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTGGACAG  >  NC_002947/772038‑772329
                                                                                                                                               |                                                                                                                                                    
 aGGCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGcc                                                                                                                                               >  1:2080192/1‑150 (MQ=255)
 aGGCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGcc                                                                                                                                               >  1:1634463/1‑150 (MQ=255)
 aGGCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGcc                                                                                                                                               >  1:4903096/1‑150 (MQ=255)
  ggCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGccc                                                                                                                                              <  2:1736444/150‑1 (MQ=255)
   gCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCa                                                                                                                                             >  1:4081301/1‑150 (MQ=255)
   gCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCa                                                                                                                                             <  1:1574687/150‑1 (MQ=255)
      tCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAggc                                                                                                                                          >  2:4452266/1‑150 (MQ=255)
      tCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAggc                                                                                                                                          >  1:3753018/1‑150 (MQ=255)
      tCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAggc                                                                                                                                          >  1:513649/1‑150 (MQ=255)
      tCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAggc                                                                                                                                          >  1:4389356/1‑150 (MQ=255)
        cAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAggcgg                                                                                                                                        <  2:1189966/150‑1 (MQ=255)
         aGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAggcggc                                                                                                                                       >  2:1847685/1‑150 (MQ=255)
          gCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGgcggcg                                                                                                                                      >  2:4647601/1‑150 (MQ=255)
           cccGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGt                                                                                                                                     >  2:3301904/1‑150 (MQ=255)
            ccGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGtt                                                                                                                                    >  1:4331879/1‑150 (MQ=255)
            ccGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGtt                                                                                                                                    <  1:3616201/150‑1 (MQ=255)
            ccGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGCGCGAACCCAGGCGGCGtt                                                                                                                                    <  2:1295607/150‑1 (MQ=255)
             cGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGCGCGAACCCAGGCGGCGTtg                                                                                                                                   <  2:515190/150‑1 (MQ=255)
             cGGAGCGGATGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTtg                                                                                                                                   >  1:4938106/1‑150 (MQ=255)
              ggAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTtgt                                                                                                                                  >  1:1337551/1‑150 (MQ=255)
                aGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGc                                                                                                                                <  2:2080192/150‑1 (MQ=255)
                  cGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGt                                                                                                                              >  1:1868948/1‑150 (MQ=255)
                  cGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGCGCGAACCCAGGCGGCGTTGTGcta                                                                                                                              <  1:3559054/150‑3 (MQ=255)
                     gTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGtttt                                                                                                                           >  2:878431/1‑150 (MQ=255)
                       gCACCAGGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGc                                                                                                                         <  1:3247930/150‑1 (MQ=255)
                                  ggaaggaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAgg                                                                                                              >  2:4158276/1‑150 (MQ=255)
                                    aaggaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTg                                                                                                            >  1:377154/1‑150 (MQ=255)
                                      ggaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTgcg                                                                                                          >  2:386635/1‑150 (MQ=255)
                                       gaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTgcgc                                                                                                         >  2:876616/1‑150 (MQ=255)
                                          aTCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCt                                                                                                      >  1:4030264/1‑150 (MQ=255)
                                            cGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCttt                                                                                                    <  1:2559733/150‑1 (MQ=255)
                                            cGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCttt                                                                                                    <  2:4792391/150‑1 (MQ=255)
                                             gCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTg                                                                                                   >  2:1406607/1‑150 (MQ=255)
                                               tGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTTGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCa                                                                                                 >  1:3079537/1‑150 (MQ=255)
                                                gCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCAc                                                                                                <  2:5002672/150‑1 (MQ=255)
                                                 ccATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACt                                                                                               >  1:3705251/1‑150 (MQ=255)
                                                  cATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTg                                                                                              <  1:4844846/150‑1 (MQ=255)
                                                      cccGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCa                                                                                          >  1:1954730/1‑150 (MQ=255)
                                                       ccGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAg                                                                                         >  2:5034508/1‑150 (MQ=255)
                                                          ggggggAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGAtt                                                                                      <  1:3618059/150‑1 (MQ=255)
                                                          gggTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGAtt                                                                                      >  1:3134487/1‑150 (MQ=255)
                                                             tGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGc                                                                                   <  1:4109001/150‑1 (MQ=255)
                                                               gagaACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCa                                                                                 <  1:3437811/150‑1 (MQ=255)
                                                                agaACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCaa                                                                                >  1:3968511/1‑150 (MQ=255)
                                                                 gaACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAATCTGGCCagt                                                                               <  2:1609333/150‑3 (MQ=255)
                                                                    cacGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCaataaa                                                                            <  2:1191804/150‑1 (MQ=255)
                                                                    cacGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCaataaa                                                                            >  1:17901/1‑150 (MQ=255)
                                                                         tGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAg                                                                       <  2:2967450/150‑1 (MQ=255)
                                                                         tGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAg                                                                       <  2:513649/150‑1 (MQ=255)
                                                                         tGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAg                                                                       <  2:3296323/150‑1 (MQ=255)
                                                                            cTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCgg                                                                    >  2:2435839/1‑150 (MQ=255)
                                                                             ttcgttcgTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGa                                                                   >  2:232664/1‑150 (MQ=255)
                                                                             ttcgttcgTCCTGGCCGCGGGGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGa                                                                   <  2:4058119/150‑1 (MQ=255)
                                                                             ttcgttcgTCCTGGCCGAGGTGGTCAGAACCGCGTCTTTGCAATGGCTAGCACAACGCCGCCTCGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTCGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGa                                                                   >  2:5012928/1‑150 (MQ=255)
                                                                                gttcgtCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTg                                                                >  1:1505711/1‑150 (MQ=255)
                                                                                gttcgtCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTg                                                                >  1:539902/1‑150 (MQ=255)
                                                                                gttcgtCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTg                                                                >  1:1011494/1‑150 (MQ=255)
                                                                                  tcgtCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGcc                                                              <  1:123863/150‑1 (MQ=255)
                                                                                  tcgtCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGcc                                                              <  1:123903/150‑1 (MQ=255)
                                                                                    gtCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCat                                                            <  1:1631340/150‑1 (MQ=255)
                                                                                     tCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCata                                                           >  2:4390510/1‑150 (MQ=255)
                                                                                       cTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATAtt                                                         >  1:771093/1‑150 (MQ=255)
                                                                                        tGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTc                                                        <  2:1132131/150‑1 (MQ=255)
                                                                                          gCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCaa                                                      <  2:3968511/150‑1 (MQ=255)
                                                                                             gcgGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTg                                                   <  2:3273209/150‑1 (MQ=255)
                                                                                                gtggtCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCa                                                <  1:1592059/150‑1 (MQ=255)
                                                                                                gtggtCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCa                                                <  1:792393/150‑1 (MQ=255)
                                                                                                  ggtCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGATTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGt                                              <  1:314760/150‑1 (MQ=255)
                                                                                                     cTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGCGCa                                           <  2:816204/150‑1 (MQ=255)
                                                                                                      tGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCAc                                          <  1:4141455/150‑1 (MQ=255)
                                                                                                       gAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACg                                         >  2:4593079/1‑150 (MQ=255)
                                                                                                        aaCCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACgg                                        >  1:3927255/1‑150 (MQ=255)
                                                                                                           cgcgTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGctc                                     <  1:3305212/150‑1 (MQ=255)
                                                                                                             cgTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCtt                                   <  1:1984754/150‑1 (MQ=255)
                                                                                                                 tttGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGc                               <  2:2973194/150‑1 (MQ=255)
                                                                                                                 tttGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGATTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGc                               <  2:2973188/150‑1 (MQ=255)
                                                                                                                   tGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTa                             <  1:3287620/150‑1 (MQ=255)
                                                                                                                   tGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTa                             <  1:3287656/150‑1 (MQ=255)
                                                                                                                      aTTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTc                          >  2:2086242/1‑150 (MQ=255)
                                                                                                                      aTTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTc                          <  1:4858311/150‑1 (MQ=255)
                                                                                                                       ttGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTcc                         <  2:3753018/150‑1 (MQ=255)
                                                                                                                          gCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCa                      <  1:49662/150‑1 (MQ=255)
                                                                                                                          gCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCa                      <  1:49683/150‑1 (MQ=255)
                                                                                                                            tAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATc                    >  1:4322814/1‑150 (MQ=255)
                                                                                                                            tAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATc                    >  1:2975279/1‑150 (MQ=255)
                                                                                                                               cacaACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGAc                 >  1:4547732/1‑150 (MQ=255)
                                                                                                                                acaACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGAcc                <  2:4389356/150‑1 (MQ=255)
                                                                                                                                   aCGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCAc             >  2:2112968/1‑150 (MQ=255)
                                                                                                                                    cgccgcCTGGGCGCGAACCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACt            >  1:4838771/1‑150 (MQ=255)
                                                                                                                                     gccgccTGGGCGCGAACCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTg           >  1:3391716/1‑150 (MQ=255)
                                                                                                                                      ccgccTGGGCGCGAACCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGc          >  1:3018350/1‑150 (MQ=255)
                                                                                                                                       cgccTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGcc         >  2:435025/1‑150 (MQ=255)
                                                                                                                                        gccTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCt        <  2:4331879/150‑1 (MQ=255)
                                                                                                                                        gccTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCt        <  1:1847685/150‑1 (MQ=255)
                                                                                                                                         ccTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGCCCTCATCGACCCACTGCCTg       <  2:1355718/150‑1 (MQ=255)
                                                                                                                                          cTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTgg      >  2:115765/1‑150 (MQ=255)
                                                                                                                                          cTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTgg      >  1:3433011/1‑150 (MQ=255)
                                                                                                                                           tGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTgga     >  2:2591401/1‑150 (MQ=255)
                                                                                                                                             ggTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTggaca   <  1:1924052/150‑1 (MQ=255)
                                                                                                                                              gTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTggacag  <  2:4903096/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GAGGCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTGGACAG  >  NC_002947/772038‑772329

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: