Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,256,133 G→A 5.9% L221L (CTG→CTA folM → bifunctional dihydrofolate reductase/dihydromonapterin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,256,1330GA 5.9% 207.3 / 2.2 68L221L (CTG→CTAfolMbifunctional dihydrofolate reductase/dihydromonapterin reductase
Reads supporting (aligned to +/- strand):  ref base G (30/34);  new base A (2/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

CCCAAGGTGCGTGTGAACGGTATCCTGCCCGGGCTGCTGATCCTCAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAGCAACAGCTCGAGCTCGAGGCCGCCGCGTTCCGGCGCCTGGTCGAACACCT  >  NC_002947/5255987‑5256279
                                                                                                                                                  |                                                                                                                                                  
cccAAGGTGCGTGTGAACGGTATCCTGCCCGGGCTGCTGATCCTCAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCt                                                                                                                                                  >  2:135609/1‑149 (MQ=255)
       tGCGTGTGAACGGTATCCTGCCCGGGCTGCTGATCCTCAAGGAAGGGGGGGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCgccagc                                                                                                                                           <  4:300266/149‑1 (MQ=255)
          gtgtGAACGGTATCCTGCCCGGGCTGCTGATCCTCAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGccagccag                                                                                                                                        <  4:282731/149‑1 (MQ=255)
            gtgAACGGTATCCTGCCCGGGCTGCTGATCCTCAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTa                                                                                                                                      >  1:131051/1‑149 (MQ=255)
             tgAACGGTATGCTGCCCGGGCTGCTGATCCTCCAGGCAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAACAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTAc                                                                                                                                     <  4:89254/149‑1 (MQ=255)
               aaCGGTATCCTGCCCGGGCTGCTGATCCTCAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACa                                                                                                                                    <  3:109705/148‑1 (MQ=255)
                         tGCCCGGGCTGCTGATCCTCAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGc                                                                                                                         >  2:152594/1‑149 (MQ=255)
                           cccGGGCTGCTGATCCTCAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGc                                                                                                                         <  1:152594/147‑1 (MQ=255)
                             cGGGCTGCTGATCCTCAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAgg                                                                                                                     <  5:2950/149‑1 (MQ=255)
                                   gctgATCCTCAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGtca                                                                                                               >  4:149475/1‑149 (MQ=255)
                                    ctgATCCTCAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCTCCAGCCAGTACAGCACCGGCAGCCAGGTGGtc                                                                                                                <  1:201651/147‑1 (MQ=255)
                                    ctgATCCTCAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGtcat                                                                                                              >  6:89140/1‑149 (MQ=255)
                                    ctgATCCTCAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTACTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGtca                                                                                                               >  7:69604/1‑148 (MQ=255)
                                    ctgATCCTCAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTACTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGtca                                                                                                               >  5:39090/1‑148 (MQ=255)
                                       aTCCTCAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCgcca                                                                                                                                             >  6:111263/1‑115 (MQ=255)
                                       aTCCTCAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCgcca                                                                                                                                             <  5:111263/115‑1 (MQ=255)
                                          cTCAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCgc                                                                                                                                               <  2:148508/110‑1 (MQ=255)
                                          cTCAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCgc                                                                                                                                               >  1:148508/1‑110 (MQ=255)
                                           tCAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGt                                                                                                       >  3:158613/1‑149 (MQ=255)
                                            cAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCCGCCAGGTGGTCCTCACCGGTg                                                                                                      >  6:154340/1‑149 (MQ=255)
                                             aaggaaggGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGTAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACggcg                                                                                                      <  8:178159/148‑3 (MQ=255)
                                             aaggaaggGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTCGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTgg                                                                                                     >  1:45744/1‑149 (MQ=255)
                                                 aaggGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGgccgc                                                                                                 <  5:89140/149‑1 (MQ=255)
                                                     gggTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCAcc                                                                                             >  1:89780/1‑149 (MQ=255)
                                                     gggTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCAcc                                                                                             >  6:147061/1‑149 (MQ=255)
                                                      ggTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCt                                                                                            >  1:192481/1‑149 (MQ=255)
                                                            gccgCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTgaagaa                                                                                      >  7:208376/3‑149 (MQ=255)
                                                            gacgCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTgaagaa                                                                                      <  1:135609/149‑1 (MQ=255)
                                                                ccTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGc                                                                                  <  2:89780/149‑1 (MQ=255)
                                                                    ccGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGa                                                                              <  3:149475/149‑1 (MQ=255)
                                                                      gCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACt                                                                            >  2:16371/1‑149 (MQ=255)
                                                                      gCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACt                                                                            >  6:156655/1‑149 (MQ=255)
                                                                       ccAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACtt                                                                           >  7:285001/1‑149 (MQ=255)
                                                                          gCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTgag                                                                        <  6:38637/149‑1 (MQ=255)
                                                                          gCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTgag                                                                        <  1:16371/149‑1 (MQ=255)
                                                                            aaaCCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTga                                                                         >  4:307447/1‑146 (MQ=255)
                                                                              aCCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGagaa                                                                    <  4:158613/149‑1 (MQ=255)
                                                                              aCCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGagaa                                                                    <  2:131051/149‑1 (MQ=255)
                                                                                ccTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAAc                                                                                     <  1:260218/130‑1 (MQ=255)
                                                                                ccTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAAc                                                                                     >  2:260218/1‑130 (MQ=255)
                                                                                    aaaaaGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGAcc                                                                <  3:307447/147‑1 (MQ=255)
                                                                                     aaaaGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCAt                                                              <  8:285001/148‑1 (MQ=255)
                                                                                       aaGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGAc                                                           <  8:208376/149‑1 (MQ=255)
                                                                                           ccTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCca                                                       <  5:147061/149‑1 (MQ=255)
                                                                                           ccTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCca                                                       >  6:214547/1‑149 (MQ=255)
                                                                                           ccTCGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCca                                                       <  2:45744/149‑1 (MQ=255)
                                                                                                cTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAGc                                                   >  8:180116/1‑148 (MQ=255)
                                                                                                    aGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAGCAACag                                              <  5:210794/149‑1 (MQ=255)
                                                                                                        cGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCg                                                                                                                                                >  1:152043/1‑47 (MQ=255)
                                                                                                        cGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCg                                                                                                                                                <  2:152043/47‑1 (MQ=255)
                                                                                                        cGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTgaa                                                                                         <  8:189415/102‑1 (MQ=255)
                                                                                                         gAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTGATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAGCAACagctc                                           >  1:101384/1‑147 (MQ=255)
                                                                                                         gAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAg                                                    >  2:252311/1‑138 (MQ=255)
                                                                                                         gAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAg                                                    <  1:252311/138‑1 (MQ=255)
                                                                                                         gAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAGCAACagctcg                                          <  1:253080/148‑1 (MQ=255)
                                                                                                              cGGCGCCGGCCCGCTGATGGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACACCACTGGTAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAGCAACagctcgagctcg                                    <  1:264741/149‑1 (MQ=255)
                                                                                                                 cgcCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAGCAACAGCTCGAGCTCGAgg                                 <  7:180116/149‑1 (MQ=255)
                                                                                                                          gCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAGCAACAGCTCGAGCTCGAGGCCGCCGCGt                        <  2:192481/149‑1 (MQ=255)
                                                                                                                              aTCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAGCAACAGCTCGAGCTCGAGGCCGCCGCGTTCCg                    <  5:214547/149‑1 (MQ=255)
                                                                                                                                  agaCCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAGCAACAGCTCGAGCTCGAGGCCGCCGCGTTCCGgcgc                >  8:293414/1‑149 (MQ=255)
                                                                                                                                      cGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAGCAACAGCTCGAGCTCGAGGCCGCCGCGTTCCGGCGCCTgg            >  4:84134/1‑149 (MQ=255)
                                                                                                                                       gTCAAGTACCTACTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAGCAACAGCTCGAGCTCGAGGCCGCCGCGTTCCGGCGCCTGGt           <  8:69604/149‑1 (MQ=255)
                                                                                                                                       gTCAAGTACCTACTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAGCAACAGCTCGAGCTCGAGGCCGCCGCGTTCCGGCGCCTGGt           <  6:39090/149‑1 (MQ=255)
                                                                                                                                           aGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAGCAACAGCTCGAGCTCGAGGCCGCCGCGTTCCGGCGCCTGGTCGaa       >  3:129160/1‑149 (MQ=255)
                                                                                                                                               ccTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAGCAACAGCTCGAGCTCGAGGCCGCCGCGTTCCGGCGCCTGGTCGAACAcc   >  4:68861/1‑149 (MQ=255)
                                                                                                                                                cTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAGCAACAGCTCGAGCTCGAGGCCGCCGCGTTCCgg                   >  3:38002/1‑132 (MQ=255)
                                                                                                                                                cTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAGCAACAGCTCGAGCTCGAGGCCGCCGCGTTCCgg                   <  4:38002/132‑1 (MQ=255)
                                                                                                                                                cTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCAACAGCAACAGCTCGAGCTCGAGGCCGCCGCGTTCCGGCGCCTGGTCGAACACCt  >  2:173242/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CCCAAGGTGCGTGTGAACGGTATCCTGCCCGGGCTGCTGATCCTCAAGGAAGGGGGTGACGACGCCTACCGCCAGCAAACCCTGAAAAAGGCCTTGCTGGAGTTCGAACCCGGCGCCGGCCCGCTGATCGAGACCGTCAAGTACCTGCTCGCCAGCCAGTACAGCACCGGCAGCCAGGTGGTCATCAACGGTGGCCGCCACCTGAAGAACCGCATGACTTGAGAGAAGACCATGACCCCACAGCAACAGCTCGAGCTCGAGGCCGCCGCGTTCCGGCGCCTGGTCGAACACCT  >  NC_002947/5255987‑5256279

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