Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,618,527 A→T 15.3% W84R (TGG→AGG)  hda ← ATPase regulatory factor involved in DnaA inactivation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,618,5270AT15.3% 220.1 / 29.0 85W84R (TGG→AGG) hdaATPase regulatory factor involved in DnaA inactivation
Reads supporting (aligned to +/- strand):  ref base A (43/29);  new base T (9/4);  total (52/33)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CGCCATCTCCCACAACTCATCGCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCG  >  NC_000913/2618429‑2618626
                                                                                                  |                                                                                                   
cGCCATCTCCCACAACTCATCGCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAg                                                                                                    <  2:2026656/100‑1 (MQ=255)
 gCCATCTCCCACAACTCATCGCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAgg                                                                                                   >  2:95160/1‑100 (MQ=255)
 gCCATCTCCCACAACTCATCGCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAgg                                                                                                   >  1:402426/1‑100 (MQ=255)
  ccATCTCCCACAACTCATCGCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGt                                                                                                  >  1:2057804/1‑100 (MQ=255)
       tcCCACAACTCATCGCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGtt                                                                                             <  2:1314056/100‑1 (MQ=255)
       tcCCACAACTCATCGCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGtt                                                                                             >  1:1522538/1‑100 (MQ=255)
        cccACAACTCATCGCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGttt                                                                                            >  1:2218220/1‑100 (MQ=255)
              aCTCATCGCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATcc                                                                                        <  1:640686/98‑1 (MQ=255)
                tCATCGCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCg                                                                                    <  1:245317/100‑1 (MQ=255)
                 cATCGCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCcggcgg                                                                                   >  2:2149497/1‑100 (MQ=255)
                 cATCGCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCgg                                                                                   <  1:347701/100‑1 (MQ=255)
                 cATCGCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCgg                                                                                   >  2:762368/1‑100 (MQ=255)
                 cATCGCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCg                                                                                    <  2:605959/99‑1 (MQ=255)
                  aTCGCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCggg                                                                                  <  2:2218194/100‑1 (MQ=255)
                   tCGCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTGCGGAACAAACCTGGTGCGTTTATCCAGCGGGa                                                                                 <  1:2279631/100‑1 (MQ=255)
                   tCGCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGa                                                                                 <  1:2006281/100‑1 (MQ=255)
                   tCGCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGa                                                                                 <  2:1812751/100‑1 (MQ=255)
                     gCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACa                                                                               >  1:1751131/1‑100 (MQ=255)
                     gCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACa                                                                               >  1:1284314/1‑100 (MQ=255)
                      ccTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACAt                                                                              >  2:301135/1‑100 (MQ=255)
                      ccTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACAt                                                                              <  2:754495/100‑1 (MQ=255)
                       cTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATa                                                                             >  2:652105/1‑100 (MQ=255)
                         gCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGc                                                                           <  2:739302/100‑1 (MQ=255)
                         gCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGc                                                                           >  2:2217957/1‑100 (MQ=255)
                           aaTACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGccc                                                                         <  2:1495132/100‑1 (MQ=255)
                           aaTACACTCAATGTTGTCGATAAAGACAAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCAAGCGGGACATAGccc                                                                         >  2:828482/1‑100 (MQ=255)
                            aTACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCa                                                                        <  1:1363859/100‑1 (MQ=255)
                                acTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCAcc                                                                      >  2:784663/1‑98 (MQ=255)
                                acTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGc                                                                    >  1:682312/1‑100 (MQ=255)
                                    aaTGTTGTCGATACAGACCAGCGACAAATGCTCCATTCCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCAt                                                                  <  2:1305486/98‑1 (MQ=255)
                                      tgttgtCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCTGGTGCGTTTATCCAGCGGGACATAGCTCACCGCATCGcc                                                              >  1:1179815/1‑100 (MQ=255)
                                       gttgtCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCa                                                             >  2:1987953/1‑100 (MQ=255)
                                       gttgtCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCa                                                             >  1:954051/1‑100 (MQ=255)
                                       gttgtCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCa                                                             >  1:2305304/1‑100 (MQ=255)
                                        ttgtCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCag                                                            <  2:1008386/100‑2 (MQ=255)
                                        ttgtCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCa                                                             <  1:668620/99‑1 (MQ=255)
                                         tgtCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACg                                                           >  2:1212410/1‑100 (MQ=255)
                                         tgtCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCAAGCGGGGCATAGCCCACCGCATCGCCa                                                             >  2:376708/1‑98 (MQ=255)
                                            cGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCt                                                         >  1:392371/1‑99 (MQ=255)
                                            cGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTg                                                        >  1:1324002/1‑100 (MQ=255)
                                             gATACAGACCAGTGCANAATGCTCCATACCGTANAGTCATTNCGTGNCCAGCCAGGTGGGTCAATCCAGCGGTGTATAGCCCACCACACCGCGACGCTGc                                                       >  2:2210209/1‑100 (MQ=255)
                                              aTACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCTGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGcn                                                      <  2:856739/100‑2 (MQ=255)
                                              aTACAGACCAGCGACAAATGCTCCATACCGTCAAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGCGACATAGCCCTCCGCATCGCCACGCTGc                                                       >  2:112136/1‑99 (MQ=255)
                                               tACAGACCATCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGa                                                     >  1:1638422/1‑100 (MQ=255)
                                               tACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCTGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCg                                                      >  1:40363/1‑99 (MQ=255)
                                               tACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGa                                                     <  2:776897/100‑1 (MQ=255)
                                                aCAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGAc                                                    <  1:968871/100‑1 (MQ=255)
                                                 cAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACa                                                   >  2:1048868/1‑100 (MQ=255)
                                                   gACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCTGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAAt                                                 >  2:763080/1‑100 (MQ=255)
                                                     ccAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTc                                               >  2:1788030/1‑100 (MQ=255)
                                                      cAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCTGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTcc                                              >  1:2070313/1‑100 (MQ=255)
                                                           aCAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCa                                         <  2:1733261/100‑1 (MQ=255)
                                                            cAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCaa                                        >  2:717787/1‑100 (MQ=255)
                                                             aaaTGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAg                                       <  2:1638402/100‑1 (MQ=255)
                                                               aTGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGcc                                     <  1:220025/100‑1 (MQ=255)
                                                                  cTCCATACCGTCGAGCACTTCCGGAACAAACCTGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCcgcg                                  <  1:1570879/100‑1 (MQ=255)
                                                                  cTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCcgcg                                  <  1:355354/100‑1 (MQ=255)
                                                                   tCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGt                                 >  2:2233119/1‑100 (MQ=255)
                                                                    ccATACCGTCGAGCACTTCCGGAACAAACCTGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTg                                <  1:1716125/100‑1 (MQ=255)
                                                                    ccATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCcgcg                                  >  1:1920684/1‑98 (MQ=255)
                                                                     cATACCGTCGAGCACTTCCGGAACAAACGTGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCTTgc                               >  1:1428111/1‑100 (MQ=255)
                                                                      aTACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTgca                              >  1:39112/1‑100 (MQ=255)
                                                                       tACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTgcag                             >  2:516797/1‑100 (MQ=255)
                                                                         ccGTCGAGCACTTCCGGAACAAACCTGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTgcagca                           >  2:1402035/1‑100 (MQ=255)
                                                                          cGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGcagcag                          >  1:2327911/1‑100 (MQ=255)
                                                                           gTCGAGCACTTCCGGAACAAACCTGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCACGCCGCGTGCAGCAGa                         >  1:2296860/1‑100 (MQ=255)
                                                                             cGAGCACTTCCGGAACAAACCTGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATg                       >  2:1122153/1‑100 (MQ=255)
                                                                               aGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGc                     <  2:1806859/100‑1 (MQ=255)
                                                                                  aCTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTgc                  >  1:1439628/1‑100 (MQ=255)
                                                                                   cTTCCGGAACAAACCTGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTgcg                 >  1:344476/1‑100 (MQ=255)
                                                                                   cTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTgcg                 >  1:2311698/1‑100 (MQ=255)
                                                                                    ttCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTgcgc                >  1:2381364/1‑100 (MQ=255)
                                                                                    ttCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTgcgc                >  2:1813944/1‑100 (MQ=255)
                                                                                     tCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGcc               >  1:383634/1‑100 (MQ=255)
                                                                                      ccGGAACAAACCAGGTGCTTTTATCCAGAGGGATATAGCCCATCTGATCGCTACGCTGCGACAATTCCGCGCAAGCCGTGTGCAGCAGATGGCTGCGccc              <  2:2229136/100‑1 (MQ=255)
                                                                                       cGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGcccc             <  2:541765/100‑1 (MQ=255)
                                                                                          atctaaCCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCcgcg          <  1:1990083/96‑1 (MQ=255)
                                                                                          aaCAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCcgcg          <  1:132651/100‑1 (MQ=255)
                                                                                           aCAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCgcgc         >  1:1408507/1‑100 (MQ=255)
                                                                                           aCAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCgcgc         >  1:1071049/1‑100 (MQ=255)
                                                                                            cAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGcc        >  2:289918/1‑100 (MQ=255)
                                                                                            cAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGcc        >  1:2094340/1‑100 (MQ=255)
                                                                                              aaCCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCcgcgag        >  2:1990061/1‑96 (MQ=255)
                                                                                                ccAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTcg    <  1:2025114/100‑1 (MQ=255)
                                                                                                ccAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTcg    <  1:1359615/100‑1 (MQ=255)
                                                                                                  aGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTcgcg  >  2:1351542/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CGCCATCTCCCACAACTCATCGCCTGCAATACACTCAATGTTGTCGATACAGACCAGCGACAAATGCTCCATACCGTCGAGCACTTCCGGAACAAACCAGGTGCGTTTATCCAGCGGGACATAGCCCACCGCATCGCCACGCTGCGACAATTCCGCGCAAGCCGCGTGCAGCAGATGGCTGCGCCCCGCGCCTTCGCG  >  NC_000913/2618429‑2618626

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