Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,024,749 G→C 7.6% S31W (TCG→TGG)  PP_2640 ← GNAT family acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,024,7490GC 7.6% 210.7 / 7.0 66S31W (TCG→TGG) PP_2640GNAT family acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (33/28);  new base C (2/3);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CCTTCCTTGAGCTCGATCAGCCCCGCGATAGCGCCTTCGACAAAACAGGCCAGCATGAGGTTGTCACCCTGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTC  >  NC_002947/3024617‑3024896
                                                                                                                                    |                                                                                                                                                   
ccttccttGAGCTCGATCAGCCCCGCGATAGCGCCTTCGACAAAACAGGCCAGCATGAGGTTGTCACCCTGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGAct                                                                                                                                    >  1:408758/1‑150 (MQ=255)
  ttccttGAGCTCGATCAGCCCCGCGATAGCGCCTTCGACAAAACAGGCCCGCATGCGGTTGTCCCCCTGCCTGCTCCCGGTAAACGGCGCTGGTGCTGGGACGGGTGTAAAGGTTTCGATACCTCGCGCCCACACTGTCGgcgcaaggta                                                                                                                                  <  2:1005706/150‑7 (MQ=255)
      ttGAGCTCGATCAGCCCCGCGATAGCGCCTTCGACAAAACAGGCCAGCATGAGGTTGTCACCCTGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGActgcctgc                                                                                                                              >  2:357768/1‑150 (MQ=255)
            tCGATCAGCCCCGCGATAGCGCCTTCGACAAAACAGGCCAGCATGAGGTTGTCACCCTGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAg                                                                                                                        <  2:1373737/150‑1 (MQ=255)
                 cAGCCCCGCGATAGCGCCTTCGACAAAACAGGCCAGCATGAGGTTGTCACCCTGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATc                                                                                                                   <  1:801620/150‑1 (MQ=255)
                  aGCCCCGCGATAGCGCCTTCGACAAAACAGGCCAGCATGAGGTTGTCACCCTGCATGCGCCCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATcc                                                                                                                  <  2:1478460/150‑1 (MQ=255)
                       cgcgATAGCGCCTTCGACAAAACAGGCCAGCATGAGGTTGTCACCCTGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATgac                                                                                                             <  2:408758/150‑3 (MQ=255)
                                gcCTTCGACAAAACAGGCCAGCATGAGGTTGTCACCCTGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCa                                                                                                    <  1:716149/150‑1 (MQ=255)
                                gcCTTCGACAAAACAGGCCAGCATGAGGTTGTCACCCTGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCa                                                                                                    <  1:357768/150‑1 (MQ=255)
                                            aCAGGCCAGCATGAGGTTGTCACCCTGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAg                                                                                        <  1:306287/150‑1 (MQ=255)
                                             aaGGCCAGCATGAGGTTGTCACCCTGCATGCGCTAGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAgg                                                                                       <  2:922909/149‑1 (MQ=255)
                                                gCCAGCATGAGGCTGTCACCCTGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGCCAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGGGgtcg                                                                                    >  2:402739/1‑150 (MQ=255)
                                                  cAGCATGAGGTTGTCACCCTGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGgtcgtc                                                                                  >  2:643228/1‑150 (MQ=255)
                                                  cAGCATGAGGTTGTCACCCTGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCCACACTGTCGgcgcaag                                                                                                                                     <  1:945968/99‑4 (MQ=255)
                                                  cAGCATGAGGTTGTCACCCTGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCCACACTGTCGgcgcaag                                                                                                                                     <  1:945421/99‑4 (MQ=255)
                                                  cAGCATGAGGTTGTCACCCTGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCCACACTGTCGgcgcaag                                                                                                                                     >  2:945421/1‑96 (MQ=255)
                                                  cAGCATGAGGTTGTCACCCTGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCCACACTGTCGgcgcaag                                                                                                                                     >  2:945968/1‑96 (MQ=255)
                                                      aTGAGGTTGTCACCCTGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAggg                                                                              <  1:1294487/150‑1 (MQ=255)
                                                                 aCCCTGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAg                                                                                        <  2:1295708/129‑1 (MQ=255)
                                                                 aCCCTGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAg                                                                                        >  1:1295708/1‑129 (MQ=255)
                                                                     tGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACt                                                               >  2:1152125/1‑150 (MQ=255)
                                                                        aTGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGgtc                                                                                     <  1:412694/125‑1 (MQ=255)
                                                                        aTGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGgtc                                                                                     >  2:412694/1‑125 (MQ=255)
                                                                            gcTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATgg                                                        >  1:520584/1‑150 (MQ=255)
                                                                            gcTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATgg                                                        >  2:92335/1‑150 (MQ=255)
                                                                            gcTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTGGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATgg                                                        >  2:229458/1‑150 (MQ=255)
                                                                             cTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATg                                                         >  1:275324/1‑148 (MQ=255)
                                                                             cTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATg                                                         <  2:275324/148‑1 (MQ=255)
                                                                                  cAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCg                                                  >  1:375184/1‑150 (MQ=255)
                                                                                   aaaCGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGa                                                                                         >  2:13306/1‑110 (MQ=255)
                                                                                   aaaCGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGa                                                                                         <  1:13306/110‑1 (MQ=255)
                                                                                   aaaCGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGt                                                 >  1:859298/1‑150 (MQ=255)
                                                                                   aaaCGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGt                                                 >  2:947627/1‑150 (MQ=255)
                                                                                    aaCGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTc                                                <  1:92335/150‑1 (MQ=255)
                                                                                    aaCGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTc                                                <  1:643228/150‑1 (MQ=255)
                                                                                          tGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAgcgtg                                                                     >  1:911771/1‑123 (MQ=255)
                                                                                          tGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAgcgtg                                                                     <  2:911771/123‑1 (MQ=255)
                                                                                              gTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGTTAATCGTCCTTGTCGGgt                                      <  1:947627/150‑1 (MQ=255)
                                                                                                 cTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTc                                   >  1:1070395/1‑150 (MQ=255)
                                                                                                 cTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTc                                   >  2:1421941/1‑150 (MQ=255)
                                                                                                  tGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTcc                                  <  2:520584/150‑1 (MQ=255)
                                                                                                       ccTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTg                             <  1:402739/150‑1 (MQ=255)
                                                                                                            gTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCGTGGGGGAGGTCGTCAGGGGGGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGcc                        >  1:110034/1‑150 (MQ=255)
                                                                                                              aaaGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCATGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCa                      >  2:1067707/1‑150 (MQ=255)
                                                                                                                  gTCTCGATACCTTGCGCCGACAGTGTGGGCGTGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATgg                                                        >  2:299736/1‑112 (MQ=255)
                                                                                                                  gTCTCGATACCTTGCGCCGACAGTGTGGGCGTGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATgg                                                        <  1:299736/112‑1 (MQ=255)
                                                                                                                      cGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAAcc              >  2:805892/1‑150 (MQ=255)
                                                                                                                       gATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCt             >  1:12681/1‑150 (MQ=255)
                                                                                                                       gATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCATGCGTACTTGCATGGTACTCGTCCTTGCCGCGTGTACTTCTGTTCCCAACTTTACCCt             >  2:581004/1‑150 (MQ=255)
                                                                                                                         tACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCt             >  2:1627001/1‑148 (MQ=255)
                                                                                                                         tACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCt             <  1:1627001/148‑1 (MQ=255)
                                                                                                                         tACCTTGCGCCGACAGTGTAGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGgtc                                                                                     <  2:1373046/76‑1 (MQ=255)
                                                                                                                         tACCTTGCGCCGACAGTGTAGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGgtc                                                                                     >  1:1373046/1‑76 (MQ=255)
                                                                                                                            cTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTg                             >  2:211206/1‑129 (MQ=255)
                                                                                                                            cTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTg                             <  1:211206/129‑1 (MQ=255)
                                                                                                                            cTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAAc               <  1:1484983/143‑1 (MQ=255)
                                                                                                                            cTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAAc               >  2:1484983/1‑143 (MQ=255)
                                                                                                                             ttGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTg                                                                                              <  2:1173968/63‑1 (MQ=255)
                                                                                                                             ttGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTg                                                                                              >  1:1173968/1‑63 (MQ=255)
                                                                                                                             ttGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATc                                                   <  1:1187000/106‑1 (MQ=255)
                                                                                                                             ttGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATc                                                   >  2:1187000/1‑106 (MQ=255)
                                                                                                                             ttGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGc                         >  1:587080/1‑132 (MQ=255)
                                                                                                                             ttGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGc                         <  2:587080/132‑1 (MQ=255)
                                                                                                                                  ccGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTc  >  2:130703/1‑150 (MQ=255)
                                                                                                                                   cGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGgt                                      <  1:1597172/113‑1 (MQ=255)
                                                                                                                                   cGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGgt                                      >  2:1597172/1‑113 (MQ=255)
                                                                                                                                    gACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTc                                   <  1:400325/115‑1 (MQ=255)
                                                                                                                                    gACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTc                                   >  2:400325/1‑115 (MQ=255)
                                                                                                                                    |                                                                                                                                                   
CCTTCCTTGAGCTCGATCAGCCCCGCGATAGCGCCTTCGACAAAACAGGCCAGCATGAGGTTGTCACCCTGCATGCGCTCGGCAAACGCCTGGGGTGCTGACACCTGGGTAAAGGTCTCGATACCTTGCGCCGACAGTGTGGGCGCGACTGCCTGCATGAAGGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTC  >  NC_002947/3024617‑3024896

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