Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,904,208 C→T 8.7% intergenic (+287/‑233) PP_t18 → / → PP_1703 tRNA‑Asp/assimilatory nitrate reductase/sulfite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,904,2080CT 8.7% 919.4 / 51.8 311intergenic (+287/‑233)PP_t18/PP_1703tRNA‑Asp/assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base C (147/137);  new base T (22/5);  total (169/142)
Fisher's exact test for biased strand distribution p-value = 3.92e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 312 total aligned reads displayed.

TCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCA  >  NC_002947/1904061‑1904357
                                                                                                                                                   |                                                                                                                                                     
  aGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCC                                                                                                                                                   >  1:537646‑M1/1‑150 (MQ=255)
  aGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCC                                                                                                                                                   >  2:1311221‑M1/1‑150 (MQ=255)
  aGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCC                                                                                                                                                   <  1:260923‑M1/150‑1 (MQ=255)
       gccggctgccccggcatgctgatgcaagcggggcaacacggcgcctttgtcgtgaagtggtgagggctgtcgacgttgacagcccgggatgaggtagatcagcaatttatgggccagctataccggccctTTCGCGGGCACGCCCGCTcc                                                                                                                                              >  2:5938291‑M1/131‑150 (MQ=255)
       cTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCC                                                                                                                                              >  2:3015172‑M1/1‑150 (MQ=255)
         cggctgccccggcatgctgatgcaagcggggcaacacggcgcctttgtcgtgaagtggtgagggctgtcgacgttgacagcccgggatgaggtagatcagcaatttatgggccagctataccggcctcTTCGCGGGCACGCCCGCTCCca                                                                                                                                            >  1:4822595‑M1/129‑150 (MQ=255)
         cAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCc                                                                                                                                            >  1:5318954‑M1/1‑149 (MQ=255)
         cAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCA                                                                                                                                            >  2:2365918‑M1/1‑150 (MQ=255)
            aCCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAG                                                                                                                                         >  1:1528869‑M1/1‑150 (MQ=255)
            aCCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAG                                                                                                                                         >  2:4937462‑M1/1‑150 (MQ=255)
            aCCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAG                                                                                                                                         <  2:489195‑M1/150‑1 (MQ=255)
            aCCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAG                                                                                                                                         <  2:187454‑M1/150‑1 (MQ=255)
               cGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTa                                                                                                                                      <  2:6476557/150‑1 (MQ=255)
                   cTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTgc                                                                                                                                  <  1:5353775/150‑1 (MQ=255)
                    tCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTgcg                                                                                                                                 <  2:3981263/150‑1 (MQ=255)
                    tCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTgcg                                                                                                                                 >  1:3003046/1‑150 (MQ=255)
                    tCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACCCCCCCACCCACAGGTaaagcg                                                                                                                                 >  1:4083990/1‑145 (MQ=255)
                      gATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGcgcc                                                                                                                               <  1:702043/150‑1 (MQ=255)
                      gATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGcgcc                                                                                                                               >  2:4731390/1‑150 (MQ=255)
                      gATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGGTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGcgcc                                                                                                                               <  1:462077/150‑1 (MQ=255)
                        tCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCACCCACAGGGACTGcgccgc                                                                                                                             >  1:4113059/1‑150 (MQ=255)
                            cGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTc                                                                                                                         >  2:4251226/1‑150 (MQ=255)
                             ggCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACTGCGCCGCCCTCa                                                                                                                        <  2:2049416/150‑1 (MQ=255)
                              gCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAg                                                                                                                       >  2:3343094/1‑150 (MQ=255)
                                   cggggcaacacggcgcctttgtcgtgaagtggtgagggctgtcgacgttgacagcccgggatgaggtagatcagcaatttatgggccagctataccggccctTTCGCGGGCATGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCtt                                                                                                                  >  2:367114‑M1/103‑150 (MQ=255)
                                   cggggcaacacggcgcctttgtcgtgaagtggtgagggctgtcgacgttgacagcccgggatgaggtagatcagcaatttatgggccagctataccggccccTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCtt                                                                                                                  >  2:2640552‑M1/103‑150 (MQ=255)
                                    ccTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTg                                                                                                                 <  1:5852230/150‑1 (MQ=255)
                                     cTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGc                                                                                                                <  1:4296814/150‑1 (MQ=255)
                                      tGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCg                                                                                                               >  2:2397865/1‑150 (MQ=255)
                                      tGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCg                                                                                                               <  1:903987/150‑1 (MQ=255)
                                        ggACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCGCCCACAGGTACTGCGCCGCCCGCAGGCCTGGCGat                                                                                                             >  1:6324530/1‑150 (MQ=255)
                                         gACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGata                                                                                                            <  1:5332905/150‑1 (MQ=255)
                                          aCATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGgatatcgcaaggcctgagggcggcgcagt                                                                                                           <  1:1369866‑M1/150‑30 (MQ=255)
                                            aTGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCc                                                                                                         >  1:2636348‑M2/1‑149 (MQ=255)
                                            aTGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACCCCCGCTCCCACCCGTACTGCGCCGCCCGCAGGCCGTGCGATATCc                                                                                                         <  1:3450457‑M2/150‑2 (MQ=255)
                                               ggAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCcctg                                                                                                      <  2:6232415‑M2/150‑5 (MQ=255)
                                                gAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCcctgt                                                                                                     >  1:6153749‑M2/1‑145 (MQ=255)
                                                gAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCcctgt                                                                                                     >  1:991283‑M2/1‑145 (MQ=255)
                                                cgcctttgtcgtgatgtggtgagggctgtcgacgtggccagcccgggatgaggtagatcagcaatttatgggccagctataccggccctTTCGCGGGCATGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCtgt                                                                                                     >  1:5265203‑M1/90‑150 (MQ=255)
                                                cgcctttgtcgtgacgtggtgagggctgtcgacgttgacagcccgggatgaggtagatcagcaatttatgggccagctataccggccctTTCGCGGGCATGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCtgt                                                                                                     >  2:2426188‑M1/90‑150 (MQ=255)
                                                cgcctttgtcgtgaagtggtgagggctgtcgacgttgacagcccgggatgaggtagatcagcaatttatgggccagctataccggccctTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCtgt                                                                                                     >  2:6300788‑M1/90‑150 (MQ=255)
                                                  cctttgtcttgatgtggtgagggctgtcgacgttgacagcccgggatgaggtagatcagcaatttatgggccagctataccggccctTTCGCGGGCATGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTgg                                                                                                   >  1:4124300‑M1/88‑150 (MQ=255)
                                                  cctttgtcgtgaagtggtgagggctgtcgacgttgacagcccgggatgaggtagatcagcaatttatgggccagctataccggccctTTCGCGGGCATGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTgg                                                                                                   >  2:6549281‑M1/88‑150 (MQ=255)
                                                   ctttgtcgtgaagtggtgagggctgtcgacgttgacagcccgggatgaggtagatcagcaatttatgggccagctataccggccctTTCGCGGGCATGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTggg                                                                                                  >  2:2909555‑M1/87‑150 (MQ=255)
                                                    gCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCcctgtggga                                                                                                 >  2:1744172‑M2/1‑141 (MQ=255)
                                                     cAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGgatatcgcaaggcctgagggcggcgcagtacctgtgggag                                                                                                <  2:4303247‑M1/150‑41 (MQ=255)
                                                      aGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTTTTCGCGGGCATGCCCGCTCCCACAGGgatatcgcaaggcctgagggcggcgcagtacctgtgggagc                                                                                               <  2:4560899‑M1/150‑42 (MQ=255)
                                                      aGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCGCAGGCCTTGCGCTATCcctgtgggagc                                                                                               <  2:860184‑M2/150‑12 (MQ=255)
                                                      aGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCCCCCGCAGGCCGTGCGCTATCcctgtgggagc                                                                                               <  1:5969541‑M2/150‑12 (MQ=255)
                                                       gTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGgatatcgcaaggcctgagggcggcgcagtacctgtgggagcg                                                                                              <  1:1238251‑M1/150‑43 (MQ=255)
                                                        tgttCAGGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCgg                                                                                             <  1:439204/147‑1 (MQ=255)
                                                        tGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCcctgtgggagcgg                                                                                             >  1:6700239‑M2/1‑137 (MQ=255)
                                                          cTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCcctgtgggagcgggc                                                                                           <  2:3492186‑M2/150‑16 (MQ=255)
                                                            cATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGATGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCAt                                                                                         <  1:1033966/150‑1 (MQ=255)
                                                             aTGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATg                                                                                        >  1:2136419/1‑150 (MQ=255)
                                                                cTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTa                                                                                                                                      >  2:5469119/1‑101 (MQ=255)
                                                                cTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCTCCCTCAGGCCTTGCGCTATCCCTGTGGGAGCGGGCATGccc                                                                                     >  1:5466459/1‑150 (MQ=255)
                                                                 tGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCcg                                                                                    >  2:2255867/1‑150 (MQ=255)
                                                                      atgctcccAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCCCGCCCACACGGACAGCGCCCCCCTCACGCCTTGCGGAAACCCAGTGGGACCCGGCATGCCCCCCaaa                                                                               >  1:2897126/6‑150 (MQ=255)
                                                                      aTTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGTTACTGCGCCGCCCTCAGGCCTTGCTTTCTCCCTGTGGGAGCGGGCATGCCCGCGaaa                                                                               <  1:3974111/150‑1 (MQ=255)
                                                                      aTTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGaaa                                                                               <  2:2483301/150‑1 (MQ=255)
                                                                       ttGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCACAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAg                                                                              >  1:4750992/1‑150 (MQ=255)
                                                                       ttGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGGGCGTACCCTGTGCAGGCCTCGTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCGCCGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCcgctccca                                                                              <  1:1726728/150‑6 (MQ=38)
                                                                        tGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAgg                                                                             <  1:4869401/150‑1 (MQ=255)
                                                                          ccccACCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCGTCGCGGGCACGCCCGCCCCCACCGGTACAGCGCCGCCCTCAGGCCTTGCTATATCCCTGTGGGAGCGGGCAGGCCCGCGAAAGGGc                                                                           <  1:2341366/150‑1 (MQ=255)
                                                                           cccAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGGGCGGGCAATCCCGCGAAAGGGcc                                                                          >  1:1492803/1‑150 (MQ=255)
                                                                           cccAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGcc                                                                          <  2:317454/150‑1 (MQ=255)
                                                                            cggacAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCg                                                                         <  1:5707825/146‑1 (MQ=255)
                                                                            ccAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCACATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCg                                                                         >  1:4137557/1‑150 (MQ=255)
                                                                             cAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCgg                                                                        <  2:5394550/150‑1 (MQ=255)
                                                                                         tGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGccc                                                            >  1:3831909/1‑150 (MQ=255)
                                                                                          gggAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTCTAGCTGGCCCa                                                           >  2:6002465/1‑150 (MQ=255)
                                                                                          gggAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCGTGCGATATCCCGGTGGGCGCGGGCATGCCCCCGCAAGGGCCGGGGTAGCTGGCCCa                                                           >  1:4193131/1‑150 (MQ=38)
                                                                                             aGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATaa                                                        <  1:2757276/150‑1 (MQ=255)
                                                                                              ccATGCGGTGGCAACCTCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATaaa                                                       >  1:1055469/2‑150 (MQ=255)
                                                                                              ccATGCGGTGGCAACCTCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCAGaaa                                                       >  1:3920694/2‑150 (MQ=37)
                                                                                                aTGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAAtt                                                     >  1:4532645/1‑150 (MQ=255)
                                                                                                aTGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCGCCCACAGGTACTTCGCCGCCCTCAGGCCGTGCGATATCCCGGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCAGaaagt                                                     >  1:5523929/1‑148 (MQ=38)
                                                                                                 tGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCGCCCCCCGGTACAGCGCCGCCCTCAGGCCGGGCGATATCCCGGGGGGAGCGGGCATGCCCCCGCAAGGGCCGGGAAAGCTGGCCCAAAAAATg                                                    >  1:1817640/1‑150 (MQ=255)
                                                                                                  gAGGTTGCCACCGCTTGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGTTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGc                                                   <  1:2000402/150‑1 (MQ=255)
                                                                                                  gAGGTTGCCACCGCTTGGTGCGTAACCTGTGAAGGCCTCTTCGCGGGCACGCCCGCTCCCACTGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGc                                                   <  1:3146715/150‑1 (MQ=255)
                                                                                                  gAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGc                                                   <  2:5850010/150‑1 (MQ=255)
                                                                                                          cACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTAc                                           >  1:5722571/1‑150 (MQ=255)
                                                                                                            ccGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACAGCGCCGCCCACAGGCCGTGCGATAACCCTGTGGGAGCGGGCCTGCCCCCCAAAGGGCCCGTATAGCTGGCCCCTAAAATTCTGAtaaacca                                         >  1:4960510/1‑144 (MQ=37)
                                                                                                             cGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTc                                        <  1:375201/150‑1 (MQ=255)
                                                                                                               cATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCAt                                      >  1:6820533/1‑150 (MQ=255)
                                                                                                                     gCGTAACCTGTGCAGGCCTCTTCGCAGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCggg                                <  2:5251104/150‑1 (MQ=255)
                                                                                                                               tgCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACTGCGCCGCCCTCAGGCCTTGCGATATCCCGGGGGGAGCGGGCATGCCCGCGAAAGGGCCGGGATAGCTGGCCCAGAAATTGCTGATCTACCTCATCCCGGGCTGTCCACgg                      >  2:3555801/1‑149 (MQ=38)
                                                                                                                                 cAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCg                    <  2:1154215/150‑1 (MQ=255)
                                                                                                                                  aGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCGTGCGAGATCCCGGGGGGAGCGGGCATGCCCGCCAAAGGGCCGGGATAGCGGGCCCATAAATTGCTGATCTACCGCATCCCCGGCGGGCAACGGCgc                   >  1:3318058/1‑149 (MQ=37)
                                                                                                                                  aGGCCTCTTCGCGGGCACGCCCGCGCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGa                   <  1:5163034/150‑1 (MQ=255)
                                                                                                                                   ggCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGAc                  >  2:4454658/1‑150 (MQ=255)
                                                                                                                                        tTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGccc             <  1:5616514‑M1/149‑1 (MQ=255)
                                                                                                                                        ctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCcctgtgggagcggg                                                                                            <  2:1341979‑M2/71‑15 (MQ=255)
                                                                                                                                        ctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCGACGTTGACAGccc             <  2:1817640/150‑1 (MQ=255)
                                                                                                                                         TTCGCGGGCACGCCCGCTCCCACCGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGGGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCGACGTTGACAgcccg            >  1:2493492‑M1/1‑149 (MQ=255)
                                                                                                                                          TCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTc           >  2:6738567‑M1/1‑150 (MQ=255)
                                                                                                                                            GCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTcac         >  2:1518117‑M1/1‑150 (MQ=255)
                                                                                                                                            GCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTcac         >  1:3343094‑M1/1‑150 (MQ=255)
                                                                                                                                                   cGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCa  <  1:6254989/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCA  >  NC_002947/1904061‑1904357

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