Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 560,276 A→G 6.3% Q130Q (CAA→CAG secY → protein translocase subunit SecY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947560,2760AG 6.3% 186.9 / 2.2 63Q130Q (CAA→CAGsecYprotein translocase subunit SecY
Reads supporting (aligned to +/- strand):  ref base A (29/30);  new base G (2/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

ATGCCGTACATTTCGGCATCGATCATCATGCAGTTGATGACGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACTTTCGTGGCAGGCGCGATGTTCATGATGTGGCTGGGCGAGCAGATCACCGAGCGCGGTGTGG  >  NC_002947/560130‑560421
                                                                                                                                                  |                                                                                                                                                     
aTGCCGTACATTTCGGCATCGATCATCATGCAGTTGATGACGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAg                                                                                                                                                      <  5:144351/148‑1 (MQ=255)
 tGCCGTACATTTCGGCATCGATCATCATGCAGTTGATGACGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAGGc                                                                                                                                                     >  2:208314/1‑148 (MQ=255)
   ccGTACATTTCGGCATCGATCATCATGCAGTTGATGACGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCAt                                                                                                                                                  >  7:51785/1‑149 (MQ=255)
     gTACATTTCGGCATCGATCATCATGCAGTTGATGACGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTg                                                                                                                                                <  5:259188/149‑1 (MQ=255)
      tACATTTCGGCATCGATCATCATGCAGTTGATGACGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTgg                                                                                                                                               <  1:122584/149‑1 (MQ=255)
      tACATTTCGGCATCGATCATCATGCAGTTGATGACGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTgg                                                                                                                                               <  2:215866/149‑1 (MQ=255)
        cATTTCGGCATCGATCATCATGCAGTTGATGACGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGc                                                                                                                                             <  8:51785/148‑1 (MQ=255)
        cATTTCGGCATCGATCATCATGCAGTTGATGACGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGc                                                                                                                                             <  2:200344/148‑1 (MQ=255)
         aTTTCGGCATCGATCATCATGCAGTTGATGACGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCAt                                                                                                                                           >  8:9274/1‑149 (MQ=255)
         aTTTCGGCATCGATCATCATGCAGTTGATGACGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCAt                                                                                                                                           >  1:118830/1‑149 (MQ=255)
                       catcatGCAGTTGATGACGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCAATggacat                                                                                                                                           <  1:174459/136‑5 (MQ=255)
                       catcatGCAGTTGATGACGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCAATggacat                                                                                                                                           >  2:174459/1‑132 (MQ=255)
                           atGCAGTTGATGACGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAGGCCAATGGACATGCCAAtggcttggacc                                                                                                                          <  1:208314/148‑4 (MQ=255)
                           atGCAGTTGATGACGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAGGCCAATGGACATGCCAAtggcttggac                                                                                                                           >  7:140850/1‑145 (MQ=255)
                              cAGTTGATGACGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAGGCCAATGGACATGCCAAtggcttggaccagt                                                                                                                       <  8:140850/148‑7 (MQ=38)
                                gTTGATGACGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAgg                                                                                                                    <  7:9274/149‑1 (MQ=255)
                                   gatgaCGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCTCTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGggc                                                                                                                  <  5:14856/148‑1 (MQ=255)
                                       aCGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGggc                                                                                                                  >  1:303171/1‑144 (MQ=255)
                                       aCGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGggc                                                                                                                  <  2:303171/144‑1 (MQ=255)
                                           cggTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGttt                                                                                                         >  2:59388/1‑149 (MQ=255)
                                            ggTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGtttt                                                                                                        <  2:107450/149‑1 (MQ=255)
                                                    cGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGc                                                                                                <  7:57669/149‑1 (MQ=255)
                                                       aaCTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGc                                                                                                                                             >  4:193438/1‑101 (MQ=255)
                                                       aaCTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGc                                                                                                                                             <  3:193438/101‑1 (MQ=255)
                                                        aCTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTg‑‑‑g                                                                                         <  8:130742/149‑1 (MQ=255)
                                                         cTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTg                                                                                                     <  5:210006/139‑1 (MQ=255)
                                                         cTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTg                                                                                                     >  6:210006/1‑139 (MQ=255)
                                                          tGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGc                                                                                       >  2:304066/1‑149 (MQ=255)
                                                              gCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTcc                                                                                   <  4:270859/149‑1 (MQ=255)
                                                               cAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGcc                                                                                                                                                    <  1:190170/87‑1 (MQ=255)
                                                               cAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGcc                                                                                                                                                    >  2:190170/1‑87 (MQ=255)
                                                                aGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTcc                                                                                   <  5:185318/147‑1 (MQ=255)
                                                                  ttGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGt                                                                                                    >  7:89147/1‑131 (MQ=255)
                                                                  ttGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGt                                                                                                    <  8:89147/131‑1 (MQ=255)
                                                                     aagaagGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATAtcgt                                                                            >  2:138934/1‑149 (MQ=255)
                                                                          ggaaggTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTGAGTGCCTTCCATAtcgtc                                                                           >  1:111181/1‑148 (MQ=255)
                                                                          ggaaggTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTg                                                                       >  1:237771/1‑149 (MQ=255)
                                                                           gaaggTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTc                                                                                    >  1:277723/1‑135 (MQ=255)
                                                                           gaaggTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTc                                                                                    <  2:277723/135‑1 (MQ=255)
                                                                                  aaGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACt                                                               >  6:187641/1‑149 (MQ=255)
                                                                                     cTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACTTTc                                                            >  6:80987/1‑149 (MQ=255)
                                                                                        gCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACTTTCGTg                                                         >  6:325689/1‑149 (MQ=255)
                                                                                        gCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACTTTCGTg                                                         >  1:136058/1‑149 (MQ=255)
                                                                                                        ccAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACttt                                                             >  5:28727/1‑129 (MQ=255)
                                                                                                        ccAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACttt                                                             <  6:28727/129‑1 (MQ=255)
                                                                                                          aGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACTTTCGTGGCAGGTGCGATGTTCatg                                       <  2:237771/149‑1 (MQ=255)
                                                                                                          aGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACTTTCGTGGCAGGCGCGATGTTCat                                        >  4:143517/1‑148 (MQ=255)
                                                                                                                cccGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACTTTCGTGGCAGGCGCGATGTTCATGATGTgg                                 <  5:80987/149‑1 (MQ=255)
                                                                                                                 ccGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATAtc                                                                              <  8:208769/103‑1 (MQ=255)
                                                                                                                 ccGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATAtc                                                                              >  7:208769/1‑103 (MQ=255)
                                                                                                                 ccGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCa                                                                 >  5:54499/1‑116 (MQ=255)
                                                                                                                 ccGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCa                                                                 <  6:54499/116‑1 (MQ=255)
                                                                                                                       cGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACTTTCGTGGCAGGCGCGATGTTCATGATGTGGCTggg                            <  2:118830/147‑1 (MQ=255)
                                                                                                                        ggCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACTTTCGTGGCAGGCGCGATGTTCATGATGTGGCTGGGc                           >  4:135181/1‑147 (MQ=255)
                                                                                                                        ggCACCGATATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTCGCCTGGCCAACCACGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACTTTCGTGGCAGGCGCGAGGTTCATGCTGTGGCTGGGc                           >  6:23300/1‑147 (MQ=255)
                                                                                                                          cACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACTTTCGTGGCAGGCGCGATGTTCATGATGTGGCTGGGCGAGc                       <  1:138934/149‑1 (MQ=255)
                                                                                                                          cACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACTTTCGTGGCAGGCGCGATGTTCATGATGTGGCTGGGCGAGc                       <  3:135181/149‑1 (MQ=255)
                                                                                                                                   ccTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACTTTCGTGGCAGGCGCGATGTTCATGATGTGGCTGGGCGAGCAGATCAct               >  1:114203/1‑147 (MQ=255)
                                                                                                                                   ccTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACTTTCGTGGCAGGCGCGATGTTCATGATGTGGCTGGGCGAGCAGATCAcc               >  6:306007/1‑148 (MQ=255)
                                                                                                                                    cTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACATTCGTGGCAGGCGCGATGTTCATGATGTGGCTGGGCGAGCAGATCACCGa             >  6:104431/1‑149 (MQ=255)
                                                                                                                                             gTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACTTTCGTGGCAGGCGCGATGTTCATGATGTGGCTGGGCGAGCAGATCACCGAGCGCGgtgt    <  5:187641/149‑1 (MQ=255)
                                                                                                                                             gTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACTTTCGTGGCAGGCGCGATGTTCATGATGTGGCTGGGCGAGCAGATCACCGAGCGCGgtgt    <  5:325689/149‑1 (MQ=255)
                                                                                                                                                cAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTACTGTAGGCCTG‑‑‑CGCTTCCATATCGTCGCTGTTTCCACTATCGTGGCAGGCGCGATCATCATGATGTGGCTTGGCGAGCAGATCACCGACCCCGGTGTgg  >  7:251913/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ATGCCGTACATTTCGGCATCGATCATCATGCAGTTGATGACGGCGGTCAGCCCGCAACTGGAGCAGTTGAAGAAGGAAGGTGAAGCTGGCCGTCGCAAGATCAGCCAGTACACCCGCTACGGCACCGTTATCCTGGCACTGGTCCAAGCCATTGGCATGTCCATTGGCCTGGCCAACCAGGGCGTGGCGTTTTCTGTAGGCCTG‑‑‑GGCTTCCATATCGTCGCTGTTTCCACTTTCGTGGCAGGCGCGATGTTCATGATGTGGCTGGGCGAGCAGATCACCGAGCGCGGTGTGG  >  NC_002947/560130‑560421

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: