Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,835,598 T→G 7.8% P110P (CCT→CCG PP_5114 → adenine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,835,5980TG 7.8% 206.5 / 6.8 77P110P (CCT→CCGPP_5114adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (34/37);  new base G (4/2);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

AGGGCGCCCGCGTGCTGGACGCCTTCACCGGTAGTGGTGCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATGCTGGAGCAGAATCAGTGGCTGGGCGA  >  NC_002947/5835450‑5835747
                                                                                                                                                    |                                                                                                                                                      
aGGGCGCCCGCGTGCTGGACGCCTTCACCGGTAGTGGTGCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTc                                                                                                                                                       >  2:19367/1‑150 (MQ=255)
    cgcCCGCGTGCTGGACGCCTTCACCGGTAGTGGTGCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCgg                                                                                                                                                   <  2:945568/150‑1 (MQ=255)
     gcCCGCGTGCTGGACGCCTTCACCGGTAGTGGTGCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCggg                                                                                                                                                  <  1:438188/150‑1 (MQ=255)
              cTGGACGCCTTCACCGGTAGTGGTGCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTg                                                                                                                                         <  2:805537/150‑1 (MQ=255)
                      cTTCACCGGTAGTGGTGCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGa                                                                                                                                 >  1:865025/1‑150 (MQ=255)
                      cTTCACCGGTAGTGGTGCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGa                                                                                                                                 >  2:72748/1‑150 (MQ=255)
                       ttCACCGGTAGTGGTGCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGAc                                                                                                                                >  2:395981/1‑150 (MQ=255)
                       ttCACCGGTAGTGGTGCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGAc                                                                                                                                >  1:932098/1‑150 (MQ=255)
                         cACCGGTAGTGGTGCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCGGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGc                                                                                                                              >  2:463826/1‑150 (MQ=255)
                         cACCGGTAGTGGTGCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGc                                                                                                                              <  1:81686/150‑1 (MQ=255)
                            cGGTAGTGGTGCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCGCGAGGGCagcgc                                                                                                                                            <  1:144576/133‑5 (MQ=255)
                            cGGTAGTGGTGCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCGCGAGGGCagcgc                                                                                                                                            >  2:144576/1‑129 (MQ=255)
                             ggTAGTGGTGCCCTGGTGCGCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCCCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGcccag                                                                                                                          >  2:510278/1‑148 (MQ=255)
                              gTAGTGGTGCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTgc                                                                                                                         <  2:932098/150‑1 (MQ=255)
                                 gtggtgCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAACTCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCt                                                                                                                      <  1:741195/150‑1 (MQ=255)
                                  tggtgCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCCCGGGCAAATCCTGCAAACCGACGCCCTGCGCTa                                                                                                                     >  2:66279/1‑150 (MQ=255)
                                         cTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCaa                                                                                                              >  2:451181/1‑150 (MQ=255)
                                          tGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACTCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAg                                                                                                             >  1:238943/1‑150 (MQ=255)
                                                tCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCgg                                                                                                       >  2:203373/1‑150 (MQ=255)
                                                tCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCgg                                                                                                       <  2:32153/150‑1 (MQ=255)
                                                  gAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCgg                                                                                                       <  2:1303680/148‑1 (MQ=255)
                                                  gAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCgg                                                                                                       >  1:1303680/1‑148 (MQ=255)
                                                  gAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGcc                                                                                                     >  1:841682/1‑150 (MQ=255)
                                                    aGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCaa                                                                                                   <  2:289131/150‑1 (MQ=255)
                                                      cccTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAagc                                                                                                 >  1:641367/1‑150 (MQ=255)
                                                       ccTGTCCCGTGGGGCCGAGGACTCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAagca                                                                                                <  2:238943/150‑1 (MQ=255)
                                                       ccTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTg                                                                                                                                         >  1:296624/1‑109 (MQ=255)
                                                       ccTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTg                                                                                                                                         <  2:296624/109‑1 (MQ=255)
                                                        cTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCaa                                                                                                              <  1:137461/135‑1 (MQ=255)
                                                        cTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCaa                                                                                                              >  2:137461/1‑135 (MQ=255)
                                                        cTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAg                                                                                                             <  2:1381864/136‑1 (MQ=255)
                                                        cTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAg                                                                                                             >  1:1381864/1‑136 (MQ=255)
                                                         tGTCCCGTGGGGCCCAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCCCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAGACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAagcagc                                                                                              <  1:1359209/150‑1 (MQ=255)
                                                             ccGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGtt                                                                                          >  1:136438/1‑150 (MQ=255)
                                                             ccGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGtt                                                                                          >  1:125502/1‑150 (MQ=255)
                                                               gTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGc                                                                                                      >  1:1202328/1‑136 (MQ=255)
                                                               gTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGc                                                                                                      <  2:1202328/136‑1 (MQ=255)
                                                                 ggggCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGAt                                                                                      >  1:1160379/1‑150 (MQ=255)
                                                                    gCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTc                                                                                   <  2:136438/150‑1 (MQ=255)
                                                                      cGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCgt                                                                                <  2:679573/150‑1 (MQ=255)
                                                                       gAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCgtg                                                                               <  2:1053088/150‑1 (MQ=255)
                                                                          gACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTc                                                                            >  2:540716/1‑150 (MQ=255)
                                                                             gCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAAGCACGGAAAAGAAGCCCGTGGATCTCgtgtccgtg                                                                         >  2:1014498/1‑145 (MQ=255)
                                                                                       tGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCACTCCCGGGCAATGCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGt                                                               <  2:44707/150‑1 (MQ=255)
                                                                                            aGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCAt                                                          <  1:400019/150‑1 (MQ=255)
                                                                                              cAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGc                                                                                                           >  1:1046204/1‑100 (MQ=255)
                                                                                              cAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGc                                                                                                           <  2:1046204/100‑1 (MQ=255)
                                                                                                      cGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACctgc                                                <  2:641367/150‑1 (MQ=255)
                                                                                                           aTCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGcc                                           <  2:736532/150‑1 (MQ=255)
                                                                                                           aTCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGcc                                           <  1:203373/150‑1 (MQ=255)
                                                                                                           aTCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCCCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGcc                                           <  1:131780/150‑1 (MQ=255)
                                                                                                            tCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCaa                                                                                                              >  2:931583/1‑83 (MQ=255)
                                                                                                            tCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCaa                                                                                                              <  1:931583/83‑1 (MQ=255)
                                                                                                                cAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAgcagca                                                                                             >  2:1423991/1‑96 (MQ=255)
                                                                                                                cAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAgcagca                                                                                             <  1:1423991/96‑1 (MQ=255)
                                                                                                                    cTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCaa                                                                                                              <  1:2899/75‑1 (MQ=255)
                                                                                                                    cTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCaa                                                                                                              >  2:2899/1‑75 (MQ=255)
                                                                                                                          aacaacCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATg                            <  2:1160379/150‑1 (MQ=255)
                                                                                                                          aacaacCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATg                            <  2:801356/150‑1 (MQ=255)
                                                                                                                            caacCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATGCt                          >  1:1373456/1‑150 (MQ=255)
                                                                                                                             aacCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCgt                                                                                >  1:1260907/1‑95 (MQ=255)
                                                                                                                             aacCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCgt                                                                                <  2:1260907/95‑1 (MQ=255)
                                                                                                                               ccTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATGCTGGa                       <  1:395981/150‑1 (MQ=255)
                                                                                                                               ccTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATGCTGGa                       <  1:72748/150‑1 (MQ=255)
                                                                                                                                cTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGCCGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGCCCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATGCTGGAg                      <  2:26302/150‑1 (MQ=255)
                                                                                                                                cTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTTCGCGTGCTGTAGGGCCCGGGAAAGCACCAGGGGGGCGCGGCTGAGCCCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTTCGGCATTCTGGAg                      >  1:679088/1‑150 (MQ=255)
                                                                                                                                 tCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATGCTGGAGc                     <  2:125502/150‑1 (MQ=255)
                                                                                                                                   gAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATGCTGGAGCAg                   >  2:539255/1‑150 (MQ=255)
                                                                                                                                       tCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATGCTGGAGCAGAATc               <  1:19367/150‑1 (MQ=255)
                                                                                                                                          tGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATGCTGGAGCAGAATCAGt            <  1:451181/150‑1 (MQ=255)
                                                                                                                                           gcgcTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATGCTGGAGCAGAATCAGTg           >  2:173125/1‑150 (MQ=255)
                                                                                                                                            atttGCCCGCGAGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATGCTGGAGCAGAATCAGTgg          >  1:1328641/4‑150 (MQ=255)
                                                                                                                                             tttGCCCGCGAGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATTCTGGAGCAGAATCAGTGGc         >  2:187298/3‑150 (MQ=255)
                                                                                                                                             tttGCCCGCGAGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATGCTg                         >  2:1323365/3‑134 (MQ=255)
                                                                                                                                             tttGCCCGCGAGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATGCTg                         <  1:1323365/132‑1 (MQ=255)
                                                                                                                                                 cccTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATTAGGACCTGCTCGCCAACACCTGCAACATGCTGGAGCAGAATCAGTGGCTggg     >  1:808372/1‑150 (MQ=255)
                                                                                                                                                  ccTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATGCTGGAGCAGAATCAGTGGCTGGGc    >  1:260848/1‑150 (MQ=255)
                                                                                                                                                    tCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATGCTGGAGCAGAATCAGTGGCTGGGCGa  >  1:661484/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
AGGGCGCCCGCGTGCTGGACGCCTTCACCGGTAGTGGTGCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATGCTGGAGCAGAATCAGTGGCTGGGCGA  >  NC_002947/5835450‑5835747

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