Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,301,533 G→A 7.0% A95T (GCT→ACT)  PP_1136 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,301,5330GA 7.0% 164.7 / 2.3 58A95T (GCT→ACT) PP_1136hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/28);  new base A (2/2);  total (27/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 = 6.83e-01

AACGATTATGGTAGGGCGTGTAAAGCATGATGGTACAGCAGATTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGTTGAAAAGATGAGTGCAAGTTCAGCGGGGCGTGGAAGAAACAATGACG  >  NC_002947/1301385‑1301680
                                                                                                                                                    |                                                                                                                                                    
aaCGATTATGGTAGGGCGTGTAAAGCATGATGGTACAGCAGATTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCg                                                                                                                                                      >  2:188263/1‑149 (MQ=255)
  cttttATGGTAGGGCTTGTAAAGCATGATTGTACAGCAGATTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTTATTCCTCAACCTGACGGAGGCGTAATCGCt                                                                                                                                                    <  8:60706/146‑1 (MQ=255)
       aTGGTAGGGCGTGTAAAGCATGATGGTACAGCAGATTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGtt                                                                                                                                                <  3:165301/148‑1 (MQ=255)
        tGGTAGGGCGTGTAAAGCATGATGGTACAGCAGATTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGttt                                                                                                                                               <  2:127140/148‑1 (MQ=255)
        tGGTAGGGCGTGTAAAGCATGATGGTACAGCAGATTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTACGCGCGCTGAAGGGTTCATACCTAAACCTGAGGGAGGCGTAAACGCTGGTTTg                                                                                                                                              >  5:205158/1‑149 (MQ=255)
         ggTAGGGCGTGTAAAGCATGATGGTACAGCAGATTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTg                                                                                                                                              >  4:82589/1‑148 (MQ=255)
          gTAGGGCGTGTAAAGCATGATGGTACAGCAGATTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGta                                                                                                                                            >  4:76044/1‑149 (MQ=255)
            aGGGCGTGTAAAGCATGATGGTACAGCAGATTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGtat                                                                                                                                           >  2:153775/1‑148 (MQ=255)
             gggCGTGTAAAGCATGATGGTACAGCAGATTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATaa                                                                                                                                         <  3:82589/149‑1 (MQ=255)
                   gtAAAGCATGATGGTACAGCAGATTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGTAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATc                                                                                                                                   >  6:74028/1‑149 (MQ=255)
                     aatgCAGGATGGTACAGCAGATTTGACGTTTGGTTTGCGAGGTTTCAGCTCCAATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCg                                                                                                                                  <  8:289973/145‑1 (MQ=255)
                      aaGCATGATGGTACAGCAGATTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCg                                                                                                                                  >  4:70327/1‑147 (MQ=255)
                               ggTACAGCAGATTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑tt                                                                                                                      >  2:277230/1‑149 (MQ=255)
                                gTACAGCAGATTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTg                                                                                                                     <  3:352803/149‑1 (MQ=255)
                                 tACTGCTGATTTGACGTTTGGTTTGCGAGGTTTCTGCTCCCTTAATATTTCTGACTTTACGGATGCTTTTCCGCGGGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTg                                                                                                                     <  6:205158/148‑1 (MQ=255)
                                    agcagATTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTcc                                                                                                                 >  7:77730/1‑149 (MQ=255)
                                       agATTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCt                                                                                                              <  3:70327/149‑1 (MQ=255)
                                        gATTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCtt                                                                                                             <  3:157646/149‑1 (MQ=255)
                                         aTTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGTAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTc                                                                                                            <  5:74028/149‑1 (MQ=255)
                                         aTTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCtt                                                                                                             >  1:178373/1‑148 (MQ=255)
                                           ttGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTc                                                                                                            >  1:264944/1‑147 (MQ=255)
                                                            aGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTc                                                                                          >  3:41816/1‑148 (MQ=255)
                                                              gTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGat                                                                                       <  8:122883/149‑1 (MQ=255)
                                                              gTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTACGCGCGCTGCAGTGTTGATTCCTCAACCTGACGGAGGCGTAATCTCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGTTCCGCTTCAAGCTCAATGGTCAGGTCGat                                                                                       <  7:270636/149‑1 (MQ=255)
                                                                  cAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATg                                                                                    <  2:264944/148‑1 (MQ=255)
                                                                   aGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTc                                                                                  <  8:77730/149‑1 (MQ=255)
                                                                         cATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTgg                                                                            >  4:210138/1‑149 (MQ=255)
                                                                           taataTTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGt                                                                                >  3:234541/1‑143 (MQ=255)
                                                                           taataTTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGt                                                                                <  4:234541/143‑1 (MQ=255)
                                                                                     gACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCAc                                                                >  2:99032/1‑149 (MQ=255)
                                                                                        tttACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGt                                                             >  3:145202/1‑149 (MQ=255)
                                                                                            cGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTc                                                         <  1:188263/149‑1 (MQ=255)
                                                                                             ggATGCTTTTCCGCGCGCTGCAGGGTTGATTCTTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTcc                                                        <  1:277230/149‑1 (MQ=255)
                                                                                             ggATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGCCTTcc                                                        >  5:192285/1‑149 (MQ=255)
                                                                                                 gCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCcaca                                                     >  8:319312/1‑148 (MQ=255)
                                                                                                  cTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACaaa                                                   <  4:41816/149‑1 (MQ=255)
                                                                                                   ttttCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAg                                                  >  7:104349/1‑149 (MQ=255)
                                                                                                    tttCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGt                                                 >  2:278136/1‑149 (MQ=255)
                                                                                                     ttCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCACTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGtt                                                >  4:314998/1‑149 (MQ=255)
                                                                                                     ttCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCACTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGtt                                                >  4:320165/1‑149 (MQ=255)
                                                                                                        cgcgcgCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGTTGaa                                             <  3:76044/149‑1 (MQ=255)
                                                                                                          cgcgcTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGTTGaaaa                                           <  1:278136/149‑1 (MQ=255)
                                                                                                           gcgcTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGTTGaa                                             >  5:368879/1‑146 (MQ=255)
                                                                                                           gcgcTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGTTGaa                                             <  6:368879/146‑1 (MQ=255)
                                                                                                              cTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGTTGAAAAGATg                                       <  8:9944/149‑1 (MQ=255)
                                                                                                                 cAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGTTGAAAAGATGAg                                     <  1:153775/148‑1 (MQ=255)
                                                                                                                   gggTTGATTCCTCAACCTGACGGAGGCGTAATCACTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGTTGAAAAGATGAGTGc                                  <  3:320165/149‑1 (MQ=255)
                                                                                                                   gggTTGATTCCTCAACCTGACGGAGGCGTAATCACTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGTTGAAAAGATGAGTGc                                  <  3:314998/149‑1 (MQ=255)
                                                                                                                       tGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGGTGAAAAGATGAGTGCAAGt                              >  2:305336/1‑149 (MQ=255)
                                                                                                                          ttCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTc                                                         >  7:146471/1‑119 (MQ=255)
                                                                                                                          ttCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTc                                                         <  8:146471/119‑1 (MQ=255)
                                                                                                                             cTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAATTTGAAAAGATGAGTGCAAGTGCAGCg                        >  3:286403/1‑149 (MQ=255)
                                                                                                                             cTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGTTGAAAAGATGAGTGCAAGTTCAGc                         <  7:319312/148‑1 (MQ=255)
                                                                                                                              tCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGTTGAAAAGATGAGTGCAAGTTCAGCgg                       >  7:91817/1‑149 (MQ=255)
                                                                                                                                aaCCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGTTGAAAAGATGAGTGCAAGTTCAGCggg                      <  4:145202/148‑1 (MQ=255)
                                                                                                                                    tGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGTTGAAAAGATGAGTGCAAGTTCAGCGGGGCGTg                 >  3:194422/1‑149 (MQ=255)
                                                                                                                                       cGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGTTGAAAAGATGAGTGCAAGTTCAGCGGGGCGTGgaa              <  2:284477/149‑1 (MQ=255)
                                                                                                                                           ggCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGTTGAAAAGATGAGTGCAAGTTCAGCGGGGCGTGGAAGaaa          <  8:104349/149‑1 (MQ=255)
                                                                                                                                           ggCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGTTGAAAAGATGAGTGCAAGTTCAGCGGGGCGTGGAAGaaa          <  3:210138/149‑1 (MQ=255)
                                                                                                                                                    gCTGGTTTGTATAATGCATCGGAGCTGGGCTTTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGTTGAAAAGATGAGTGCAAGTTCAGCGGGGCGTGGAAGAAACAATGACg  >  1:43976/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AACGATTATGGTAGGGCGTGTAAAGCATGATGGTACAGCAGATTTGACGTTTGGTTTGCGAGGTTTCAGCTCCCATAATATTTCTGACTTTACGGATGCTTTTCCGCGCGCTGCAGGGTTGATTCCTCAACCTGACGGAGGCGTAATCGCTGGTTTGTATAATGCATCGGAGCTGGGC‑TTGGTCCGCTTCAAGCTCAATGGTGAGGTCGATATGTCGTTTGGGCAGGCCGGCACGGTCTTCCACAAAGTTGAAAAGATGAGTGCAAGTTCAGCGGGGCGTGGAAGAAACAATGACG  >  NC_002947/1301385‑1301680

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