Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,560,900 T→A 7.9% E231D (GAA→GAT PP_4048 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,560,9000TA 7.9% 264.5 / 10.6 87E231D (GAA→GATPP_4048hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/41);  new base A (3/4);  total (42/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 = 1.00e+00

GCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTGGT  >  NC_002947/4560753‑4561040
                                                                                                                                                   |                                                                                                                                            
gCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTc                                                                                                                                            <  1:360857/150‑1 (MQ=255)
gCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGATACGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTc                                                                                                                                            <  1:1071324/150‑1 (MQ=255)
 ccAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTcc                                                                                                                                           >  1:557520/1‑150 (MQ=255)
    gagGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGAAATCCCGc                                                                                                                                        <  2:387590/150‑1 (MQ=255)
       gTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAg                                                                                                                                     <  1:1059513/150‑1 (MQ=255)
       gTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAg                                                                                                                                     <  2:1933617/150‑1 (MQ=255)
        tCACCCGCCCAGCACTCATGGCCATGGGCTCCTGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAgg                                                                                                                                    <  1:1879339/150‑1 (MQ=255)
          aCCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCa                                                                                                                                  >  2:275451/1‑150 (MQ=255)
            ccGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATg                                                                                                                                >  2:1172781/1‑150 (MQ=255)
            ccGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATg                                                                                                                                >  2:1307625/1‑150 (MQ=255)
            ccGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGAAATCCCGCatgc                                                                                                                                    <  1:1024287/146‑4 (MQ=255)
            ccGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGAAATCCCGCatgc                                                                                                                                    >  2:1024287/1‑143 (MQ=255)
                         aTGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCg                                                                                                                   >  2:506436/1‑150 (MQ=255)
                         aTGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCg                                                                                                                   <  2:63220/150‑1 (MQ=255)
                            gCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATACGGGATTTCCCGCAAGGCATGCACCTGCCCGGCGTaa                                                                                                                >  2:740479/1‑149 (MQ=255)
                              cATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGaa                                                                                                              <  2:1809173/150‑1 (MQ=255)
                               aTGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGAATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAAc                                                                                                             >  1:773930/1‑150 (MQ=255)
                                  ggCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTgc                                                                                                          <  1:1176176/150‑1 (MQ=255)
                                   gCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTgcg                                                                                                         >  1:1283986/1‑150 (MQ=255)
                                    cTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTgcgc                                                                                                        >  2:1010775/1‑150 (MQ=255)
                                        ggCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTcc                                                                                                    <  1:422400/150‑1 (MQ=255)
                                           ttttGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATg                                                                                                 >  2:775818/1‑150 (MQ=255)
                                               gcgcATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCt                                                                                             <  1:1867121/150‑1 (MQ=255)
                                                   aTGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCaa                                                                                         >  2:1745370/1‑150 (MQ=255)
                                                       aCCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTcc                                                                                                                                           <  1:1887445/96‑1 (MQ=255)
                                                       aCCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTcc                                                                                                                                           >  2:1887445/1‑96 (MQ=255)
                                                         cgccggcTGGGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCAc                                                                                   >  2:484926/7‑150 (MQ=255)
                                                         cccGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCAc                                                                                   <  1:275451/150‑1 (MQ=255)
                                                               cTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAg                                                                             >  1:1179553/1‑150 (MQ=255)
                                                                         tGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAggac                                                                   <  1:1536915/150‑3 (MQ=255)
                                                                              gggCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTg                                                              <  2:299740/150‑1 (MQ=255)
                                                                                gCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGt                                                            <  2:557520/150‑1 (MQ=255)
                                                                                gCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGAATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGt                                                            <  2:773930/150‑1 (MQ=255)
                                                                                  aGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGGGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATg                                                                                                                                <  2:679946/80‑1 (MQ=255)
                                                                                  aGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGGGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATg                                                                                                                                >  1:679946/1‑80 (MQ=255)
                                                                                  aGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGcc                                                                                              <  1:1033866/114‑1 (MQ=255)
                                                                                  aGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGcc                                                                                              >  2:1033866/1‑114 (MQ=255)
                                                                                     tCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACa                                                                              >  1:1048403/1‑127 (MQ=255)
                                                                                     tCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACa                                                                              <  2:1048403/127‑1 (MQ=255)
                                                                                        cTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACAGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGc                                                    <  2:999019/150‑1 (MQ=255)
                                                                                          ggCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAGCGGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCa                                                  <  1:1810162/150‑1 (MQ=255)
                                                                                            ccTGCTCCAGTACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAGCAGGGGCTGGGCCGCCTGGTTGAACGGCCAGt                                                >  2:1699265/1‑150 (MQ=255)
                                                                                            ccTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGc                                                                                               >  1:114622/1‑103 (MQ=255)
                                                                                            ccTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGc                                                                                               <  2:114622/103‑1 (MQ=255)
                                                                                             cTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGtt                                               >  2:197784/1‑150 (MQ=255)
                                                                                              tGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACg                                                                                                                      <  1:78000/78‑1 (MQ=255)
                                                                                              tGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACg                                                                                                                      >  2:78000/1‑78 (MQ=255)
                                                                                                  ccAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCa                                          <  1:1172781/150‑1 (MQ=255)
                                                                                                    aGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGGCAGTTGCCCAgg                                        <  2:1529146/150‑1 (MQ=255)
                                                                                                     gCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAgg                                        >  1:765553/1‑149 (MQ=255)
                                                                                                     gCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAgg                                        <  2:765553/149‑1 (MQ=255)
                                                                                                        ccTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGc                                    >  2:568426/1‑150 (MQ=255)
                                                                                                        ccTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGc                                    >  2:609370/1‑150 (MQ=255)
                                                                                                         cTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATg                                                                                                                                <  1:149083/57‑1 (MQ=255)
                                                                                                         cTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATg                                                                                                                                >  2:149083/1‑57 (MQ=255)
                                                                                                                  ccTGGCATCGAGCTGCATCAGCGCATGCGGCATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCGTCCATGACCTGCAACCCCACCAACAGTTGCTATGCCGCCTGCTCGAACGGCTAGCTGTCTAGGTAGTTGCGAAGCtt                          >  2:100731/1‑150 (MQ=255)
                                                                                                                   cTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCttt                         >  1:385726/1‑150 (MQ=255)
                                                                                                                   cTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCAAGTTGCCCAGGTAGCTGCGAAGCttt                         <  1:506436/150‑1 (MQ=255)
                                                                                                                      gCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTttgcc                      >  2:58433/1‑150 (MQ=255)
                                                                                                                      gCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTttgcc                      <  1:740479/150‑1 (MQ=255)
                                                                                                                        aTCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTttgccgt                    >  1:324933/1‑150 (MQ=255)
                                                                                                                        aTCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTttgccgt                    >  1:1026555/1‑150 (MQ=255)
                                                                                                                        aTCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTttgccgt                    >  1:1016126/1‑150 (MQ=255)
                                                                                                                              cTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCcgc              <  2:1283986/150‑1 (MQ=255)
                                                                                                                              cTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCcgc              <  1:1638900/150‑1 (MQ=255)
                                                                                                                                gCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCga            <  2:385726/150‑1 (MQ=255)
                                                                                                                                gCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCga            <  1:1699265/150‑1 (MQ=255)
                                                                                                                                gCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCACACCAACAGGTGc                                                                         >  1:1275697/1‑89 (MQ=255)
                                                                                                                                gCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCACACCAACAGGTGc                                                                         <  2:1275697/89‑1 (MQ=255)
                                                                                                                                  aTCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCgaga          >  1:657126/1‑150 (MQ=255)
                                                                                                                                   tCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCa                                          <  2:1808617/117‑1 (MQ=255)
                                                                                                                                   tCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCa                                          >  1:1808617/1‑117 (MQ=255)
                                                                                                                                   tCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGaa         >  1:592136/1‑150 (MQ=255)
                                                                                                                                   tCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGaa         >  2:1738345/1‑150 (MQ=255)
                                                                                                                                    cAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAAc        >  2:644153/1‑150 (MQ=255)
                                                                                                                                      tgccttGCGGGAAATCCCGCATGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAGGCt                           >  2:749429/6‑129 (MQ=255)
                                                                                                                                      tgccttGCGGGAAATCCCGCATGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAGGCt                           <  1:749429/124‑1 (MQ=255)
                                                                                                                                      gcgcATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCt      <  2:1179553/150‑1 (MQ=255)
                                                                                                                                       cgcATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCtt     <  2:1016126/150‑1 (MQ=255)
                                                                                                                                         cttGCGGGAAATCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTttgccgttgc                 >  1:401786/3‑136 (MQ=255)
                                                                                                                                         cttGCGGGAAATCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTttgccgttgc                 <  2:401786/134‑1 (MQ=255)
                                                                                                                                         cATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTgg   <  1:568426/150‑1 (MQ=255)
                                                                                                                                          aTGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTGGt  <  1:1307625/150‑1 (MQ=255)
                                                                                                                                             cGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAAc                                                                                                             <  1:1580308/40‑1 (MQ=255)
                                                                                                                                             cGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAAc                                                                                                             >  2:1580308/1‑40 (MQ=255)
                                                                                                                                                gATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGc                                                    >  1:1123348/1‑94 (MQ=255)
                                                                                                                                                gATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGc                                                    <  2:1123348/94‑1 (MQ=255)
                                                                                                                                                  tttCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGaa                              <  2:535339/114‑1 (MQ=255)
                                                                                                                                                  tttCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGaa                              >  1:535339/1‑114 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
GCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTGGT  >  NC_002947/4560753‑4561040

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: