Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,669,965 C→G 9.6% L148V (CTG→GTG)  metF → 5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,669,9650CG 9.6% 172.6 / 7.0 61L148V (CTG→GTG) metF5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (32/23);  new base G (3/3);  total (35/26)
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.08e-01

CTGGCCGAGTACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGCTACTTCTACTTCGTC  >  NC_002947/5669818‑5670111
                                                                                                                                                   |                                                                                                                                                  
ctgGCCGAGTACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCt                                                                                                                                                   >  3:214252/1‑149 (MQ=255)
   gCCGAGTACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTgg                                                                                                                                                 >  7:126710/1‑148 (MQ=255)
      gAGTACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGt                                                                                                                                             <  7:118203/149‑1 (MQ=255)
         tgcAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGc                                                                                                                                          <  1:212850/147‑1 (MQ=255)
           cAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGTGCATGGGCATGGCCAGTGGTGAACTGCGCTATTCCAGTGACCTGTTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCAGGTGGAAGTgatca                                                                                                                                        <  7:185948/149‑5 (MQ=255)
           cAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCAGGTGGAAGTgatca                                                                                                                                        <  3:178043/149‑5 (MQ=255)
           cAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTg                                                                                                                                        >  4:203707/1‑149 (MQ=255)
            aaGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGc                                                                                                                                       >  2:39183/1‑149 (MQ=255)
                cAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTAc                                                                                                                                   >  6:17749/1‑149 (MQ=255)
                   ccGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCg                                                                                                                                >  3:154392/1‑149 (MQ=255)
                       cATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGaaa                                                                                                                            >  7:105151/1‑149 (MQ=255)
                             aCGCATTGTCGCGCTGCTTGGGTACCTGCCGGC‑GGCATGGGCATGGCCAGTGGTGAACCGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCAccc                                                                                                                     <  4:129441/149‑1 (MQ=255)
                              cGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCAccc                                                                                                                     <  5:17749/149‑1 (MQ=255)
                                     tCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGgc                                                                                                               <  8:126710/148‑1 (MQ=255)
                                       gcgcTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGgcgcg                                                                                                            >  4:213943/1‑149 (MQ=255)
                                           tGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAAc                                                                                                        >  1:293933/1‑149 (MQ=255)
                                               tggtgACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTc                                                                                                     >  7:85342/1‑148 (MQ=255)
                                                ggtgACCTGCCGTCGGGCATGGGCATGGCCTGTGGTGAACTGCGCTATGCAAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTc                                                                                                     >  1:28274/1‑147 (MQ=255)
                                                 gtgACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGa                                                                                                   <  7:189477/148‑1 (MQ=255)
                                                    aCCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCg                                                                                               >  5:208395/1‑149 (MQ=255)
                                                      cTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGa                                                                                             >  5:192933/1‑149 (MQ=255)
                                                          cGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTg                                                                                         <  8:85342/149‑1 (MQ=255)
                                                            tCGGGCATGGGCATGGCCAGTGGTGATCTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGc                                                                                       <  8:261428/149‑1 (MQ=255)
                                                              gggcatgggcatGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCAGGTGGAAGTgatca                                                                                                                                        >  5:12158/1‑94 (MQ=255)
                                                              gggcatgggcatGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCAGGTGGAAGTgatca                                                                                                                                        <  6:12158/98‑5 (MQ=255)
                                                                gcatgggcatggCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAAc                                                                                   >  8:225104/1‑149 (MQ=255)
                                                                gcatgggcatggCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACACGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAAc                                                                                   >  3:213884/1‑149 (MQ=255)
                                                                     ggcatggCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTtgt                                                                              <  3:213943/149‑1 (MQ=255)
                                                                                  gtgAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAg                                                                 >  4:55635/1‑149 (MQ=255)
                                                                                     aaCTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTg                                                                                         >  2:200335/1‑122 (MQ=255)
                                                                                     aaCTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTg                                                                                         <  1:200335/122‑1 (MQ=255)
                                                                                     aaCTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGcc                                                              >  6:45208/1‑149 (MQ=255)
                                                                                      aCTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCg                                                             <  2:290933/149‑1 (MQ=255)
                                                                                          cgcTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGc                                                         <  4:214252/149‑1 (MQ=255)
                                                                                           gcTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGcc                                                        >  7:33708/1‑149 (MQ=255)
                                                                                           gcTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGcc                                                        <  4:112813/149‑1 (MQ=255)
                                                                                                       aCCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCAt                                             >  4:285818/1‑148 (MQ=255)
                                                                                                            gTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCAcc                                         >  1:133665/1‑147 (MQ=255)
                                                                                                            gTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCAcc                                         >  1:300131/1‑147 (MQ=255)
                                                                                                               gAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTa                                    <  8:33708/149‑1 (MQ=255)
                                                                                                               gAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTa                                    >  4:149000/1‑149 (MQ=255)
                                                                                                                  ttCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTActt                                 <  1:39183/149‑1 (MQ=255)
                                                                                                                     aTCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCAcc                                         <  1:28841/138‑1 (MQ=255)
                                                                                                                     aTCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTActtctt                              >  5:272134/1‑149 (MQ=255)
                                                                                                                     aTCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGAGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCAcc                                         >  2:28841/1‑138 (MQ=255)
                                                                                                                         gcgAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCACGGCCGGTGCCGCCATCGCCATCCCCCAGTACTGCTTCAAc                          >  8:140772/1‑149 (MQ=255)
                                                                                                                         gcgAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAAc                          <  6:208395/149‑1 (MQ=255)
                                                                                                                           gaggagACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGt                                     >  1:149415/1‑136 (MQ=255)
                                                                                                                           gaggagACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGt                                     <  2:149415/136‑1 (MQ=255)
                                                                                                                            aggagACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAGATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGcc                       >  1:244742/1‑149 (MQ=255)
                                                                                                                                gaCCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACa                   <  8:105151/149‑1 (MQ=255)
                                                                                                                                  ccGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGc                 >  7:68361/1‑149 (MQ=255)
                                                                                                                                       agccACTTCCAGGTGGAAGTGGCTGCCTAACCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctactt            >  1:22092/4‑149 (MQ=255)
                                                                                                                                        aTCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctacttc           <  2:133665/149‑1 (MQ=255)
                                                                                                                                        aTCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctacttc           <  3:203707/149‑1 (MQ=255)
                                                                                                                                         ccACTTCCAGGTGGAAGTGGCTGCCTACCCTGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAATGCCGGTGCCGACAGCGCCATCAACCAGTACTTCTTCAACGCCGACAGcta               >  5:203870/2‑144 (MQ=255)
                                                                                                                                            cTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCg                                                                                               <  2:160101/61‑1 (MQ=255)
                                                                                                                                            cTTACACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCg                                                                                               >  1:160101/1‑61 (MQ=255)
                                                                                                                                             ttCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGc                 <  7:197011/138‑1 (MQ=255)
                                                                                                                                             ttCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGc                 >  8:197011/1‑138 (MQ=255)
                                                                                                                                                 aCCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCT‑GCCAACTTTGTGCCCAAGGTCAAGGCCGTTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGCTACTTCTACTtcatc  >  8:116852/1‑148 (MQ=255)
                                                                                                                                                   cTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctacttct          >  6:99253/1‑139 (MQ=255)
                                                                                                                                                   cTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctacttct          <  5:99253/139‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTGGCCGAGTACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGCTACTTCTACTTCGTC  >  NC_002947/5669818‑5670111

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