Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 772,181 T→C 8.8% 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 8.8% 430.5 / 22.1 147intergenic (+39/‑160)PP_0663/PP_0664AsnC family transcriptional regulator/homoserine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (74/60);  new base C (8/5);  total (82/65)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

Only 100 of 148 total aligned reads displayed.

GCGTCGAGGCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTGGACAG  >  NC_002947/772033‑772329
                                                                                                                                                    |                                                                                                                                                    
gcgTCGAGGCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCAGGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGtt                                                                                                                                                     >  1:2131102/1‑150 (MQ=255)
  gTCGAGGCCTCCAGCCCGGAGCGGGTGCACCATGTGGGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTGCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTcg                                                                                                                                                   >  2:4166877/1‑150 (MQ=255)
     gAGGCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCgcgc                                                                                                                                                <  1:2683671/150‑1 (MQ=255)
     gAGGCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCgcgc                                                                                                                                                <  1:622539/150‑1 (MQ=255)
       ggCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGccc                                                                                                                                              >  1:4612932/1‑150 (MQ=255)
        gCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCa                                                                                                                                             <  2:3769193/150‑1 (MQ=255)
         ccTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAg                                                                                                                                            <  1:4930819/150‑1 (MQ=255)
         ccTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAg                                                                                                                                            >  2:4376006/1‑150 (MQ=255)
          cTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAgg                                                                                                                                           >  1:1483106/1‑150 (MQ=255)
           tCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAggc                                                                                                                                          <  1:2988942/150‑1 (MQ=255)
           tCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAggc                                                                                                                                          >  1:941829/1‑150 (MQ=255)
           tCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAggc                                                                                                                                          >  1:1605593/1‑150 (MQ=255)
           tCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAggc                                                                                                                                          >  1:4215209/1‑150 (MQ=255)
           tCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAggc                                                                                                                                          >  1:3738536/1‑150 (MQ=255)
            ccAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAggcg                                                                                                                                         <  2:4669396/150‑1 (MQ=255)
               gCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGgcggcg                                                                                                                                      <  1:2334078/150‑1 (MQ=255)
                cccGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGt                                                                                                                                     >  2:4892270/1‑150 (MQ=255)
                 ccGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGtt                                                                                                                                    >  2:341380/1‑150 (MQ=255)
                 ccGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGt                                                                                                                                     <  2:3325841/149‑1 (MQ=255)
                 ccGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGt                                                                                                                                     >  1:3325841/1‑149 (MQ=255)
                 ccGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGCGCGAACCCAGGCGGCGtt                                                                                                                                    <  1:3997097/150‑1 (MQ=255)
                   ggAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTtgt                                                                                                                                  >  2:3826573/1‑150 (MQ=255)
                   ggAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTtgt                                                                                                                                  <  2:567146/150‑1 (MQ=255)
                     aGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTGGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGc                                                                                                                                <  2:3674053/150‑1 (MQ=255)
                         ggTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCAGTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGttt                                                                                                                            <  2:4586098/150‑1 (MQ=255)
                          gTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGtttt                                                                                                                           <  2:2417373/150‑1 (MQ=255)
                            gCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAATGCCGCCTGGGCGCGAACCCAGGCGGCGTtg                                                                                                                                   >  2:4405880/1‑140 (MQ=39)
                            gCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAATGCCGCCTGGGCGCGAACCCAGGCGGCGTtg                                                                                                                                   <  1:4405880/140‑1 (MQ=39)
                            gCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGctagccaa                                                                                                                         <  2:2381840/150‑8 (MQ=255)
                                      tgGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAg                                                                                                               <  2:2920682/150‑1 (MQ=255)
                                        gaaggaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGt                                                                                                             >  2:755365/1‑150 (MQ=255)
                                        gaaggaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGt                                                                                                             >  2:1234298/1‑150 (MQ=255)
                                          aggaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTgc                                                                                                           >  2:905308/1‑150 (MQ=255)
                                          aggaaATCGCTGCCAGGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTgc                                                                                                           >  1:303051/1‑150 (MQ=255)
                                                tCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCtt                                                                                                     >  2:5032287/1‑150 (MQ=255)
                                                tCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCtt                                                                                                     <  1:324216/150‑1 (MQ=255)
                                                 cGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCttt                                                                                                    <  1:2473795/150‑1 (MQ=255)
                                                    tGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCGTTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCa                                                                                                 >  2:1082864/1‑150 (MQ=255)
                                                    tGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCa                                                                                                 >  2:2221542/1‑150 (MQ=255)
                                                      ccATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACt                                                                                               >  1:1723141/1‑150 (MQ=255)
                                                       cATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTg                                                                                              <  2:3557780/150‑1 (MQ=255)
                                                               gggTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGAtt                                                                                      >  1:1874232/1‑150 (MQ=255)
                                                                  tGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGAATTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGc                                                                                   <  1:4653553/150‑1 (MQ=255)
                                                                    gagaACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCa                                                                                 >  1:874940/1‑150 (MQ=255)
                                                                      gaACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTgg                                                                                                                    >  2:4399761/1‑113 (MQ=255)
                                                                      gaACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTgg                                                                                                                    <  1:4399761/113‑1 (MQ=255)
                                                                       aaCACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCaata                                                                              <  1:2662552/150‑1 (MQ=255)
                                                                          acGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCATTGGCCAATCTGGCCagtgcaa                                                                           <  1:4892270/150‑7 (MQ=255)
                                                                          acGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCaataaat                                                                           >  2:3144118/1‑150 (MQ=255)
                                                                             gTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACTACGCCGCCTGGGTTCGCGCCCAGGCGTCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCGAATAAATACa                                                                        <  1:5010490/150‑1 (MQ=255)
                                                                               gACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGc                                                                      <  2:2401096/150‑1 (MQ=255)
                                                                                  ttcgttcgTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGa                                                                   >  1:1595286/1‑150 (MQ=255)
                                                                                     gttcgtCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTg                                                                >  2:5062118/1‑150 (MQ=255)
                                                                                        cgtCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCa                                                             <  1:899284/150‑1 (MQ=255)
                                                                                          tCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTgcgcgc                                                                                                       <  2:4338065/106‑1 (MQ=255)
                                                                                          tCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCata                                                           <  1:5084129/150‑1 (MQ=255)
                                                                                          tCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCata                                                           <  2:2200695/150‑1 (MQ=255)
                                                                                          tCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCGGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTgcgcgc                                                                                                       >  1:4338065/1‑106 (MQ=255)
                                                                                           ccTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCatat                                                          >  1:1317421/1‑150 (MQ=255)
                                                                                             tGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGCTTGGCCAATAAATACAGCGGATTGCCCTATTc                                                        >  1:3276342/1‑150 (MQ=255)
                                                                                                 cgcgGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACt                                                    >  2:3454121/1‑150 (MQ=255)
                                                                                                 cgcgGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACt                                                    >  2:3454329/1‑150 (MQ=255)
                                                                                                  gcgGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTg                                                   >  2:3169426/1‑150 (MQ=255)
                                                                                                   cgGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTgg                                                  <  2:4703454/150‑1 (MQ=255)
                                                                                                    ggtggtCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGc                                                 >  2:3121883/1‑150 (MQ=255)
                                                                                                           tGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGc                                                                                   >  1:4809399/1‑109 (MQ=255)
                                                                                                           tGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGc                                                                                   <  2:4809399/109‑1 (MQ=255)
                                                                                                           tGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCAc                                          <  2:874940/150‑1 (MQ=255)
                                                                                                           tGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCAc                                          >  2:3673799/1‑150 (MQ=255)
                                                                                                             aaCCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGc                                                                                   <  1:343849/107‑1 (MQ=255)
                                                                                                             aaCCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACgg                                        >  1:2680014/1‑150 (MQ=255)
                                                                                                                cgcgTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGctc                                     >  2:4790383/1‑150 (MQ=255)
                                                                                                                  cgTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGt                                                                                                             >  1:1487947/1‑76 (MQ=255)
                                                                                                                  cgTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGt                                                                                                             <  2:1487947/76‑1 (MQ=255)
                                                                                                                  cgTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCGACTGGCAGTGCACGGCTCtt                                   <  1:5377143/150‑1 (MQ=255)
                                                                                                                   gTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTa                                  >  2:2610948/1‑150 (MQ=255)
                                                                                                                    tCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTAc                                 <  2:267262/150‑1 (MQ=255)
                                                                                                                           aTTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTc                          <  1:607226/150‑1 (MQ=255)
                                                                                                                           aTTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTc                          >  2:3188175/1‑150 (MQ=255)
                                                                                                                            ttGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTcc                         <  1:3359760/150‑1 (MQ=255)
                                                                                                                              aaaacgCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTc                       >  1:4366/6‑150 (MQ=255)
                                                                                                                              aaaacgCACAACGCCGCCTGGGCGCGAACCCAGGCGGCGTTGTGCTTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTc                       >  1:3400684/6‑150 (MQ=38)
                                                                                                                               gCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCa                      >  2:4672005/1‑150 (MQ=255)
                                                                                                                                cTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCAt                     >  1:3486563/1‑150 (MQ=255)
                                                                                                                                 tAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCTTc                    <  1:2533812/150‑1 (MQ=255)
                                                                                                                                  aGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCg                   >  1:5340509/1‑150 (MQ=255)
                                                                                                                                     acaACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGAcc                >  2:1125771/1‑150 (MQ=255)
                                                                                                                                     acaACGCCGCCTGGGCGCGAACCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGAt                                                                  >  2:1054340/1‑100 (MQ=37)
                                                                                                                                     acaACGCCGCCTGGGCGCGAACCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGAt                                                                  <  1:1054340/100‑1 (MQ=37)
                                                                                                                                       aaCGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCa              >  1:3582721/1‑150 (MQ=255)
                                                                                                                                       aaCGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCa              >  1:3582617/1‑150 (MQ=255)
                                                                                                                                       aaCGCCGCCTGGGCGCGAACCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCa              >  2:165242/1‑150 (MQ=255)
                                                                                                                                            cgccTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGcc         >  2:2701817/1‑150 (MQ=255)
                                                                                                                                             gccTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGa                                                                   <  1:1018214/91‑1 (MQ=255)
                                                                                                                                             gccTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGa                                                                   >  2:1018214/1‑91 (MQ=255)
                                                                                                                                             gccTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCt        >  1:240276/1‑150 (MQ=255)
                                                                                                                                             gccTGGGCGCGAACCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCt        >  2:2099187/1‑150 (MQ=255)
                                                                                                                                              ccTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTg       >  2:4850770/1‑150 (MQ=255)
                                                                                                                                              ccTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTg       <  2:1294847/150‑1 (MQ=255)
                                                                                                                                                   gTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGACAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTggacag  >  2:4699241/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGTCGAGGCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTGGACAG  >  NC_002947/772033‑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: