Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 933,993 C→T 6.6% P257L (CCG→CTG)  ftsK → DNA translocase at septal ring sorting daughter chromsomes

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913933,9930CT 6.6% 149.4 / 2.8 61P257L (CCG→CTG) ftsKDNA translocase at septal ring sorting daughter chromsomes
Reads supporting (aligned to +/- strand):  ref base C (25/31);  new base T (2/2);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

GATGATACCTGGGTCGATGAAGATGAGTATGAAGACGACGAAGAGTATGAAGATGAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACG  >  NC_000913/933845‑934139
                                                                                                                                                    |                                                                                                                                                  
gatgatACCAGGGTCGATGAAGATGAGTATGAAGACGACGAAGAGTATGAAGATGAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTATTCGCGGCGCGCTAGCGCGTCGTAAACTGTTGGCGGAAAAATTCATTAATcc                                                                                                                                                    >  1:810995/1‑149 (MQ=255)
catgatACCGGGCTCGAAGACCATGAGGATGAAGACCCCGAAGCGTAGGAAGATGTAAAGCACGGCAAACAGCCTGAATCACGCCGCGCCCGTAATCTTCGCGGCGCGCTCGCGCGGCGTCAAGGGTTGGCGGAAAAATTCATTAATcc                                                                                                                                                    <  2:597378/148‑1 (MQ=255)
  tgatACCTGGGTCGATGAAGATGAGTATGAAGACGACGAAGAGTATTAAGATTAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTTTTCTTCGCGGCGCGCTAGCGCGTCTTAAACTGTTGGCGGTAAAATTCATTAATCCGa                                                                                                                                                  <  1:520994/149‑1 (MQ=255)
   gatACCTGGGTCGTGGAATATGATTATGAAGACGACGAAGAGTATGAAGTTGAAAATCACGGCAAAGAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGAt                                                                                                                                                 <  1:158250/149‑1 (MQ=255)
        cTGGGTCGATGAAGATGAGTATGAAGACGACGAAGAGTATGAAGATGAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGggc                                                                                                                                            >  1:436521/1‑149 (MQ=255)
          gggTCGATGAAGATGAGTATGAAGACGACGAAGAGTATGAAGATGAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTTAGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGAAAAAATTCATTAATACGATGGggcgg                                                                                                                                          >  1:996118/1‑149 (MQ=255)
           ggTCGATGAAGATGAGTATGAAGACGACGAAGAGTATGAAGATGAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGggcggc                                                                                                                                         <  2:378293/149‑1 (MQ=255)
               gATGAAGATGAGTATGAAGACGACGAAGAGTATGAAGATGAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTCGCGCGTCGTAAACCGTTGGCGGAAAAATTCAATAATCCGATGGGGCGGCAAAc                                                                                                                                     >  1:948658/1‑149 (MQ=255)
               gATGAAGATGAGTAGGAAGACGACGAAGAGTATGAAGATGAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAAc                                                                                                                                     <  1:444687/149‑1 (MQ=255)
                  taaGATGAGTATGAAGACGACGAAGAGTATGAAGATGAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGa                                                                                                                                  <  1:877070/148‑1 (MQ=255)
                        gAGTATGAAGACGACGAAGAGTATGAAGATGAAAATCACGGCAAACCGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGc                                                                                                                            >  2:737009/1‑149 (MQ=255)
                            aTGCAGACGACGAAGAGTATGAAGATGAACATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGCTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCgttg                                                                                                                        <  1:751554/149‑1 (MQ=255)
                            aTGAAGACGACGAAGAGTATGAAGATGAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCgttg                                                                                                                        <  1:510680/149‑1 (MQ=255)
                            aTGAAGACGACGAAGAGTATGAAGATGAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGGTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCgttg                                                                                                                        <  1:391220/149‑1 (MQ=255)
                                  acgacgAAGAGTATGAAGATGAAAATCACGGCAAACAGCATGCATCACGCCGTGACCGTATTCTTCGCGTCGCGCTAGCGCGTCGTAAACGGGTGGCGTAAAAAATCATTAATCACATGCGGCGTCAAACAGAATCTGCGTTGTACTcc                                                                                                                  >  2:345934/1‑149 (MQ=255)
                                   cgacgaAGAGTATGAAGATGAAAATGACGGCAAACTGCATGAATCACGCAGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCg                                                                                                                 <  1:186351/149‑1 (MQ=255)
                                     acgaAGAGTATGAAGATGAAAATCATGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTACGCGGCGCGCTAGCGCGTCGTAAACGGGTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCgga                                                                                                               >  2:308732/1‑148 (MQ=255)
                                            gTATGAAGATGAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCgga                                                                                                        >  1:938174/1‑149 (MQ=255)
                                            gTATGAAGATGAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCgga                                                                                                        <  1:891080/149‑1 (MQ=255)
                                                gAAGATGAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCggatgga                                                                                                    >  2:799063/1‑149 (MQ=255)
                                                gAAGATGAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCggatgga                                                                                                    <  1:450181/149‑1 (MQ=255)
                                                 aaGATGAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCggatggat                                                                                                   >  1:511030/1‑149 (MQ=255)
                                                  aGATGAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGATGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGCCGCTGCGATGTTCTCCGGTAAGCGgatggatg                                                                                                  >  2:451439/1‑149 (MQ=255)
                                                      gAAAATCACGGCAAACAGCATGTATCACGCCGTGCCCGGTTTGTTCGCGGCGCGCTAGCGCGTCGTACTCGGTTGGCGTAAAAATTAATTAATCCGTTGGGGCTTAAAACAGACGATGCGTTGTTCTCCGGTTAGCTGATTGaggatgt                                                                                              >  2:580835/1‑148 (MQ=255)
                                                      gAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTACTCTGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAtgatga                                                                                              >  1:698877/1‑149 (MQ=255)
                                                      gAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCTGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGAtgatga                                                                                              >  1:159354/1‑149 (MQ=255)
                                                          aTCACGGCAAACAGAATGAATCACGCCGTGCCAGTATTCTTAGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAAATAATTAATACGATGGGGCGGCAAACAGACGCTCCGTTGTTCTCCGGTAAGCGGATGGATGATGCCGaa                                                                                          >  2:720881/1‑149 (MQ=255)
                                                            cACGGCAATCAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGATAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGTTGGATGATGACGAaga                                                                                        >  2:404713/1‑149 (MQ=255)
                                                            cACGGCAAACAGCATGCATCACGCCGGGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAATCGGGTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCGGCGTTGTTCTCCGGTTAGCGGATGGATGATGACGCaga                                                                                        >  2:150454/1‑149 (MQ=255)
                                                            cACGGCAAACAGCATGAATCACGCCGAGCCCGTATTCTTCGCCGCGCGCTAGCGCGTCGTAAACGGTTGGCGTAAAAATTCATAAATCCGATGGGGCGGCAAACAGCCGCTGCGTTGTTCTCCGGTAACCGGATGGATGATGACGAaga                                                                                        >  2:763251/1‑149 (MQ=255)
                                                               ggCAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGCTGACGAAGAGAt                                                                                     >  2:212144/1‑149 (MQ=255)
                                                                 cAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTa                                                                                   <  1:53001/149‑1 (MQ=255)
                                                                 cAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGAGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTa                                                                                   <  1:118464/149‑1 (MQ=255)
                                                                      aGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTac                                                                              <  1:838664/149‑1 (MQ=255)
                                                                      aGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTACTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTac                                                                              >  1:948550/1‑149 (MQ=255)
                                                                         aTGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACt                                                                           >  2:645401/1‑149 (MQ=255)
                                                                              tCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCCAACAGACGCTGCGGTGTTCTCCGGTAAGCGGATGGATGAGGCCGAAGAGATGACCTACCCTGCAcg                                                                      >  2:1011726/1‑149 (MQ=255)
                                                                                aCGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGCCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCAcacg                                                                    >  1:764181/1‑149 (MQ=255)
                                                                                      tGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTgg                                                              <  2:5798/149‑1 (MQ=255)
                                                                                      tGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGCGTGGCGGAAAACTTCATTAATCCGATGGGGCCGCAAACAGACGCAGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTgg                                                              <  2:541122/149‑1 (MQ=255)
                                                                                      tGCCCGTATTCTTCGCGGCGCGCTAGCGCGGCGTAAACGGTTGGCGGAAAAATTCAATAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTgg                                                              <  2:515637/149‑1 (MQ=255)
                                                                                         ccGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTg                                                           >  1:204770/1‑149 (MQ=255)
                                                                                             aTTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGAGGATGGATGATGACGAAGAGATTAGCTACACTGCACGCGGTGTGGCTGCCgc                                                       >  1:585114/1‑148 (MQ=255)
                                                                                                 ttcGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGAGGTGTGGCTGCCGACCCg                                                   >  1:991152/1‑149 (MQ=255)
                                                                                                   cgcgGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCTGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGga                                                 <  2:159354/149‑1 (MQ=255)
                                                                                                   cgcgGCGCGCTAGCGCGGCGTAAACGGTGGGCGGAAAAATTCATTAATCTGATGGGGCGGCAAACAGACGCTGCGTTGTGCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGga                                                 <  2:698877/149‑1 (MQ=255)
                                                                                                      ggcgCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGgacga                                              <  2:13139/149‑1 (MQ=255)
                                                                                                          cgcTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGGTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTc                                          <  1:132331/149‑1 (MQ=255)
                                                                                                             tAGCGCGGCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCGGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTa                                       <  1:1009878/149‑1 (MQ=255)
                                                                                                                   gtcgtAAACGGTGGGCGGAAACATTCATTAATCCGCTGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCg                                 <  1:102110/149‑1 (MQ=255)
                                                                                                                   ggcggaaaCGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCg                                 <  1:44762/144‑1 (MQ=255)
                                                                                                                    tcgtAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCTGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCgg                                <  1:645401/149‑1 (MQ=255)
                                                                                                                    gcgtAAACGGTGGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCgg                                <  2:1668/148‑1 (MQ=255)
                                                                                                                     cgtAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGCCGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCggg                               <  2:524038/149‑1 (MQ=255)
                                                                                                                            ggTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCAGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGt                        <  2:902046/149‑1 (MQ=255)
                                                                                                                             gttttctgAATAATTCATTACTACTATGCGGCGGCATACAGCGGCTGCGTTGTTCTCCGGTAAGCGGTTGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTg                       <  1:345934/142‑1 (MQ=255)
                                                                                                                                ggCGGAAAACTTCATAAATCCGACGGGGCGGCAAACAGACGCGGCGGTGTTCTCCGGCAAGCGGAAGGATGATGACGAAGACACGACCTACACTGCACGCGGTGGGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCaa                    <  2:948658/149‑1 (MQ=255)
                                                                                                                                ggCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCaa                    <  1:737009/149‑1 (MQ=255)
                                                                                                                                            cATTAATCCGATGGGGCGGCAAACAGACGCTGCGTGGTTCTCAGGTAAGCGGATTGATGATGACGAAGAGATTACCTACACTGCACGCGGTTTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAat        <  2:810995/149‑1 (MQ=255)
                                                                                                                                             aTTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCCCGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAata       >  2:797301/1‑149 (MQ=255)
                                                                                                                                             aTTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGCCGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAata       >  2:242903/1‑149 (MQ=255)
                                                                                                                                              ttAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAatat      <  2:511030/149‑1 (MQ=255)
                                                                                                                                                  tCCGATGGGGCGGCAAACAGACGCTGAGTTGTTCTCCGGTAAGCGGATGGATGAAGACGAAGAGGTTACCTACACTGCACGCGGTGTGGCTGCCGACCCGTACGACGTCCTATTTTCGGGCATTCTTGCAACGCAGCCAGAATATGAcc  >  1:89744/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GATGATACCTGGGTCGATGAAGATGAGTATGAAGACGACGAAGAGTATGAAGATGAAAATCACGGCAAACAGCATGAATCACGCCGTGCCCGTATTCTTCGCGGCGCGCTAGCGCGTCGTAAACGGTTGGCGGAAAAATTCATTAATCCGATGGGGCGGCAAACAGACGCTGCGTTGTTCTCCGGTAAGCGGATGGATGATGACGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACG  >  NC_000913/933845‑934139

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