Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,520,513 T→C 65.3% D64G (GAT→GGT)  hofQ ← DNA catabolic putative fimbrial transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,520,5130TC65.3% 95.0 / 71.8 85D64G (GAT→GGT) hofQDNA catabolic putative fimbrial transporter
Reads supporting (aligned to +/- strand):  ref base T (16/13);  new base C (29/26);  total (45/40)
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.99e-01

TGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGG  >  NC_000913/3520414‑3520610
                                                                                                   |                                                                                                 
tGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACAt                                                                                                   >  2:33330/1‑100 (MQ=255)
 gCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATc                                                                                                  <  1:1910832/100‑1 (MQ=255)
  ccAGGCAATGGAATGCACTGAGAGAATGTTGCCTCCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCTAGGCAACACCt                                                                                                 <  2:1178733/100‑1 (MQ=255)
   cAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTg                                                                                                <  1:1561278/100‑1 (MQ=255)
     ggCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGtt                                                                                              <  1:33330/100‑1 (MQ=255)
      gCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTa                                                                                             <  2:747782/100‑1 (MQ=255)
      gCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACGACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTa                                                                                             <  2:1203051/100‑1 (MQ=255)
        aaTGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGa                                                                                           <  1:2206511/100‑1 (MQ=255)
                gCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGa                                                                                   >  1:1852937/1‑100 (MQ=255)
                 cACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGac                                                                                  >  1:2007463/1‑100 (MQ=255)
                 cACTGAGAGAAGGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGac                                                                                  >  2:1170757/1‑100 (MQ=255)
                  aCTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGaca                                                                                 >  1:2028985/1‑100 (MQ=255)
                   cTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGacac                                                                                >  1:1724764/1‑100 (MQ=255)
                     gagagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGTTTCACTTGTTAGATGTAAGGACACCg                                                                              <  2:2085729/100‑1 (MQ=255)
                      agagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACCGTCTGTAGTGCCTCCTTCCAGGGAACATCTGTTAGATGTAACGACACCGt                                                                             >  2:1888360/1‑100 (MQ=255)
                         gaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTTGTGCCTGCTTACAGGGAACACCTGTTAGATGTAACGACACCGTAc                                                                           >  2:553127/1‑99 (MQ=255)
                            tGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCt                                                                       >  1:718501/1‑100 (MQ=255)
                              ttGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGa                                                                     <  2:1112713/100‑1 (MQ=255)
                                gCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACg                                                                   >  1:1023233/1‑100 (MQ=255)
                                gCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACg                                                                   >  1:1562906/1‑100 (MQ=255)
                                gCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACg                                                                   >  2:2126343/1‑100 (MQ=255)
                                      cTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGc                                                             >  1:1387573/1‑100 (MQ=255)
                                       tGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGc                                                             >  1:1139710/1‑99 (MQ=255)
                                       tGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCg                                                            >  1:782709/1‑100 (MQ=255)
                                       tGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCg                                                            <  1:1031450/100‑1 (MQ=255)
                                       tGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCg                                                            <  1:2126357/100‑1 (MQ=255)
                                        gCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGa                                                           >  2:754618/1‑100 (MQ=255)
                                        gCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGa                                                           >  2:1073019/1‑100 (MQ=255)
                                        gCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGa                                                           >  1:669937/1‑100 (MQ=255)
                                          cccgTTATCAGTCCGGCGCTTTGCACTACAGTTTGTAGTGCCCGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGaca                                                         <  1:1338894/98‑1 (MQ=255)
                                           gcgTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGacac                                                        <  1:2126178/100‑1 (MQ=255)
                                              ttATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTGCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGAc                                                     <  1:1799306/100‑1 (MQ=255)
                                               tATCAGTCCGGCGCTTTTCACAACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGAc                                                     >  2:649466/1‑99 (MQ=255)
                                                  cAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAgg                                                 <  2:182840/100‑1 (MQ=255)
                                                  cAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAgg                                                 >  1:1540813/1‑100 (MQ=255)
                                                  cAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAgg                                                 <  2:2007449/100‑1 (MQ=255)
                                                  cAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAgg                                                 <  1:995946/100‑1 (MQ=255)
                                                   aGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGt                                                <  1:1621550/100‑1 (MQ=255)
                                                    gTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGtt                                               >  1:2081656/1‑100 (MQ=255)
                                                    gTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGtt                                               >  1:656876/1‑100 (MQ=255)
                                                    gTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGtt                                               <  1:423080/100‑1 (MQ=255)
                                                      ccGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCa                                             >  2:1221443/1‑100 (MQ=255)
                                                        ggCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAAc                                           <  1:490544/100‑1 (MQ=255)
                                                        ggCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAAc                                           <  2:1622232/100‑1 (MQ=255)
                                                           gcTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTtc                                        >  1:2247089/1‑100 (MQ=255)
                                                            gcaggcACTACAAACTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACAGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTtct                                       >  2:47622/6‑100 (MQ=255)
                                                              tttCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTcc                                     <  2:1947380/100‑1 (MQ=255)
                                                                tctcTACAGTTTGTATTACCAGCTTTAAGGGAAAATCTGTTAGATGTAACGACACCGTTCCGGTGAAGTCTGCCGACACAACGAGTATAAACTTCTcc                                     >  2:1991588/4‑98 (MQ=255)
                                                                tCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTg                                   <  1:2249806/100‑1 (MQ=255)
                                                                 cACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTg                                   >  1:2238209/1‑99 (MQ=255)
                                                                  aCTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGtt                                 >  1:1697174/1‑100 (MQ=255)
                                                                   cTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTc                                >  1:544229/1‑100 (MQ=255)
                                                                   cTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTc                                <  2:656876/100‑1 (MQ=255)
                                                                    tACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTca                               <  2:799254/100‑1 (MQ=255)
                                                                     aCAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcag                              >  2:2082730/1‑100 (MQ=255)
                                                                     aCAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcag                              <  2:1879600/100‑1 (MQ=255)
                                                                      cAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcag                              <  1:155628/99‑1 (MQ=255)
                                                                        gTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcagcca                           >  1:1098267/1‑100 (MQ=255)
                                                                        gTTTGTAGTGCCTGCTTCCAGGGAACACATGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcagcca                           >  1:97719/1‑100 (MQ=255)
                                                                           tGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGATGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAgcg                        >  2:82679/1‑100 (MQ=255)
                                                                           tGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAgcg                        >  1:1216416/1‑100 (MQ=255)
                                                                             tAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGcc                      <  2:968305/100‑1 (MQ=255)
                                                                              aGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCt                     >  1:587117/1‑100 (MQ=255)
                                                                                tGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGc                   <  1:1414700/100‑1 (MQ=255)
                                                                                 gCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCa                  >  2:94626/1‑100 (MQ=255)
                                                                                   cTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAac                <  2:1756268/100‑1 (MQ=255)
                                                                                    tGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAaca               <  2:313884/100‑1 (MQ=255)
                                                                                      cTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGATGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACAcc             >  2:531184/1‑100 (MQ=255)
                                                                                      cTTCCAGGGAACAGCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGGGACCCGACCAGGGTCAACGTCTCCTGTTCAGCCAGCGCCTGCa                  <  2:606637/95‑1 (MQ=255)
                                                                                      cTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACAcc             <  2:1784233/100‑1 (MQ=255)
                                                                                      cTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACAcc             >  2:1615434/1‑100 (MQ=255)
                                                                                      cTTCCAGGGAACACCTGTTAAATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACAcc             >  2:1387096/1‑100 (MQ=255)
                                                                                       ttCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCt            >  1:159601/1‑100 (MQ=255)
                                                                                         ccAGGGAACATCTATTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGa          >  1:1524096/1‑100 (MQ=255)
                                                                                          cAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAg         <  2:1911474/100‑1 (MQ=255)
                                                                                          cAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAg         <  1:942757/100‑1 (MQ=255)
                                                                                          cAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGa          <  2:2154431/99‑1 (MQ=255)
                                                                                          cAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAg         <  2:1852925/100‑1 (MQ=255)
                                                                                          cAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAg         >  2:486308/1‑100 (MQ=255)
                                                                                          cAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAg         >  2:236704/1‑100 (MQ=255)
                                                                                           aGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGc        >  2:677848/1‑100 (MQ=255)
                                                                                            gggAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCt       <  2:1023231/100‑1 (MQ=255)
                                                                                            gggAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCt       <  2:2247074/100‑1 (MQ=255)
                                                                                            gggAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCt       >  1:1394157/1‑100 (MQ=255)
                                                                                               aaCACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTAcc    <  2:481416/100‑1 (MQ=255)
                                                                                                 cATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCgg  >  1:1358746/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
TGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGG  >  NC_000913/3520414‑3520610

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: