Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,051,047 C→G 5.1% F360L (TTC→TTG oprB‑III → carbohydrate‑selective porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,051,0470CG 5.1% 257.4 / 3.4 79F360L (TTC→TTGoprB‑IIIcarbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base C (34/41);  new base G (2/2);  total (36/43)
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.59e-01

CGCTTCGCCGTTCACCAAGTGGTATGGCGCCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGAC  >  NC_002947/4050903‑4051194
                                                                                                                                                   |                                                                                                                                                   
cGCTTCGCCGTTCACCAAGTGGTATGGCGCCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAAcg                                                                                                                                                 >  6:6423/1‑149 (MQ=255)
     cGCCGTTCACCAAGTGGTATGGCGCCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCgg                                                                                                                                              <  6:152935/147‑1 (MQ=255)
     cGCCGTTCACCAAGTGGTATGGCGCCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCgg                                                                                                                                              <  4:96252/147‑1 (MQ=255)
      gCCGTTCACCAAGTGGTATGGCGCCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGc                                                                                                                                            <  7:15798/148‑1 (MQ=255)
       ccGTTCACCAAGTGGTATGGCGCCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCca                                                                                                                                          <  3:385184/149‑1 (MQ=255)
         gTTCACCAAGTGGTATGGCGCCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCcag                                                                                                                                         >  5:72045/1‑148 (MQ=255)
         gTTCACCAAGTGGTATGGAGCCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGCGCcag                                                                                                                                         >  7:239807/1‑148 (MQ=255)
             aCCAAGTGGTATGGCGCCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGGCCCCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTGAACGCggcatcttc                                                                                                                                       <  1:348088/148‑8 (MQ=255)
             aCCAAGTGGTATGGCGCCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGtt                                                                                                                                    >  4:409122/1‑149 (MQ=255)
              ccAAGTGGTATGGCGCCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGTCGCGTTCAACGCGGGCCAGCAGTTc                                                                                                                                   >  1:107128/1‑149 (MQ=255)
              ccAAGTGGTATGGCGCCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTc                                                                                                                                   >  8:368456/1‑149 (MQ=255)
               cAAGTGGTATGGCGCCGGCGTGGTGCTCTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCTTGGCACTGGGTTACGG‑‑CCGCACAGAACCCAACCAGCGTACCCTCGATGTGCTGCAAGATGCCGCGTGCAACGACGTCCAGCAtttac                                                                                                                                  >  4:324847/1‑147 (MQ=255)
                 aGTGGTATGGCGCCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCg                                                                                                                                <  2:257627/149‑1 (MQ=255)
                 aGTGGTATGGCGCCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCg                                                                                                                                <  5:6423/149‑1 (MQ=255)
                     gTATGGCGCCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCTATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACa                                                                                                                            <  6:418585/149‑1 (MQ=255)
                        tGGCGCCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACACAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATc                                                                                                                          >  7:276943/1‑148 (MQ=255)
                         ggCGCCGGCGTGGAGCTGGACCAGCCGGTCGACGGCCGCCGGCACGCTACCGTTGCACGGGGGGACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGAGGTGCGGGAAGAAGGCGCGTTCCACGCGGGCCAGCAGTTCCCGGAGATCGa                                                                                                                        <  8:207903/149‑1 (MQ=255)
                          gtgcCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAAc                                                                                                                                                  <  7:188408/120‑1 (MQ=255)
                          gtgcCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAAc                                                                                                                                                  >  8:188408/3‑122 (MQ=255)
                          gcgcCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGa                                                                                                                        >  4:419971/1‑148 (MQ=255)
                            gcCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAg                                                                                                                     <  3:75139/149‑1 (MQ=255)
                             ccGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAgc                                                                                                                    >  2:321342/1‑149 (MQ=255)
                              cGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAgcg                                                                                                                   <  6:296778/149‑1 (MQ=255)
                              cGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAgcg                                                                                                                   <  3:409122/149‑1 (MQ=255)
                                    ggtgCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGaa                                                                                                              >  8:64531/1‑148 (MQ=255)
                                        cggTACAAGCCGTTCGAAGGCCGCGAGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGGGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGTGCCGAACAGt                                                                                                          <  8:239807/146‑1 (MQ=255)
                                        cTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGtt                                                                                                         >  3:81772/1‑149 (MQ=255)
                                        cTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGtt                                                                                                         >  8:227763/1‑149 (MQ=255)
                                        cTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGt                                                                                                          <  7:396156/148‑1 (MQ=255)
                                             cAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATc                                                                                                     <  3:172878/148‑1 (MQ=255)
                                             cAAGCAGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCg                                                                                                    <  3:324847/149‑1 (MQ=255)
                                               aGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGaa                                                                                                  <  3:419971/149‑1 (MQ=255)
                                                gCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGaa                                                                                                  >  5:121206/1‑148 (MQ=255)
                                                  cGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCCCACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACt                                                                                                <  1:331529/149‑1 (MQ=255)
                                                   gTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCg                                                                                                    >  5:125887/1‑143 (MQ=255)
                                                   gTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCg                                                                                                    <  6:125887/143‑1 (MQ=255)
                                                    ttCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAAc                                                                                                             >  7:72104/1‑133 (MQ=255)
                                                    ttCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAAc                                                                                                             <  8:72104/133‑1 (MQ=255)
                                                      cGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGc                                                                                            >  7:116097/1‑148 (MQ=255)
                                                      cGAAGGCCGCCCGAAGGATACCGTGGCACTGAGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGc                                                                                            >  1:265092/1‑148 (MQ=255)
                                                          ggCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACgg                                                                                       <  8:68790/149‑1 (MQ=255)
                                                             cgccCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTa                                                                                    <  1:321342/149‑1 (MQ=255)
                                                                ccGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTAc                                                                                   >  2:429300/1‑147 (MQ=255)
                                                                ccGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCa                                                                                 <  4:101736/149‑1 (MQ=255)
                                                                   aaGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGc                                                                              <  8:276943/149‑1 (MQ=255)
                                                                             tGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGTCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATgg                                                                    <  2:107128/149‑1 (MQ=255)
                                                                             tGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATgg                                                                    <  6:443256/149‑1 (MQ=255)
                                                                             tGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATgg                                                                    <  2:265092/149‑1 (MQ=255)
                                                                                       ttACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTgc                                                          <  6:72045/149‑1 (MQ=255)
                                                                                          cGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTgcgc                                                        <  4:81772/148‑1 (MQ=255)
                                                                                          cGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGcc                                                       >  3:380016/1‑149 (MQ=255)
                                                                                          cGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATCGCTGAACCTGCGcc                                                       <  7:368456/149‑1 (MQ=255)
                                                                                               ccGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCCGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCgg                                                    <  7:227763/149‑1 (MQ=255)
                                                                                                cgcgCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGa                                                   >  5:291109/1‑149 (MQ=255)
                                                                                                 gcgcCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGAc                                                  <  1:429300/149‑1 (MQ=255)
                                                                                                     cGTACCCAACCCTCTTATCCGCGATGTGCTGTAAGATTCCGCGTTCAAGGCGGTCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACATGCGCCCGGACGTGc                                              <  2:318547/149‑1 (MQ=255)
                                                                                                         cccAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATa                                          >  1:21290/1‑149 (MQ=255)
                                                                                                         ccaaaCACGCGAAGCCGCGAAGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACAACGACAGCACCGAAAACTTGATCCAACTCAGCTACGCCTACCAGGCCACGCCAAGGCAGAACCTCCGCACGGAGGTGCAATa                                          >  8:118634/4‑149 (MQ=255)
                                                                                                          ccAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACGAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATAc                                         <  6:121206/149‑1 (MQ=255)
                                                                                                          ccAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATa                                          >  5:444756/1‑148 (MQ=255)
                                                                                                          ccAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATa                                          >  4:94699/1‑148 (MQ=255)
                                                                                                               ccGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATAcatcat                                    >  8:381319/1‑149 (MQ=255)
                                                                                                                 gcgTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGt                                                <  7:124155/135‑1 (MQ=255)
                                                                                                                 gcgTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGt                                                >  8:124155/1‑135 (MQ=255)
                                                                                                                 gcgTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACatcatc                                   >  5:390727/1‑148 (MQ=255)
                                                                                                                             tgtgCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTTCGCCCGGACGTGCAATACATCATCg                                  <  2:289759/137‑1 (MQ=255)
                                                                                                                             tgtgCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTACGCCCGGACGTGCAATACATCATCg                                  >  1:289759/1‑137 (MQ=255)
                                                                                                                                 cTGGAAGATGCCGCGTTCTACGCGGGCCAGCTGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGATCTGAGCTACGGCTACCAGGCCACGCCTTGGCTGAAGCTGCTCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTtctc                  <  7:118634/149‑1 (MQ=255)
                                                                                                                                    gAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCgg                <  7:33139/148‑1 (MQ=255)
                                                                                                                                    gAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCg                 <  8:38708/147‑1 (MQ=255)
                                                                                                                                       tggcccGCGTTGAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGa                                 <  4:120394/124‑1 (MQ=255)
                                                                                                                                       tggcccGCGTTGAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGa                                 >  3:120394/5‑128 (MQ=255)
                                                                                                                                        aTGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGc                        >  4:26156/1‑136 (MQ=255)
                                                                                                                                        aTGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGc                        <  3:26156/136‑1 (MQ=255)
                                                                                                                                          gCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGa         >  7:375708/1‑149 (MQ=255)
                                                                                                                                           ccGCGTTGAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGa         >  6:295161/1‑148 (MQ=255)
                                                                                                                                             gcgTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACAt      >  6:347637/1‑149 (MQ=255)
                                                                                                                                               gTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACatcatc                                   <  4:396246/118‑1 (MQ=255)
                                                                                                                                               gTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACatcatc                                   >  3:396246/1‑118 (MQ=255)
                                                                                                                                                   cAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGAc  >  2:23761/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGCTTCGCCGTTCACCAAGTGGTATGGCGCCGGCGTGGTGCTGTACAAGCCGTTCGAAGGCCGCCCGAAGGATACCGTGGCACTGGGTTACGG‑‑CCGCGCCGTACCCAACCCGCGTAGCCGCGATGTGCTGGAAGATGCCGCGTTCAACGCGGGCCAGCAGTTCCCGGACATCGACAGCGCCGAACAGTTGATCGAACTGAGCTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGAC  >  NC_002947/4050903‑4051194

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