Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,484,771 A→T 8.1% intergenic (+371/‑98) PP_3976 → / → PP_3977 isochorismatase family hydrolase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,484,7710AT 8.1% 265.7 / 12.4 86intergenic (+371/‑98)PP_3976/PP_3977isochorismatase family hydrolase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (38/41);  new base T (3/4);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.39e-01

TCAAACCCTGGCAGGTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCC  >  NC_002947/4484624‑4484920
                                                                                                                                                   |                                                                                                                                                     
tCAAACCCTGGCAGGTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGCCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATc                                                                                                                                                     >  1:1132834/1‑150 (MQ=255)
  aaaCCCTGGCAGGTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAg                                                                                                                                                   <  2:1299866/150‑1 (MQ=255)
  aaaCCCTGGCAGGTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAg                                                                                                                                                   <  1:938732/150‑1 (MQ=255)
      ccTGGCAGGTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGGGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCt                                                                                                                                               >  2:333638/1‑150 (MQ=255)
      ccTGGCAGGTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCt                                                                                                                                               <  1:2388708/150‑1 (MQ=255)
       cTGGCAGGTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGATTATTCAGGGCtt                                                                                                                                              <  2:2083657/150‑1 (MQ=255)
            aGGTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAat                                                                                                                                         <  1:2007440/150‑1 (MQ=255)
               tCGCAACTGCCAAGCAGATTATGTTGGCGGCAAGTCGGTCACAGCATGCGGAAATGAATATAGAAATGCGGGCGTGATCCGTATGAAAAAGGTAGAAAATCTGAGCCCCAAAAACGCAAAAAGCCCTAAATAATCAGGGCTTTGAATATg                                                                                                                                      >  2:833526/1‑150 (MQ=37)
               tCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATg                                                                                                                                      <  2:181858/150‑1 (MQ=255)
                  cAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAgg                                                                                                                                                  >  2:714854/1‑135 (MQ=255)
                  cAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAgg                                                                                                                                                  <  1:714854/135‑1 (MQ=255)
                   aaCGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCgg                                                                                                                                  >  1:603097/1‑150 (MQ=255)
                   aaCGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCgg                                                                                                                                  >  1:856551/1‑150 (MQ=255)
                     cGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGaa                                                                                                                                <  2:818803/150‑1 (MQ=255)
                      ggCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAg                                                                                                                               >  1:1408612/1‑150 (MQ=255)
                         cAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGATTATTCAGGGCtttt                                                                                                                                            >  1:1529266/1‑133 (MQ=255)
                         cAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGATTATTCAGGGCtttt                                                                                                                                            <  2:1529266/134‑2 (MQ=255)
                         cAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGt                                                                                                                            <  1:1297566/150‑1 (MQ=255)
                          aaGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGc                                                                                                                              >  1:456264/1‑147 (MQ=255)
                          aaGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGc                                                                                                                              <  2:456264/147‑1 (MQ=255)
                            gccgatcatGTTGGCGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTaga                                                                                                                         >  2:2089122/4‑150 (MQ=255)
                              aGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAgaga                                                                                                                       <  1:1479841/150‑1 (MQ=255)
                                   atgTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGCTTCac                                                                                                                  <  2:603097/150‑2 (MQ=255)
                                    tgTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCaaa                                                                                                                 <  1:2057299/150‑1 (MQ=255)
                                      ttGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAAcc                                                                                                               >  1:1782707/1‑150 (MQ=255)
                                              atatCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCa                                                                                                                                                    >  1:939614/1‑105 (MQ=255)
                                              atatCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCa                                                                                                                                                    <  2:939614/105‑1 (MQ=255)
                                                 tCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAgg                                                                                                    >  1:1742196/1‑150 (MQ=255)
                                                    gTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCaaa                                                                                                                 >  1:1475636/1‑134 (MQ=255)
                                                    gTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCaaa                                                                                                                 <  2:1475636/134‑1 (MQ=255)
                                                      cacaGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTcct                                                                                               <  2:1509439/150‑1 (MQ=255)
                                                         aGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTcctacc                                                                                            <  2:1813823/150‑1 (MQ=255)
                                                           cATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTcctaccta                                                                                          >  2:739134/1‑150 (MQ=255)
                                                           cATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTcctaccta                                                                                          >  1:2129032/1‑150 (MQ=255)
                                                             tGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTTAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACa                                                                                        >  1:838225/1‑150 (MQ=255)
                                                               cGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGa                                                                                      <  1:1755947/150‑1 (MQ=255)
                                                                  aaaTGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACtt                                                                                   >  2:752944/1‑150 (MQ=255)
                                                                     tGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGGGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAAt                                                                                <  1:333638/150‑1 (MQ=255)
                                                                     tGAACAGAGAAAGGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAAt                                                                                <  2:834431/150‑1 (MQ=255)
                                                                       aaCAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTTAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAAttt                                                                              <  2:838225/150‑1 (MQ=255)
                                                                          agagAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAGAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTg                                                                           >  1:203231/1‑150 (MQ=255)
                                                                          agagAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTg                                                                           >  1:368813/1‑150 (MQ=255)
                                                                                     ggCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTc                                                                <  2:856551/150‑1 (MQ=255)
                                                                                           aTCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTc                                                          >  2:1671911/1‑150 (MQ=255)
                                                                                            tCCGTATGACAAGGGTAGAAAATCGGAGCCTCAAAAACGCAAAAAGGCCTGAATAATCAGGGCGTTGAATATGGCGGAAGCTTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTAATCGGTAGCTCGACGccc                                                         <  2:660017/150‑1 (MQ=255)
                                                                                            tCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTcc                                                         <  2:1742196/150‑1 (MQ=255)
                                                                                               gTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCttt                                                      >  1:596262/1‑150 (MQ=255)
                                                                                               gTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCGAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCttt                                                      >  1:2104837/1‑150 (MQ=255)
                                                                                                tATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTa                                                                                                                           <  1:1453646/80‑1 (MQ=255)
                                                                                                tATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTa                                                                                                                           >  2:1453646/1‑80 (MQ=255)
                                                                                                tATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTa                                                     >  1:2122266/1‑150 (MQ=255)
                                                                                                 aTGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTAt                                                    <  2:2104837/150‑1 (MQ=255)
                                                                                                      aaGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACt                                               <  2:2129032/150‑1 (MQ=255)
                                                                                                        gggTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATc                                                                                                         >  1:359009/1‑90 (MQ=255)
                                                                                                        gggTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATc                                                                                                         <  2:359009/90‑1 (MQ=255)
                                                                                                        gggTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGt                                             >  1:1822880/1‑150 (MQ=255)
                                                                                                         ggTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTa                                            >  1:469204/1‑150 (MQ=255)
                                                                                                          gTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTaa                                           >  2:1579351/1‑150 (MQ=255)
                                                                                                           tAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAg                                          <  2:1782707/150‑1 (MQ=255)
                                                                                                           tAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAg                                          <  1:2089122/150‑1 (MQ=255)
                                                                                                            aGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGc                                         <  2:596262/150‑1 (MQ=255)
                                                                                                             gAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCt                                        <  2:2122266/150‑1 (MQ=255)
                                                                                                               aaaTCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAACATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAg                                      <  2:734803/150‑1 (MQ=255)
                                                                                                                        ccccAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCgg                             <  1:739134/150‑1 (MQ=255)
                                                                                                                        ccccAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCgg                             <  1:804032/150‑1 (MQ=255)
                                                                                                                         cccAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTcctac                                                                                             <  1:1325640/85‑1 (MQ=255)
                                                                                                                         cccAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTcctac                                                                                             >  2:1325640/1‑85 (MQ=255)
                                                                                                                         cccAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCggg                            <  1:752944/150‑1 (MQ=255)
                                                                                                                          ccAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGt                           >  2:461656/1‑150 (MQ=255)
                                                                                                                          ccAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGt                           >  2:1239707/1‑150 (MQ=255)
                                                                                                                            aaaaaCGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGa                                                 >  1:1075822/1‑126 (MQ=255)
                                                                                                                            aaaaaCGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGa                                                 <  2:1075822/126‑1 (MQ=255)
                                                                                                                                    aaaaaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACt                                               <  2:371577/120‑1 (MQ=255)
                                                                                                                                    aaaaaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACt                                               >  1:371577/1‑120 (MQ=255)
                                                                                                                                     aaaaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCaa                >  2:896772/1‑150 (MQ=255)
                                                                                                                                      aaaGCCCTGATTATTCAGGGCTTTGAATATGGCGGAAGCGTAGAGCTTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGc                                         >  2:2218812/1‑124 (MQ=255)
                                                                                                                                      aaaGCCCTGATTATTCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGc                                         <  1:2218812/124‑1 (MQ=255)
                                                                                                                                      aaaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGAc                                                             >  1:472523/1‑104 (MQ=255)
                                                                                                                                      aaaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGAc                                                             <  2:472523/104‑1 (MQ=255)
                                                                                                                                      aaaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAg               >  2:1407513/1‑150 (MQ=255)
                                                                                                                                       aaGCCCTGATTATTCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGAc                                                             <  1:1405453/103‑1 (MQ=255)
                                                                                                                                       aaGCCCTGATTATTCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGAc                                                             >  2:1405453/1‑103 (MQ=255)
                                                                                                                                       aaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTaa                                           <  2:1995892/121‑1 (MQ=255)
                                                                                                                                       aaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTaa                                           >  1:1995892/1‑121 (MQ=255)
                                                                                                                                           ccTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCa          <  1:1173678/150‑1 (MQ=255)
                                                                                                                                            cTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAg         >  1:2076149/1‑150 (MQ=255)
                                                                                                                                                 taatCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGAttt    >  2:212089/1‑150 (MQ=255)
                                                                                                                                                   atCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTcc  <  2:203231/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCAAACCCTGGCAGGTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAAACGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCC  >  NC_002947/4484624‑4484920

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