Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,945,340 C→T 6.6% R6H (CGT→CAT)  ruvB ← ATP‑dependent DNA helicase, component of RuvABC resolvasome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,945,3400CT 6.6% 186.1 / 2.7 60R6H (CGT→CAT) ruvBATP‑dependent DNA helicase, component of RuvABC resolvasome
Reads supporting (aligned to +/- strand):  ref base C (30/26);  new base T (2/2);  total (32/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 = 4.23e-01

GAGGGCATCGCCGCGCAGTTTCGCTGCTTTGATGAAAATCTCCATCTGTGAACGAACCTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGGCTTGCTTCTTGTGGTTTATAGCCCAGCGCCACCAGCGCGG  >  NC_000913/1945192‑1945483
                                                                                                                                                    |                                                                                                                                               
gAGGGCATCGCCGCGCAGTTTCGCTGCTTTGATGAAAATCTCCATCTGTGAACGAACCTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGAc                                                                                                                                                 <  2:245346/149‑1 (MQ=255)
 aGGGCATCGCCGCGCAGTTTCGCTGCTTTGATGAAAATCTCCATCTGTGAACGAACCTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACg                                                                                                                                                <  4:294507/149‑1 (MQ=255)
                    tCGCTGCTTTGATGAAAATCTCCATCTGTGAACGAACCTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCtt                                                                                                                             >  2:61326/1‑149 (MQ=255)
                      gctgctTTGATGAAAATCTCCATCTGTGAACGAACCTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCttt                                                                                                                            <  4:232477/148‑1 (MQ=255)
                         gctTTGATGAATATCTCCATCTGTGAACGAACCTGCGGCTGACCAACATACTCTTCCAGTAAATTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTCCCGGCAGAAGTCAGACGGTCTGCTTCAATCATCCTTTACCt                                                                                                                        >  2:155486/1‑149 (MQ=255)
                                          cATCTGTGAACGAACCTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAg                                                                                                       >  5:143573/1‑149 (MQ=255)
                                           aTCTGTGAACGAAACTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAgg                                                                                                      <  2:228766/149‑1 (MQ=255)
                                                tgAACGAACCTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCt                                                                                                   <  2:20690/147‑1 (MQ=255)
                                                  aacgaacCTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGCCGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTcgc                                                                                               >  1:214835/1‑149 (MQ=255)
                                                        ccTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAAtt                                                                                          >  2:225583/1‑148 (MQ=255)
                                                          tGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTa                                                                                         >  4:94991/1‑147 (MQ=255)
                                                             ggCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGtt                                                                                    >  7:322909/1‑149 (MQ=255)
                                                              gCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGttt                                                                                   <  8:322909/149‑1 (MQ=255)
                                                                  aCCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCa                                                                                 <  3:355246/147‑1 (MQ=255)
                                                                    cAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAg                                                                                                       >  8:376947/1‑123 (MQ=255)
                                                                    cAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAg                                                                                                       <  7:376947/123‑1 (MQ=255)
                                                                     aaCATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCActgct                                                                            <  1:65838/149‑1 (MQ=255)
                                                                      aCATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCActgctg                                                                           >  8:100642/1‑149 (MQ=255)
                                                                          aCTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGt                                                                       >  3:287452/1‑149 (MQ=255)
                                                                          aCTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGt                                                                       >  3:198266/1‑149 (MQ=255)
                                                                          aCTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGt                                                                       >  7:214360/1‑149 (MQ=255)
                                                                           ctctTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTc                                                                      <  8:214360/149‑1 (MQ=255)
                                                                            tctTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGt                                                                       >  3:193456/1‑147 (MQ=255)
                                                                              ttCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAgg                                                                   >  2:263780/1‑149 (MQ=255)
                                                                              ttCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCCCGTAGGCCTTCGCGAATTAAAGATTCACTGCTTGCGTCAgt                                                                   >  2:185335/1‑148 (MQ=255)
                                                                                ccAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAg                                                                    <  7:186547/146‑1 (MQ=255)
                                                                                ccAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGgc                                                                 >  4:348898/1‑149 (MQ=255)
                                                                                 cAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGgc                                                                 >  1:114457/1‑148 (MQ=255)
                                                                                  aGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGgcga                                                               <  4:198266/149‑1 (MQ=255)
                                                                                  aGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGgc                                                                 >  3:331086/1‑147 (MQ=255)
                                                                                   gTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGATGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCa                                                                                 >  1:22138/1‑130 (MQ=255)
                                                                                   gTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGATGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCa                                                                                 >  1:307227/1‑130 (MQ=255)
                                                                                   gTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGATGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCa                                                                                 <  2:307227/130‑1 (MQ=255)
                                                                                   gTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGATGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCa                                                                                 <  2:22138/130‑1 (MQ=255)
                                                                                          ggggCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGAttt                                                        <  2:114457/148‑1 (MQ=255)
                                                                                           gggCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTg                                                      <  4:193456/149‑1 (MQ=255)
                                                                                             gCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCt                                                    <  1:280868/149‑1 (MQ=255)
                                                                                              cGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTc                                                   <  3:348898/149‑1 (MQ=255)
                                                                                                     cgcgATCTTGTTCATGTTACGTCAAAGTGGTACAGGCATATATCTGACGGTCTTCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGTCTTCGCGAATTAAAGTTTCACTGCTGGCGTCATGGCGAGCGATTTTGCTCACCATGc                                            <  1:185335/149‑1 (MQ=255)
                                                                                                     cgcgATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGc                                            <  1:61326/149‑1 (MQ=255)
                                                                                                              cTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGgcttgctt                                   >  4:276687/1‑149 (MQ=255)
                                                                                                              cTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGgc                                         <  4:331086/143‑1 (MQ=255)
                                                                                                                  aTCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGGCTTGCTTCTtg                               >  4:112910/1‑149 (MQ=255)
                                                                                                                  aTCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGGCTTGCTTCTtg                               >  3:335228/1‑149 (MQ=255)
                                                                                                                       ccGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGGCTTGCTTCTtgtg                             <  7:261681/146‑1 (MQ=255)
                                                                                                                       ccGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGGCTTGCTTCTtgtg                             <  7:291293/146‑1 (MQ=255)
                                                                                                                       ccGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGGCTTGCTTCTTGTgg                            <  4:335228/147‑1 (MQ=255)
                                                                                                                       ccGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGGCTTGCTTCTTGTgg                            <  7:100642/147‑1 (MQ=255)
                                                                                                                        cGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGATTTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGGCTTGCTTCTTGTGGtt                          >  2:374088/1‑148 (MQ=255)
                                                                                                                         ggCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGGCTTGCTTCTTGTGGTTta                        <  4:287452/149‑1 (MQ=255)
                                                                                                                                    aCCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGGCTTGCTTCTTGTGGTTTATAGCCCAgcgc             >  2:364539/1‑149 (MQ=255)
                                                                                                                                        gCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGGCTTGCTTCTTGTGGTTTATAGCCCAGCGccacc         >  8:216970/1‑149 (MQ=255)
                                                                                                                                        gCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGGCTTGCTTCTTGTGGTTTATAGCCCAGCGccac          >  8:153832/1‑148 (MQ=255)
                                                                                                                                           gAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGGCTTGCTTCTTGTGGTTTATAGCCCAGCGCCACCAgc      <  4:28621/149‑1 (MQ=255)
                                                                                                                                            aaaTCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGgcgagc                                                             <  3:11201/93‑1 (MQ=255)
                                                                                                                                            aaaTCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGgcgagc                                                             >  4:11201/1‑93 (MQ=255)
                                                                                                                                            aaaTCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGgcttg                                      <  6:134436/116‑1 (MQ=255)
                                                                                                                                            aaaTCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGGCTTGCTTCTTGTGGTTTATAGCCCAGCGCCACCAgcg     >  4:161695/1‑149 (MQ=255)
                                                                                                                                            aaaTCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGAGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGgcttg                                      >  5:134436/1‑116 (MQ=255)
                                                                                                                                               tCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGGCTTGCTTCTTGTGGTTTATAGCCCAGCGCCACCAGCGCgg  >  3:174609/1‑149 (MQ=255)
                                                                                                                                                cAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGGCTTGCTTCTTGTGGTTTATAGCCCAGCGCCACCAGCGCgg  >  3:333164/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GAGGGCATCGCCGCGCAGTTTCGCTGCTTTGATGAAAATCTCCATCTGTGAACGAACCTGCGGCTGACCAACATACTCTTCCAGTAATTTGGGGCGAATGGCGCGATCTGCTACATCTTCCGGCAAAGTGGTACCGGCAGAAATCAGACGGTCTGCTTCAATCATCCTTTACCTCATAACGCGGCGCGTAGGGCTTCGCGAATTAAAGTTTCACTGCTGGCGTCAGGGCGAGCGATTTTGCTCACCATGCGGCTTGCTTCTTGTGGTTTATAGCCCAGCGCCACCAGCGCGG  >  NC_000913/1945192‑1945483

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