Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,835,601 C→A 7.6% R111R (CGC→CGA PP_5114 → adenine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,835,6010CA 7.6% 242.5 / 4.1 77R111R (CGC→CGAPP_5114adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (34/37);  new base A (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 = 9.98e-01

CGCCCGCGTGCTGGACGCCTTCACCGGTAGTGGTGCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATGCTGGAGCAGAATCAGTGGCTGGGCGA  >  NC_002947/5835454‑5835747
                                                                                                                                                   |                                                                                                                                                   
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                                                                                                                                >  1:932098/1‑150 (MQ=255)
                   ttCACCGGTAGTGGTGCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGAc                                                                                                                                >  2:395981/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:32153/150‑1 (MQ=255)
                                            tCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCgg                                                                                                       >  2:203373/1‑150 (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                                                                                                              >  2:137461/1‑135 (MQ=255)
                                                    cTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCaa                                                                                                              <  1:137461/135‑1 (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:125502/1‑150 (MQ=255)
                                                         ccGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGtt                                                                                          >  1:136438/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                                                                                                           <  2:1046204/100‑1 (MQ=255)
                                                                                          cAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGc                                                                                                           >  1:1046204/1‑100 (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                                                                                                              >  2:2899/1‑75 (MQ=255)
                                                                                                                cTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCaa                                                                                                              <  1:2899/75‑1 (MQ=255)
                                                                                                                      aacaacCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATg                            <  2:1160379/150‑1 (MQ=255)
                                                                                                                      aacaacCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATg                            <  2:801356/150‑1 (MQ=255)
                                                                                                                        caacCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATGCt                          >  1:1373456/1‑150 (MQ=255)
                                                                                                                         aacCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCgt                                                                                <  2:1260907/95‑1 (MQ=255)
                                                                                                                         aacCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCgt                                                                                >  1:1260907/1‑95 (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)
                                                                                                                                                   |                                                                                                                                                   
CGCCCGCGTGCTGGACGCCTTCACCGGTAGTGGTGCCCTGGTGCTCGAAGCCCTGTCCCGTGGGGCCGAGGACGCCGTGGCGCTGGACAGCAACCCTGCGGCCATCGCCAACCTGAAGAACAACCTCGAAATCCTGCGCTGCCCTCGCGGGCAAATCCTGCAAACCGACGCCCTGCGCTACCTGCAAGGCCCGGCCAAGCAGCAGTTCGATGTCGTGTTCCTCGACCCGCCGTTCCATCAGGACCTGCTGGCCAACACCTGCAACATGCTGGAGCAGAATCAGTGGCTGGGCGA  >  NC_002947/5835454‑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: