Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,546,866 T→C 9.0% intergenic (+22/+18) PP_2238 → / ← PP_5504 metallopeptidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,546,8660TC 9.0% 218.2 / 9.5 78intergenic (+22/+18)PP_2238/PP_5504metallopeptidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base C (3/4);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

GAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGTGCTGCAGGCGCT  >  NC_002947/2546719‑2547015
                                                                                                                                                   |                                                                                                                                                      
gAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGc                                                                                                                                                      <  1:188163/150‑1 (MQ=255)
         aCCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTCGCAAACTGccc                                                                                                                                             <  2:1827986/150‑1 (MQ=255)
         aCCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTCGCAAACTGccc                                                                                                                                             <  2:508257/150‑1 (MQ=255)
          ccAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTCGCAAACTGCCCt                                                                                                                                            <  2:530212/150‑1 (MQ=255)
           cAAGGAGGCGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCtt                                                                                                                                           >  2:2135236/1‑150 (MQ=255)
           cAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCtt                                                                                                                                           >  1:1606745/1‑150 (MQ=255)
           cAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATCAATAGTTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCtt                                                                                                                                           <  1:1969521/150‑1 (MQ=255)
             aggaggaggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCtttt                                                                                                                                         <  2:463410/150‑1 (MQ=255)
              ggaggaggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGc                                                                                                                                                      >  1:848937/1‑136 (MQ=255)
              ggaggaggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGc                                                                                                                                                      <  2:848937/136‑1 (MQ=255)
               gaggaggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTtct                                                                                                                                       <  2:1758048/150‑1 (MQ=255)
                    ggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTAtt                                                                                                                                  <  2:996434/150‑1 (MQ=255)
                     gTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAAc                                                                                                                                                  >  2:391192/1‑133 (MQ=255)
                     gTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAAc                                                                                                                                                  <  1:391192/133‑1 (MQ=255)
                       ggAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTAc                                                                                                                               <  2:994805/150‑1 (MQ=255)
                           tGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgat                                                                                                                           <  2:1606745/150‑1 (MQ=255)
                           tGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGAACGGCACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgat                                                                                                                           >  1:1971852/1‑150 (MQ=255)
                             gCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgatga                                                                                                                         >  2:2057477/1‑150 (MQ=255)
                               tGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTGTTTTCAgttgatt                                                                                                                       >  2:2352357/1‑149 (MQ=255)
                                 aaCGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgatgatgat                                                                                                                     >  2:1576164/1‑150 (MQ=255)
                                  ccGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATc                                                                                                                    <  1:356526/149‑1 (MQ=255)
                                      cTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTAc                                                                                                                               >  1:1853998/1‑135 (MQ=255)
                                      cTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTAc                                                                                                                               <  2:1853998/135‑1 (MQ=255)
                                               cTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCt                                                                                                                                            <  2:859261/113‑1 (MQ=255)
                                               cTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCt                                                                                                                                            >  1:859261/1‑113 (MQ=255)
                                                     aCATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGc                                                                                                <  1:2135236/150‑1 (MQ=255)
                                                            cgcCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAAAGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCgg                                                                                         <  1:2306822/150‑1 (MQ=255)
                                                             gcCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCggg                                                                                        >  1:6307/1‑150 (MQ=255)
                                                               cTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGtt                                                                                      <  1:861111/150‑1 (MQ=255)
                                                                    gcCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCg                                                                                                           <  1:949543/124‑1 (MQ=255)
                                                                    gcCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCg                                                                                                           >  2:949543/1‑124 (MQ=255)
                                                                    gcCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTg                                                                                                       >  2:1119727/1‑128 (MQ=255)
                                                                    gcCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTg                                                                                                       <  1:1119727/128‑1 (MQ=255)
                                                                     cccTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGcc                                                                                >  1:1237612/1‑150 (MQ=255)
                                                                        ttGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGaa                                                                             >  2:2138701/1‑150 (MQ=255)
                                                                            tCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACaaa                                                                         >  2:318251/1‑150 (MQ=255)
                                                                            tCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACaaa                                                                         >  2:2014397/1‑150 (MQ=255)
                                                                              cAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCTCCGTGTTGATGCCGAACaaaga                                                                       >  2:1723423/1‑148 (MQ=255)
                                                                              cAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATa                                                                       >  1:1055394/1‑150 (MQ=255)
                                                                              cAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTCGCAAACTGcc                                                                                                                                              >  1:1646425/1‑80 (MQ=39)
                                                                              cAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTCGCAAACTGcc                                                                                                                                              <  2:1646425/80‑1 (MQ=39)
                                                                                   ggATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCAc                                                                  <  1:246921/150‑1 (MQ=255)
                                                                                    gATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCAcc                                                                 >  1:2070257/1‑150 (MQ=255)
                                                                                     aTGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCg                                                                <  2:1055394/150‑1 (MQ=255)
                                                                                         gcTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCct                                                            >  2:2209746/1‑150 (MQ=255)
                                                                                           tAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCctct                                                          >  1:809371/1‑150 (MQ=255)
                                                                                              ccTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCa                                                       <  1:2014397/150‑1 (MQ=255)
                                                                                                    tggAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTa                                                 >  2:412021/1‑150 (MQ=255)
                                                                                                     ggAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGtt                                                                                      >  1:2102634/1‑112 (MQ=255)
                                                                                                     ggAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGtt                                                                                      <  2:2102634/112‑1 (MQ=255)
                                                                                                        gaTGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCGGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGc                                                   <  2:1477254/144‑1 (MQ=255)
                                                                                                        gaTGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCGGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGc                                                   >  1:1477254/1‑144 (MQ=255)
                                                                                                             gTANTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAg                                        >  2:1497376/1‑150 (MQ=255)
                                                                                                               aCTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTa                                                                                                  >  1:1121841/1‑90 (MQ=255)
                                                                                                               aCTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTa                                                                                                  <  2:1121841/90‑1 (MQ=255)
                                                                                                                 tGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATa                                                                       <  2:984206/115‑1 (MQ=255)
                                                                                                                 tGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATa                                                                       >  1:984206/1‑115 (MQ=255)
                                                                                                                    aCCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCaa                                 >  1:2170918/1‑150 (MQ=255)
                                                                                                                         tGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGAcc                            >  1:912159/1‑150 (MQ=255)
                                                                                                                         tGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGAcc                            >  2:1178267/1‑150 (MQ=255)
                                                                                                                         tGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGAcc                            <  1:1218248/150‑1 (MQ=255)
                                                                                                                          gTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGAccc                           <  2:2322945/150‑1 (MQ=255)
                                                                                                                          gTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGAccc                           >  1:2322945/1‑150 (MQ=255)
                                                                                                                            aaCGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTATCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCgg                         <  1:34995/150‑1 (MQ=255)
                                                                                                                             aCGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTa                                                 <  1:1849106/125‑1 (MQ=255)
                                                                                                                             aCGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTa                                                 >  2:1849106/1‑125 (MQ=255)
                                                                                                                              cGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTa                       <  1:978391/150‑1 (MQ=255)
                                                                                                                               gCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCa                                                                                                                   >  1:2427680/1‑58 (MQ=255)
                                                                                                                               gCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCa                                                                                                                   <  2:2427680/58‑1 (MQ=255)
                                                                                                                                cAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGc                     <  1:1178267/150‑1 (MQ=255)
                                                                                                                                 aGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGcc                    <  2:1120598/150‑1 (MQ=255)
                                                                                                                                    gagaaaaGGGCAGTTCGCAAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCcgcg                 >  1:1601002/4‑150 (MQ=255)
                                                                                                                                        aaaGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATTATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGtg             >  1:503833/1‑150 (MQ=255)
                                                                                                                                         aaGGGCAGTTCGCAAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGtgc            >  2:384679/1‑150 (MQ=255)
                                                                                                                                           ggacAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGtgctg          >  2:1674103/4‑150 (MQ=255)
                                                                                                                                             gCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCCTTGTACTGGAGGGTGATCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGTGCTGCa        <  2:62815/150‑1 (MQ=255)
                                                                                                                                              cAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGTGCTGCAg       <  2:1971852/150‑1 (MQ=255)
                                                                                                                                                gTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGTTCTGCAGgc     >  2:381279/1‑150 (MQ=255)
                                                                                                                                                 tttGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGc                                                                                 >  1:1511747/1‑73 (MQ=255)
                                                                                                                                                 tttGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGc                                                                                 <  2:1511747/73‑1 (MQ=255)
                                                                                                                                                 tttGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGTGCTGCAGgcg    <  1:1674103/150‑1 (MQ=255)
                                                                                                                                                   tGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGTGCTGCAGGCGCt  <  2:2170918/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGTGCTGCAGGCGCT  >  NC_002947/2546719‑2547015

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