Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,758,946 G→A 6.8% D71N (GAC→AAC)  PP_3323 → CobW/P47K family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,758,9460GA 6.8% 170.0 / 2.4 60D71N (GAC→AAC) PP_3323CobW/P47K family protein
Reads supporting (aligned to +/- strand):  ref base G (29/26);  new base A (2/2);  total (31/29)
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.82e-01

TCGGAGCCGGCAAGACTACCTTGCTCAACCATATTCTGAACAACCGCGAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGACCAATTGGTCATCGAATCGACGGGCATCTCCGAGCCACTGCCGGTGGC  >  NC_002947/3758803‑3759094
                                                                                                                                               |                                                                                                                                                    
tCGGAGCCGGCAAGACTACCTTGCTCAACCATATTCTGAACAACCGCGAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACcgacga                                                                                                                                                  >  4:278800/1‑148 (MQ=255)
  ggAGCCGGCAAGTCTACCTTGCTCATCCATATTCTTAACAACCGCGTAGGCCGGCGAGTCGCGGTGATCGTGAACGACATGAGCGATGTGAACATCGATGCCTCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGaaaa                                                                                                                                               <  4:86263/149‑1 (MQ=255)
     gCCGGCAAGACTACCTTGCTCAACCATATTCTGAACAACCGCGAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCt                                                                                                                                            <  3:100453/149‑1 (MQ=255)
      ccGGCAAGACTACCTTGCTCAACCATATTCTGAACAACCGCGAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCtg                                                                                                                                           <  1:360327/149‑1 (MQ=255)
        ggCAAGACTACCTTGCTCAACCATATTCTGAACAACCGCGAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCtggt                                                                                                                                         <  7:37097/149‑1 (MQ=255)
          cAAGACTACCTTGCTCAACCATATTCTGAACAACCGCGAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCtggtgg                                                                                                                                       <  7:292761/149‑1 (MQ=255)
             gACTACCTTGCTCAACCATATTCTGAACAACCGCGAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCtggtg                                                                                                                                        >  5:288407/1‑145 (MQ=255)
                    ttGCTCAACCATATTCTGAACAACCGCGAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCaa                                                                                                                             >  1:407283/1‑149 (MQ=255)
                        tCAACCATATTCTGAACAACCGCGAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGc                                                                                                                         >  4:115598/1‑149 (MQ=255)
                          aaCCATATTCTGAACAACCGCGAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATAGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTg                                                                                                                       >  2:375722/1‑149 (MQ=255)
                            ccATATTCTGAACAACCGCGAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCa                                                                                                                     >  1:303420/1‑149 (MQ=255)
                                 ttCTGAACAACCGCGAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTtgc                                                                                                                >  5:8442/1‑149 (MQ=255)
                                   cTGAACAACCGCGAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTtgct                                                                                                               <  3:278800/148‑1 (MQ=255)
                                      aacaacCGCGAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAg                                                                                                                                <  6:199379/128‑1 (MQ=255)
                                      aacaacCGCGAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAg                                                                                                                                >  5:199379/1‑128 (MQ=255)
                                      aacaacCGCGAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCAc                                                                                                           >  3:381411/1‑149 (MQ=255)
                                      aacaacCGCGAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCAc                                                                                                           <  6:303632/149‑1 (MQ=255)
                                           ccGCGAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCCATGCCGCACTGATACGCACGGGTGGCGCGGATCTGTCACGCACCGACGAAAAGCTTGTGGAGAGGAGCAATGGCTGCATTTGAGGGACCTTgc                                                                                                      >  8:289699/1‑149 (MQ=255)
                                            cgcgAAGGCCGGCGAGTCGCGGTGATCGTCTACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTgcg                                                                                                     >  1:33436/1‑149 (MQ=255)
                                            cgcgAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTtg                                                                                                                 <  8:273673/137‑1 (MQ=255)
                                            cgcgAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTtg                                                                                                                 >  7:273673/1‑137 (MQ=255)
                                             gcgAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTgcgc                                                                                                    >  7:189356/1‑149 (MQ=255)
                                                aaGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTgcgact                                                                                                  >  5:225490/1‑145 (MQ=255)
                                                aaGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGaa                                                                                                 <  1:196387/149‑1 (MQ=255)
                                                     cGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCt                                                                                            >  8:208020/1‑149 (MQ=255)
                                                       gCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATACCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTc                                                                                          >  1:294102/1‑149 (MQ=255)
                                                        cGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTc                                                                                          <  4:142074/148‑1 (MQ=255)
                                                        cGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTc                                                                                          >  3:142074/1‑148 (MQ=255)
                                                                 gTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGt                                                                                <  6:8442/149‑1 (MQ=255)
                                                                    aTCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGa                                                                             >  6:216603/1‑149 (MQ=255)
                                                                     tCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGAc                                                                            <  2:303420/149‑1 (MQ=255)
                                                                       gTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCg                                                                          >  1:244106/1‑149 (MQ=255)
                                                                         cAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTc                                                                        >  1:90469/1‑149 (MQ=255)
                                                                          aaCGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTc                                                                        >  4:231016/1‑148 (MQ=255)
                                                                           aCGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTa                                                                      <  1:375722/149‑1 (MQ=255)
                                                                           aCGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCCCGAAGTGGACCGTCTa                                                                      <  2:407283/149‑1 (MQ=255)
                                                                               cATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTc                                                                                          >  8:77346/1‑125 (MQ=255)
                                                                               cATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTc                                                                                          <  7:77346/125‑1 (MQ=255)
                                                                                     cGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAgg                                                            <  6:288407/149‑1 (MQ=255)
                                                                                       aTGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAaga                                                              >  5:34575/1‑145 (MQ=255)
                                                                                       aTGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAaga                                                              <  6:34575/145‑1 (MQ=255)
                                                                                              cATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTc                                                    >  4:396102/1‑148 (MQ=255)
                                                                                               aTCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGa                                                  <  8:189356/149‑1 (MQ=255)
                                                                                                tCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGAc                                                 <  7:208020/149‑1 (MQ=255)
                                                                                                tCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGATGGTGGAGATGAGCAATGGCTGCATTTGCTTCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCTTCTAGCCAAAGAGGGTCGTTTCGAc                                                 <  7:289699/149‑1 (MQ=255)
                                                                                                      ccGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGACCAAtt                                            <  3:396102/148‑1 (MQ=255)
                                                                                                           cTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGACCAATTGGTCat                                      <  2:244106/149‑1 (MQ=255)
                                                                                                              aTACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAg                                                                     <  6:173636/115‑1 (MQ=255)
                                                                                                              aTACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAg                                                                     >  5:173636/1‑115 (MQ=255)
                                                                                                              aTACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGACCAATTGGTCatcg                                    >  2:56539/1‑148 (MQ=255)
                                                                                                               accGCAAGGTGGCCGCCGAGCTGGCACCCACCCACTAAAAGCTGGTGGAGATGACCAATGCCTTCATTTGATGAACCTTGCGCGAACACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGACCAATTGGTCatcgaa                                  <  1:105952/147‑1 (MQ=255)
                                                                                                                 cGCAAGGGTGGCGCCGAGCTGTCACGCACCAACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGACCAATTGGTCatcgaat                                 <  5:397328/148‑1 (MQ=255)
                                                                                                                 cGCAAGGGTGGCGCCGAGCTGTCACGCACCAACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGACCAATTGGTCatcgaat                                 >  6:397328/1‑148 (MQ=255)
                                                                                                                 cGCAAGGGTGGCGCCGAGCTGTCACGCACCAACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGACCAATTGGTCatcgaat                                 >  6:391598/1‑148 (MQ=255)
                                                                                                                 cGCAAGGGTGGCGCCGAGCTGTCACGCACCAACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGACCAATTGGTCatcgaat                                 <  5:391598/148‑1 (MQ=255)
                                                                                                                  gCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGACCAATTGGTCatcgaatcg                               <  4:381411/149‑1 (MQ=255)
                                                                                                                    aaGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGACCAATTGGTCATCGAATCGAc                             >  3:403855/1‑149 (MQ=255)
                                                                                                                    aaGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGACCAATTGGTCATCGAATCGAc                             >  8:1462/1‑149 (MQ=255)
                                                                                                                         tGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGACCAATTGGTCATCGAATCGACGGGCa                        >  5:164094/1‑149 (MQ=255)
                                                                                                                           gcgcCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGACCAATTGGTCATCGAATCGACGGGCAtc                      <  5:216603/149‑1 (MQ=255)
                                                                                                                                      tCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGACCAATTGGTCATCGAATCGACGGGCATCTCCGAGCCACt           >  4:297449/1‑149 (MQ=255)
                                                                                                                                              cgacgaAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGACCAATTGGTCATCGAATCGACGGGCATCTCCGAGCCACTGCCGGTgg   <  1:56539/149‑1 (MQ=255)
                                                                                                                                               gacgaAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGACCAATTGGTCATCGAATCGACGGGCATCTCCGAGCCACTGCCGGTGGc  >  2:333779/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TCGGAGCCGGCAAGACTACCTTGCTCAACCATATTCTGAACAACCGCGAAGGCCGGCGAGTCGCGGTGATCGTCAACGACATGAGCGATGTGAACATCGATGCCGCACTGATACGCAAGGGTGGCGCCGAGCTGTCACGCACCGACGAAAAGCTGGTGGAGATGAGCAATGGCTGCATTTGCTGCACCTTGCGCGAAGACCTTCTGCTCGAAGTGGACCGTCTAGCCAAAGAGGGTCGTTTCGACCAATTGGTCATCGAATCGACGGGCATCTCCGAGCCACTGCCGGTGGC  >  NC_002947/3758803‑3759094

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: ‑