Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,445,295 C→G 7.0% G117R (GGG→CGG) ‡ nicE ← maleate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,445,2950CG 7.0% 221.9 / 7.6 72G117R (GGG→CGG) ‡nicEmaleate isomerase
Reads supporting (aligned to +/- strand):  ref base C (31/36);  new base G (2/3);  total (33/39)
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.78e-01

TCGGGGATTTCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGCGTCGGAC  >  NC_002947/4445147‑4445440
                                                                                                                                                    |                                                                                                                                                  
tCGGGGATTTCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCg                                                                                                                                                   >  2:1694346/1‑150 (MQ=255)
   gggATTTCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCTGTTAGCGGACTCATCTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGgcg                                                                                                                                                >  1:611382/1‑150 (MQ=255)
      aTTTCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTg                                                                                                                                             <  2:1103728/150‑1 (MQ=255)
        ttCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCGGGGCGCTGGt                                                                                                                                           <  1:985760/150‑1 (MQ=255)
        ttCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGt                                                                                                                                           >  1:392914/1‑150 (MQ=255)
          ccAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGa                                                                                                                                         >  1:763822/1‑150 (MQ=255)
              cgcgcgcCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCAc                                                                                                                                     <  1:1137261/150‑1 (MQ=255)
                cgcgcCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCg                                                                                                                                   <  1:1566401/150‑1 (MQ=255)
                 gcgcCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCgg                                                                                                                                  <  1:168637/150‑1 (MQ=255)
                    ccAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGc                                                                                                                               >  2:445323/1‑150 (MQ=255)
                    ccAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGc                                                                                                                               <  2:35208/150‑1 (MQ=255)
                    ccAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCATCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGc                                                                                                                               <  2:1459339/150‑1 (MQ=255)
                        tCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGt                                                                                                                           >  2:55305/1‑150 (MQ=255)
                           tGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGca                                                                                                                        >  2:1063096/1‑150 (MQ=255)
                                cTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCa                                                                                                                                      >  2:1539018/1‑131 (MQ=255)
                                cTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCa                                                                                                                                      <  1:1539018/131‑1 (MQ=255)
                                  cGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCTCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGt                                                                                                                 <  1:421673/150‑1 (MQ=255)
                                       ccTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTcggcg                                                                                                            <  1:1594720/150‑1 (MQ=255)
                                       ccTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTcggcg                                                                                                            >  2:908368/1‑150 (MQ=255)
                                             tCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCGGGGCGCTgg                                                                                                                                            >  1:984883/1‑112 (MQ=255)
                                             tCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCGGGGCGCTgg                                                                                                                                            <  2:984883/112‑1 (MQ=255)
                                             tCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGcc                                                                                                     >  1:613339/1‑150 (MQ=255)
                                               gcgGATGTCGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCctg                                                                                                   <  1:15757/150‑1 (MQ=255)
                                               gcgGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCctg                                                                                                   >  1:1898864/1‑150 (MQ=255)
                                               gcgGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCctg                                                                                                   >  1:1898809/1‑150 (MQ=255)
                                               gcgGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTAGATGTAGTCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCctg                                                                                                   >  1:19051/1‑150 (MQ=255)
                                                cgGATGTCGTTCACCCACACACAGGTCCCCGGCTCGATGTACGTCGGGTCTATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCctgc                                                                                                  <  2:1647893/150‑1 (MQ=255)
                                                cgGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCctgc                                                                                                  <  1:908368/150‑1 (MQ=255)
                                                   aTGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCtgctgc                                                                                               <  2:1510166/150‑1 (MQ=255)
                                                     gTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAg                                                                                             <  2:392914/150‑1 (MQ=255)
                                                        gTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGAAATGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGgcg                                                                                          >  1:1563189/1‑150 (MQ=255)
                                                         ttCACCACCAACTCGGTCAGTGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGgcgc                                                                                         <  2:111222/150‑1 (MQ=255)
                                                           caccacCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCAtt                                                                                                                   <  2:1762408/123‑1 (MQ=255)
                                                           caccacCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCAtt                                                                                                                   <  2:1762551/123‑1 (MQ=255)
                                                           caccacCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCAtt                                                                                                                   >  1:1762551/1‑123 (MQ=255)
                                                           caccacCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCAtt                                                                                                                   >  1:1762408/1‑123 (MQ=255)
                                                            accaccAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGGCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCttt                                                                                      >  1:372071/1‑150 (MQ=255)
                                                              caccaACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTGCAGAATCTGGACGCCTCTGTTGATCACCGGTGCCCGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTc                                                                                    >  1:798465/1‑150 (MQ=255)
                                                              caccaACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTc                                                                                    >  2:1775058/1‑150 (MQ=255)
                                                              caccaACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGATTTTc                                                                                    >  1:915532/1‑150 (MQ=255)
                                                              caccaACTCGGGCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTc                                                                                    <  2:1633757/150‑1 (MQ=255)
                                                                        gTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCgg                                                                          <  1:1775058/150‑1 (MQ=255)
                                                                         tCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGt                                                                         <  1:1090384/150‑1 (MQ=255)
                                                                              ggCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCACCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATag                                                                    <  1:124429/150‑1 (MQ=255)
                                                                               gCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGgcgc                                                                                         <  2:771800/128‑1 (MQ=255)
                                                                               gCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGgcgc                                                                                         >  1:771800/1‑128 (MQ=255)
                                                                                cTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATagcc                                                                  <  1:1288009/150‑1 (MQ=255)
                                                                                          ggCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATg                                                        <  2:611382/150‑1 (MQ=255)
                                                                                           gcgcGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATgg                                                       <  1:1063096/150‑1 (MQ=255)
                                                                                               gacgaTGGCAATGCGCTTGGCCTTCATCACATGCAAGGGCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCAt                                                   <  2:372071/150‑1 (MQ=255)
                                                                                                  gaTGGCAATGCGCTTGGCCTTCTTCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGAt                                                <  2:1115909/150‑1 (MQ=255)
                                                                                                   aTGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGTCAGTGCATCCTTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATc                                               <  2:1494545/150‑1 (MQ=255)
                                                                                                   aTGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCCTTGTCGGCGG‑TGGCCTGCTGCCGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATc                                               <  2:980181/150‑1 (MQ=255)
                                                                                                       cAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGGTCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCAATGATCGCca                                           >  1:1673613/1‑150 (MQ=255)
                                                                                                       cAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCca                                           <  2:798465/150‑1 (MQ=255)
                                                                                                        aaTGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCcac                                          >  2:1734488/1‑150 (MQ=255)
                                                                                                         aTGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCcaca                                         >  2:100396/1‑150 (MQ=255)
                                                                                                           gcgcTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGt                                                                         >  2:867625/1‑116 (MQ=255)
                                                                                                           gcgcTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGt                                                                         <  1:867625/116‑1 (MQ=255)
                                                                                                             gcTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAaggc                                     >  2:1770534/1‑150 (MQ=255)
                                                                                                               ttGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCAATGATCGCCACAaggcag                                   <  2:1673613/150‑1 (MQ=255)
                                                                                                                   ccTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCgg                                                                          <  2:1285444/107‑1 (MQ=255)
                                                                                                                   ccTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCgg                                                                          >  1:1285444/1‑107 (MQ=255)
                                                                                                                    cTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATagccc                                                                 <  1:1548781/115‑1 (MQ=255)
                                                                                                                    cTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATagccc                                                                 >  2:1548781/1‑115 (MQ=255)
                                                                                                                         tcaCATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCa                         >  2:1567886/1‑150 (MQ=255)
                                                                                                                            caTGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGca                      <  2:613339/150‑1 (MQ=255)
                                                                                                                                      cTCGACCAGCGCGCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGGCAGGGCGGC‑GAGAGCGCTTTTCAAATTGGCGGTTATAGCCCAGCCCCAGGGCCTTGATCGCCACAATGCAGGCATAGCCCAGCACATCGACCtt            <  1:1875668/150‑1 (MQ=255)
                                                                                                                                      cTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCtt            >  1:554459/1‑150 (MQ=255)
                                                                                                                                       atcaCCAGCGCCGGGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGa                                                 >  1:1702471/4‑112 (MQ=255)
                                                                                                                                       atcaCCAGCGCCGGGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGa                                                 <  2:1702471/109‑1 (MQ=255)
                                                                                                                                                gcCNCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGCGTCGGAc  >  2:1199473/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCGGGGATTTCCAGCGCGCGCCAGTCCTGCACCTCGAAGCCTTCCTCGCGGATGTAGTTCACCACCAACTCGGTCAGCGGCTTCATGTACGGCGCGACGATGGCAATGCGCTTGGCCTTCATCACATGCAAGGCCTCGACCAGCGCCCCGGCGCTGGTGATCACCGGTGCCAGTGCATCATTGTCGGCGG‑TGGCCTGCTGCAGGCGCTTTTCAGATTGGCGGTGATAGCCCAGCCCCATGGCCATGATCGCCACAAGGCAGGCATAGCCCAGCACATCGACCTTGGCGTCGGAC  >  NC_002947/4445147‑4445440

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