Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,812,662 T→A 6.7% I329N (ATT→AAT)  PP_3368 → MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,812,6620TA 6.7% 180.6 / 2.8 60I329N (ATT→AAT) PP_3368MFS transporter
Reads supporting (aligned to +/- strand):  ref base T (27/29);  new base A (2/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

GCTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGCGCCGCGACTG  >  NC_002947/3812521‑3812805
                                                                                                                                             |                                                                                                                                               
gCTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCAt                                                                                                                                          <  7:32744/149‑1 (MQ=255)
    ggaaCTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAgg                                                                                                                                      <  6:244535/146‑1 (MQ=255)
     gTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGGCAATGTCCATGAgg                                                                                                                                      <  1:259763/148‑1 (MQ=255)
           aTGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAg                                                                                                                               <  5:351538/149‑1 (MQ=255)
                    gtcgCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAgc                                                                                                                      >  6:239663/1‑149 (MQ=255)
                     tcgCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAgcg                                                                                                                     <  7:350075/149‑1 (MQ=255)
                        ccGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGACGCCTAGCCTGGCAgcggcg                                                                                                                  >  7:19413/1‑149 (MQ=255)
                           tgctgcTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTcc                                                                                                               <  3:71673/149‑1 (MQ=255)
                              tgctgTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAAc                                                                                                            <  7:308126/149‑1 (MQ=255)
                               gctgTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAAc                                                                                                            >  6:26927/1‑148 (MQ=255)
                                ctgTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAAc                                                                                                            <  6:160887/147‑1 (MQ=255)
                                 tgTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCAt                                                                                                         >  6:391235/1‑149 (MQ=255)
                                     cTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATgcg                                                                                                      <  3:323216/148‑1 (MQ=255)
                                      tttACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATgcggc                                                                                                    >  7:130391/1‑149 (MQ=255)
                                        tACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATgcggc                                                                                                    <  3:329134/147‑1 (MQ=255)
                                           cTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCCGCGGCGTCCAACCATGCGGCGTTc                                                                                                >  3:376321/1‑148 (MQ=255)
                                            ttCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCaa                                                                                              >  3:318743/1‑149 (MQ=255)
                                               gccgccCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAAcc                                                                                            <  2:131253/148‑1 (MQ=255)
                                                ccgccCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTg                                                                                          >  6:81275/1‑149 (MQ=255)
                                                   cccAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGcc                                                                                       <  5:239663/149‑1 (MQ=255)
                                                    ccAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTTGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGcc                                                                                       >  8:284231/1‑148 (MQ=255)
                                                     cAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCaa                                                                                     >  4:238415/1‑149 (MQ=255)
                                                      aGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAAt                                                                                    <  8:130391/149‑1 (MQ=255)
                                                              tgGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGGCTTGGAGGCAGCGCTACAGATGCGCCTGATGGACTTTGCCCATGAGGCGCCTTGCCTGGCAGCGGGGTCCAACCATGCGGCGTTCAAGCTGGCCAATGCCCttg                                                                             <  6:314538/148‑3 (MQ=255)
                                                              tgGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCCCCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCATTGCCCTggg                                                                            <  5:391235/149‑1 (MQ=255)
                                                              tgGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTggg                                                                            <  5:26927/149‑1 (MQ=255)
                                                              tgGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTg                                                                              <  6:101799/147‑1 (MQ=255)
                                                                gAGCNTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGc                                                                          >  5:109889/1‑149 (MQ=255)
                                                                      gctgctgGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGc                                                                    >  5:262551/1‑149 (MQ=255)
                                                                       ctgctgGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCt                                                                   <  4:262577/149‑1 (MQ=255)
                                                                        tgctgGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTggg                                                                            <  1:151236/139‑1 (MQ=255)
                                                                        tgctgGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTggg                                                                            >  2:151236/1‑139 (MQ=255)
                                                                              gCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGc                                                            <  2:214237/149‑1 (MQ=255)
                                                                                         gtggGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCTGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCAc                                                 >  1:123110/1‑149 (MQ=255)
                                                                                           gggCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAg                                                              >  3:218792/1‑134 (MQ=255)
                                                                                           gggCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAg                                                              <  4:218792/134‑1 (MQ=255)
                                                                                              cACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCg                                            >  1:151431/1‑149 (MQ=255)
                                                                                               aCCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGGCTAGCCTGTCAGCGGCGTCCAACCTTGCGGCGTTACAACTGGCCAATCGCCTGCGGCGGTGGCTGGGAGGCATGGCGATCACTGCCgg                                           >  5:327141/1‑149 (MQ=255)
                                                                                               aCCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGGCGTGGCTGGGAGGCAAGGCGATCACTGCCgg                                           >  5:175397/1‑149 (MQ=255)
                                                                                               aCCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCgg                                           >  5:171228/1‑149 (MQ=255)
                                                                                               aCCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCgg                                           >  4:222790/1‑149 (MQ=255)
                                                                                                ccATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCgg                                           >  2:200/1‑148 (MQ=255)
                                                                                                 cATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCa                                         >  7:293893/1‑149 (MQ=255)
                                                                                                     aTCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATgg                                      <  1:346585/148‑1 (MQ=255)
                                                                                                             ggcggcACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGAc                               <  2:151431/147‑1 (MQ=255)
                                                                                                              gcggcACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCa                             <  3:3930/148‑1 (MQ=255)
                                                                                                                ggcACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCa                          >  7:69504/1‑149 (MQ=255)
                                                                                                                ggcACAGCTACAGCTCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGAGGCGTACAACCTGGCCAAAGAACTGGGGCCGTCGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCa                          >  8:162511/1‑149 (MQ=255)
                                                                                                                     cGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGAcc                              >  4:135361/1‑140 (MQ=255)
                                                                                                                     cGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGAcc                              <  3:135361/140‑1 (MQ=255)
                                                                                                                     cGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGa                     <  6:109889/149‑1 (MQ=255)
                                                                                                                     cGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGa                     <  5:81275/149‑1 (MQ=255)
                                                                                                                      gCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGAt                    <  4:376321/149‑1 (MQ=255)
                                                                                                                      gCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGAt                    <  4:318743/149‑1 (MQ=255)
                                                                                                                            gATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCg              <  3:138315/149‑1 (MQ=255)
                                                                                                                                 gCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGcgcc         >  5:181767/1‑149 (MQ=255)
                                                                                                                                  ccTCATGGGCAATGTCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGcgccg        >  4:307045/1‑149 (MQ=255)
                                                                                                                                   cTCATGGGCAATGTCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGcgccgc       <  3:307045/149‑1 (MQ=255)
                                                                                                                                    tCATGGGCAATGTCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGcgccgc       >  3:9970/1‑148 (MQ=255)
                                                                                                                                        ggACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGCGCCGCGACTg  >  4:37456/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                               
GCTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGCGCCGCGACTG  >  NC_002947/3812521‑3812805

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: