Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,649,531 T→C 7.7% H524R (CAT→CGT)  PP_1446 ← TonB‑dependent receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,649,5310TC 7.7% 150.6 / 2.1 52H524R (CAT→CGT) PP_1446TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base T (22/26);  new base C (2/2);  total (24/28)
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.72e-01

GTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGTCGAAGCGCAGC  >  NC_002947/1649383‑1649677
                                                                                                                                                    |                                                                                                                                                  
gTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCtcatg                                                                                                                                                   >  1:1195946/1‑150 (MQ=255)
 tCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCtcatga                                                                                                                                                  >  1:186617/1‑150 (MQ=255)
   gAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCtcatgata                                                                                                                                                >  1:525131/1‑149 (MQ=255)
         cTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCtcatgatcatgatc                                                                                                                                          <  1:186004/150‑1 (MQ=255)
              cTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTcatgatcatgatcatgat                                                                                                                                     >  1:1214451/1‑150 (MQ=255)
               tCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCatgatcatgatcatgatc                                                                                                                                    >  2:887348/1‑150 (MQ=255)
               tCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCatgatcatgatcatgatc                                                                                                                                    >  2:360969/1‑150 (MQ=255)
                   gTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGa                                                                                                                                <  2:1214451/150‑1 (MQ=255)
                      taGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCAGCAAGCTGCGAGCGTTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGAtcc                                                                                                                             <  2:749038/150‑1 (MQ=255)
                         tCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGAtcctcc                                                                                                                          >  1:452869/1‑150 (MQ=255)
                           aCCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGAtcctcctc                                                                                                                        >  2:188044/1‑150 (MQ=255)
                                      ccAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCACGAtcatgatcatgatcatgagg                                                                                                                              >  2:1050082/1‑131 (MQ=255)
                                      ccAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCACGAtcatgatcatgatcatgagg                                                                                                                              <  1:1050082/133‑3 (MQ=255)
                                       cAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCATCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTCGTATTCCGGGAAGCCCTTATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATg                                                                                                            <  2:236588/150‑1 (MQ=255)
                                         gcgcGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCt                                                                                                          >  2:157665/1‑150 (MQ=255)
                                                      gTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGgt                                                                                             >  1:1319893/1‑150 (MQ=255)
                                                       tgaCGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGgtg                                                                                            <  2:946826/148‑1 (MQ=255)
                                                           ggTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGc                                                                                        >  2:993348/1‑150 (MQ=255)
                                                            gTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGcc                                                                                       >  1:245470/1‑150 (MQ=255)
                                                             ttCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCg                                                                                      <  2:539920/150‑1 (MQ=255)
                                                             ttCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCg                                                                                      >  1:983259/1‑150 (MQ=255)
                                                                tcGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGc                                                                                        >  1:338829/1‑145 (MQ=255)
                                                                tcGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGc                                                                                        <  2:338829/145‑1 (MQ=255)
                                                                tcGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTg                                                                                   >  1:125434/1‑150 (MQ=255)
                                                                  ggCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGcc                                                                                 >  1:1268002/1‑150 (MQ=255)
                                                                  ggCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGcc                                                                                 >  2:1010218/1‑150 (MQ=255)
                                                                        cTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAg                                                                           <  2:891670/150‑1 (MQ=255)
                                                                         tGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAgg                                                                          <  1:160454/150‑1 (MQ=255)
                                                                           ccAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCACGAtcatgatcatgatcatgag                                                                                                                               >  1:356243/1‑94 (MQ=255)
                                                                           ccAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCACGAtcatgatcatgatcatgag                                                                                                                               <  2:356243/95‑2 (MQ=255)
                                                                             aGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGa                                                                      >  1:61165/1‑150 (MQ=255)
                                                                             aGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGa                                                                      <  2:125434/150‑1 (MQ=255)
                                                                               cGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATg                                                                    <  2:452869/150‑1 (MQ=255)
                                                                                  tCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGa                                                                                                                                >  1:863976/1‑87 (MQ=255)
                                                                                  tCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGa                                                                                                                                <  2:863976/87‑1 (MQ=255)
                                                                                  tCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAg                                                                 <  1:188044/150‑1 (MQ=255)
                                                                                    cTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGtt                                                               >  1:512976/1‑150 (MQ=255)
                                                                                    cTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGtt                                                               <  2:245470/150‑1 (MQ=255)
                                                                                         ccTCGATGCCGGAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGaa                                                          >  2:1308679/1‑150 (MQ=255)
                                                                                                   gTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTa                                                >  2:1100146/1‑150 (MQ=255)
                                                                                                   gTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAACCAGGCCGATGTAGTTGCGAAAGTGGCTGTa                                                <  1:887348/150‑1 (MQ=255)
                                                                                                          cgcgcCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGaacaac                                         <  1:993348/150‑1 (MQ=255)
                                                                                                              cccGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATAAGGCCGATGTAGTTGCCATTGTGGCTGTTAaacaacttta                                     <  2:525131/150‑5 (MQ=255)
                                                                                                                  cACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACt                                 <  1:1308679/150‑1 (MQ=255)
                                                                                                                   aCGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTg                                <  2:1195946/150‑1 (MQ=255)
                                                                                                                     gCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGcc                              >  1:310389/1‑150 (MQ=255)
                                                                                                                      ccTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGccc                             <  2:26925/150‑1 (MQ=255)
                                                                                                                          ggTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTtgt                         <  2:1268002/150‑1 (MQ=255)
                                                                                                                            tACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTgg                       <  1:501954/150‑1 (MQ=255)
                                                                                                                               tGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCt                                                   >  1:1233863/1‑119 (MQ=255)
                                                                                                                               tGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCt                                                   <  2:1233863/119‑1 (MQ=255)
                                                                                                                                 gTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGcc                  <  2:310389/150‑1 (MQ=255)
                                                                                                                                 gTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGcc                  >  1:522133/1‑150 (MQ=255)
                                                                                                                                         ggAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTa                                                                  >  2:972405/1‑94 (MQ=255)
                                                                                                                                         ggAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTa                                                                  <  1:972405/94‑1 (MQ=255)
                                                                                                                                            aGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGTCGAagc       <  1:1010218/150‑1 (MQ=255)
                                                                                                                                             gCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGTCGAagcg      >  2:1287084/1‑150 (MQ=255)
                                                                                                                                               ccTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGc                                                    <  1:1276786/102‑1 (MQ=255)
                                                                                                                                               ccTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGc                                                    >  2:1276786/1‑102 (MQ=255)
                                                                                                                                                 tcatgatcatgatcatgaTCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGTCGAagcgcagc  <  2:61165/150‑1 (MQ=255)
                                                                                                                                                 tcacgatcatgatcatgaTCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGTCGAagc       >  1:1272358/1‑145 (MQ=255)
                                                                                                                                                   atgatcatgatcatgatcGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGTCGAagcgcagc  <  2:593530/148‑1 (MQ=255)
                                                                                                                                                   atgatcatgatcatgatcGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGTCGAagcgcagc  >  1:593530/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGTCGAAGCGCAGC  >  NC_002947/1649383‑1649677

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