Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 772,181 T→C 7.7% 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 7.7% 897.8 / 41.7 287intergenic (+39/‑160)PP_0663/PP_0664AsnC family transcriptional regulator/homoserine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (114/150);  new base C (12/10);  total (126/161)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 291 total aligned reads displayed.

GCGTCGAGGCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTGGACAG  >  NC_002947/772033‑772329
                                                                                                                                                    |                                                                                                                                                    
   tCGAGGCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTcgc                                                                                                                                                  <  2:6915742/150‑1 (MQ=255)
   tCGAGGCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTcgc                                                                                                                                                  <  2:5908565/150‑1 (MQ=255)
   tCGAGGCCTCCAGCCCGGAGCGGGTGCACAATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTcgc                                                                                                                                                  >  2:2941093/1‑150 (MQ=255)
        gCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCa                                                                                                                                             >  1:4066852/1‑150 (MQ=255)
        gCCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCa                                                                                                                                             <  1:6850514/150‑1 (MQ=255)
          cTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAgg                                                                                                                                           >  2:2382941/1‑150 (MQ=255)
          cTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGATTCGCGCCCAgg                                                                                                                                           >  1:429454/1‑150 (MQ=255)
           tCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAggc                                                                                                                                          <  1:3599988/150‑1 (MQ=255)
           tCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAggc                                                                                                                                          >  1:5049222/1‑150 (MQ=255)
           tCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGCGCGAACCCAggc                                                                                                                                          >  1:2794646/1‑150 (MQ=255)
            cccgCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGCGCGAACCCAggcg                                                                                                                                         <  2:5331364/147‑1 (MQ=255)
            ccAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGCGCGAACCCAggcg                                                                                                                                         <  1:2077816/150‑1 (MQ=255)
               gCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGgcggcg                                                                                                                                      >  1:1263649/1‑150 (MQ=255)
                 ccGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGtt                                                                                                                                    <  2:2866158/150‑1 (MQ=255)
                 ccGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGtt                                                                                                                                    >  2:6450040/1‑150 (MQ=255)
                 ccGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGtt                                                                                                                                    >  1:1872062/1‑150 (MQ=255)
                   ggAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTtgt                                                                                                                                  >  2:4523382/1‑150 (MQ=255)
                          gTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGtttt                                                                                                                           >  2:2230030/1‑150 (MQ=255)
                            gCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGc                                                                                                                         <  1:7344962/150‑1 (MQ=255)
                              aCCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGccc                                                                                                                       >  1:2461865/1‑150 (MQ=255)
                              aCCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGccc                                                                                                                       >  2:941952/1‑150 (MQ=255)
                                 aTGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTgg                                                                                                                    <  2:3100063/150‑1 (MQ=255)
                                 aTGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTgg                                                                                                                    >  1:590864/1‑150 (MQ=255)
                                      tgGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAg                                                                                                               <  1:3040577/150‑1 (MQ=255)
                                         aaggaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTg                                                                                                            >  2:4322486/1‑150 (MQ=255)
                                          aggaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTgc                                                                                                           >  1:2344505/1‑150 (MQ=255)
                                               aTCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCt                                                                                                      <  2:2539080/150‑1 (MQ=255)
                                                 cGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCttt                                                                                                    <  1:3793719/150‑1 (MQ=255)
                                                 cGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCttt                                                                                                    <  2:5617709/150‑1 (MQ=255)
                                                  gCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTg                                                                                                   >  2:6707381/1‑150 (MQ=255)
                                                  gCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTg                                                                                                   <  1:6899675/150‑1 (MQ=255)
                                                    tGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCa                                                                                                 <  2:1113572/150‑1 (MQ=255)
                                                    tGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCa                                                                                                 <  2:6151451/150‑1 (MQ=255)
                                                    tGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCa                                                                                                 <  2:1976753/150‑1 (MQ=255)
                                                    tGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCa                                                                                                 >  1:600101/1‑150 (MQ=255)
                                                    tGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCa                                                                                                 >  2:5187111/1‑150 (MQ=255)
                                                      ccATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACt                                                                                               >  1:5635901/1‑150 (MQ=255)
                                                       cATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTg                                                                                              <  1:5740454/150‑1 (MQ=255)
                                                            ccGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAg                                                                                         <  2:5969070/150‑1 (MQ=255)
                                                               gggTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGAATTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGAtt                                                                                      <  2:2010751/150‑1 (MQ=255)
                                                                ggTGGAGAACCCGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTg                                                                                     <  2:7115182/150‑1 (MQ=255)
                                                                  tGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACCACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGc                                                                                   <  2:41427/150‑1 (MQ=255)
                                                                  tGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGc                                                                                   <  1:1098117/150‑1 (MQ=255)
                                                                    gagaACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCCGATTGGCCa                                                                                 >  1:1029047/1‑150 (MQ=255)
                                                                    gagaACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCa                                                                                 <  2:1317866/150‑1 (MQ=255)
                                                                       aaCACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCaata                                                                              <  2:2797942/150‑1 (MQ=255)
                                                                       aaCACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAATCTGGCCagtg                                                                              <  2:1509101/150‑4 (MQ=255)
                                                                        acacGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCaataa                                                                             <  2:4838217/150‑1 (MQ=255)
                                                                         cacGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCaataaa                                                                            <  2:4274451/150‑1 (MQ=255)
                                                                         cacGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCaataaa                                                                            <  2:1449522/150‑1 (MQ=255)
                                                                            ggTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATAc                                                                         >  1:6176652/1‑150 (MQ=255)
                                                                                  ttcgttcgTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACCGCGGa                                                                   <  1:5333209/150‑1 (MQ=255)
                                                                                  ttcgttcgTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGa                                                                   >  2:510419/1‑150 (MQ=255)
                                                                                    cgttcgtCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGAtt                                                                 >  1:981072/1‑150 (MQ=255)
                                                                                     gttcgtCCTGCCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCAGGGTTAGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCGTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTg                                                                >  2:3734781/1‑150 (MQ=255)
                                                                                       tcgtCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGcc                                                              <  2:1197049/150‑1 (MQ=255)
                                                                                            cTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATAtt                                                         >  1:6314363/1‑150 (MQ=255)
                                                                                             tGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATAtt                                                         >  1:724983/1‑149 (MQ=255)
                                                                                             tGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTc                                                        >  1:2912243/1‑150 (MQ=255)
                                                                                             tGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTc                                                        >  2:392671/1‑150 (MQ=255)
                                                                                                    ggtggtCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGc                                                 <  1:872346/150‑1 (MQ=255)
                                                                                                        gtCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTg                                             >  2:6191668/1‑150 (MQ=255)
                                                                                                           tGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGTGGATTGCCATATTCAACTGGCAGTGCAc                                          <  2:201700/150‑1 (MQ=255)
                                                                                                           tGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCAc                                          <  2:5731378/150‑1 (MQ=255)
                                                                                                           tGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCAc                                          <  1:5187111/150‑1 (MQ=255)
                                                                                                             aaCCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACgg                                        >  1:1997022/1‑150 (MQ=255)
                                                                                                               ccGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGct                                      >  2:4473875/1‑150 (MQ=255)
                                                                                                                cgcgTCTTTGCATTGGCTAGCACAAAGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCATGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGcta                                     >  1:427493/1‑149 (MQ=255)
                                                                                                                cgcgTCTTTGCATTGGCTAGCACAAAGCCGCCTGGGTTCGCGCCCAGGCGGCGATGTGCGTTATGCCCTGCCGCAGGTGCGCGCTTTGCACAGGCCAGATTGGCCAATAAATACAGCGGATTGCCATAATCAACTGGCAGTGCACGGctc                                     >  1:1568363/1‑150 (MQ=255)
                                                                                                                    tCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTAc                                 <  1:5594355/150‑1 (MQ=255)
                                                                                                                        tGCATTGGCTAGCACACCGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTa                             <  1:4701944/150‑1 (MQ=255)
                                                                                                                        tGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTa                             <  2:3620339/150‑1 (MQ=255)
                                                                                                                        tGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTa                             >  1:4513376/1‑150 (MQ=255)
                                                                                                                        tGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTa                             <  2:6232194/150‑1 (MQ=255)
                                                                                                                         gCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAg                            <  2:3569385/150‑1 (MQ=255)
                                                                                                                            ttGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTcc                         >  2:6875770/1‑150 (MQ=255)
                                                                                                                               gCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCa                      <  1:5452292/150‑1 (MQ=255)
                                                                                                                                  cgCACAACGCCGCCTGGGCGCGAACCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCg                   >  2:6308456/2‑150 (MQ=255)
                                                                                                                                   gCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGa                  >  1:6399995/1‑150 (MQ=255)
                                                                                                                                    cacaACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGcc                                                                                  <  1:2714250/85‑1 (MQ=255)
                                                                                                                                    cacaACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGcc                                                                                  >  2:2714250/1‑85 (MQ=255)
                                                                                                                                     acaACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGAcc                >  2:7327896/1‑150 (MQ=255)
                                                                                                                                     acaACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGAcc                >  1:2535173/1‑150 (MQ=255)
                                                                                                                                     acaACGCCGCCTGGGCGCGAACCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCAGCGAcc                >  2:6515455/1‑150 (MQ=39)
                                                                                                                                       aaCGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCa              >  1:4807225/1‑150 (MQ=255)
                                                                                                                                       aaCGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGAGTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCa              >  1:1071464/1‑150 (MQ=255)
                                                                                                                                       aaCGCCGCCTGGGCGCGAACCCAGGCGGCGTTGTGCGGTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCCGAGTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCa              >  2:4897031/1‑150 (MQ=39)
                                                                                                                                        aCGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCAc             <  1:5117174/150‑1 (MQ=255)
                                                                                                                                          gccgccTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTg           >  2:6014210/1‑150 (MQ=255)
                                                                                                                                          gccgccTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTg           >  2:3265736/1‑150 (MQ=255)
                                                                                                                                          gccgccTGGGCGCGAACCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTCCGCTAGTCCTCATCGACCCACTg           >  2:6343378/1‑150 (MQ=255)
                                                                                                                                          gccgccTGGGCGCGAACCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACt            <  2:536907/149‑1 (MQ=255)
                                                                                                                                          gccgccTGGGCGCGAACCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACt            >  1:536907/1‑149 (MQ=255)
                                                                                                                                           ccgccTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGc          <  2:1258801/150‑1 (MQ=255)
                                                                                                                                             gccTGGGCGCGAACCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCt        >  1:3285818/1‑150 (MQ=255)
                                                                                                                                              ccTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTg       <  1:5848434/150‑1 (MQ=255)
                                                                                                                                              ccTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTg       <  2:7175491/150‑1 (MQ=255)
                                                                                                                                                  ggTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTggaca   <  2:1071464/150‑1 (MQ=255)
                                                                                                                                                  ggTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTggaca   >  1:1416086/1‑150 (MQ=255)
                                                                                                                                                    ttCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCaata                                                                              <  1:6460232/73‑1 (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: