Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,867,576 G→T intergenic (‑154/‑152) ibpA ← / → yidQ small heat shock protein IbpA/DUF1375 domain‑containing putative lipoprotein YidQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,867,5760GT98.2% 206.2 / NA 55intergenic (‑154/‑152)ibpA/yidQsmall heat shock protein IbpA/DUF1375 domain‑containing putative lipoprotein YidQ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (29/25);  minor base C (1/0);  total (30/25)
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

ACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTCGGAGGATGCTCCAGCGTGATGTCCCATACCGGAGGGAAA  >  NC_000913/3867333‑3867820
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
aCTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGTGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGttt                                                                                                                                                                                                                                                <  1:180766/250‑1 (MQ=255)
      gctctggTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCt                                                                                                                                                                                                                                         <  1:439045/251‑1 (MQ=255)
       ctctggTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCtt                                                                                                                                                                                                                                        >  2:551303/1‑250 (MQ=255)
          tgggtgttatcgaaatGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAAGCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGCCCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGCACTTTTTTTTTCAGTTTATGTTCTAtt                                                                                                                                                                                                                                      <  2:241116/235‑1 (MQ=255)
                  ttCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGt                                                                                                                                                                                                                               >  1:427942‑M1/1‑248 (MQ=255)
                   tCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtag                                                                                                                                                                                                                             <  2:330348‑M1/250‑4 (MQ=255)
                       aGTGGTTAAACAAACGGTCAAATCCACTAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGAAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCCTAGtagatttg                                                                                                                                                                                                                        >  2:349795‑M1/1‑243 (MQ=255)
                              aaacaaacGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatg                                                                                                                                                                                                                 >  1:301754‑M1/1‑236 (MQ=255)
                                         aaaatCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCACCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcgg                                                                                                                                                                                                       <  2:187838‑M1/249‑26 (MQ=255)
                                              ccAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgc                                                                                                                                                                                                  >  1:529531‑M1/1‑220 (MQ=255)
                                                aaTAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggatt                                                                                                                                                                                                    >  1:18580‑M1/1‑218 (MQ=255)
                                                aaTAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggatt                                                                                                                                                                                                    <  2:18580‑M1/246‑29 (MQ=255)
                                                              aaaGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcgg                                                                                                                                                                                 <  2:73274‑M1/251‑48 (MQ=255)
                                                                  cGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGATCTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagat                                                                                                                                                                             <  1:586072‑M1/251‑52 (MQ=255)
                                                                  cGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCCGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcgg                                                                                                                                                                                 >  1:550689‑M1/1‑200 (MQ=255)
                                                                        aaaatCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctc                                                                                                                                                                       <  2:283628‑M1/250‑58 (MQ=255)
                                                                        aaaatCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGGCCGCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctc                                                                                                                                                                       <  1:444442‑M1/250‑58 (MQ=255)
                                                                         aaaTCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctcc                                                                                                                                                                      >  1:338390‑M1/1‑193 (MQ=255)
                                                                           aTCAAAGTTACGCATAATCAATAGCTCCTGCAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgc                                                                                                                                                                    >  2:444178‑M1/1‑191 (MQ=255)
                                                                           aTCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgc                                                                                                                                                                    >  1:122091‑M1/1‑191 (MQ=255)
                                                                              aaaGTTACACATAATCAATAGCTCCTAAAATCAGCGAGCATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTAAAGCCTTATCAGCGGACTTTTTTTTCCAGTTTATGGTCTTTTGCATAGtagatttgagggatgatacttgcgattggaagaaataaaaagtgcggagatgaccctagccct                                                                                                                                                                 >  2:430149‑M1/1‑188 (MQ=255)
                                                                               aaGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgccctt                                                                                                                                                                >  2:188820‑M1/1‑187 (MQ=255)
                                                                                 gTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgccctttt                                                                                                                                                              >  2:84049‑M1/1‑185 (MQ=255)
                                                                                    aCGCATAATCAATGGCCCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGTTCAGGTGGCCAGAACACCCATCAGACGCGTCTTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttctt                                                                                                                                                           <  1:106986‑M1/251‑70 (MQ=255)
                                                                                       cATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAACATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaaggtcggagatgatctccgcccttttcttat                                                                                                                                                         >  2:3920‑M1/1‑179 (MQ=255)
                                                                                       cATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttat                                                                                                                                                         >  1:49322‑M1/1‑179 (MQ=255)
                                                                                           cgttgctcttccgatcttcagcgagaatgtaagacCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagctt                                                                                                                                                    <  1:444034‑M1/234‑77 (MQ=255)
                                                                                                 aGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgcc                                                                                                                                              >  1:279087‑M1/1‑169 (MQ=255)
                                                                                                  gCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccg                                                                                                                                             >  2:396568‑M1/1‑168 (MQ=255)
                                                                                                            tCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagctt                                                                                                                                                    >  2:444034‑M1/1‑158 (MQ=255)
                                                                                                            tCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagc                                                                                                                                                      >  1:116258‑M1/1‑158 (MQ=255)
                                                                                                            tCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagc                                                                                                                                                      <  2:116258‑M1/232‑75 (MQ=255)
                                                                                                               gCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtga                                                                                                                                 >  1:466918‑M1/1‑155 (MQ=255)
                                                                                                                 gagaATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCAACAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacg                                                                                                                              <  2:122091‑M1/251‑99 (MQ=255)
                                                                                                                                tCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggccta                                                                                                                <  1:188820‑M1/250‑113 (MQ=255)
                                                                                                                                tCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggccta                                                                                                                <  1:66224‑M1/250‑113 (MQ=255)
                                                                                                                                        ggACAGGTAAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcat                                                                                                       >  1:389857‑M1/1‑130 (MQ=255)
                                                                                                                                                  ggtAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaa                                                                                             <  2:184806‑M1/251‑132 (MQ=255)
                                                                                                                                                    tAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTCTGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaat                                                                                            >  1:15459‑M1/1‑118 (MQ=255)
                                                                                                                                                             cacCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAACGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaag                                                                                  <  1:349795‑M1/251‑143 (MQ=255)
                                                                                                                                                               cccATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagat                                                                                <  1:550446‑M1/251‑145 (MQ=255)
                                                                                                                                                                ccATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagatt                                                                               >  1:526588‑M1/1‑106 (MQ=255)
                                                                                                                                                                        cgcgTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtagg                                                                       >  1:14755‑M1/1‑98 (MQ=255)
                                                                                                                                                                           gTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcct                                                                    <  2:523971‑M1/251‑157 (MQ=255)
                                                                                                                                                                                                           ggAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcat                                     <  1:18346‑M1/250‑188 (MQ=255)
                                                                                                                                                                                                           ggAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcacgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatg                                    >  2:382084‑M1/1‑63 (MQ=255)
                                                                                                                                                                                                            gAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatgg                                   >  1:525596‑M1/1‑62 (MQ=255)
                                                                                                                                                                                                                    cTTTTCAAGCATTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattacagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgacacccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcctaggcctgataagcgtagcgcatcaggcacgcccgcatggttagcgcc                           >  2:462167‑M1/1‑54 (MQ=255)
                                                                                                                                                                                                                                tATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataatggcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcg                <  2:370130‑M1/250‑209 (MQ=255)
                                                                                                                                                                                                                                 aTCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgat              >  1:199238‑M1/1‑41 (MQ=255)
                                                                                                                                                                                                                                   cAGCGGACTTGTTGTTTCAGTGTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgagctccgcccttttcttagagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgata             <  2:120612‑M1/250‑212 (MQ=255)
                                                                                                                                                                                                                                     gCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgatattc          >  1:220995‑M1/1‑37 (MQ=255)
                                                                                                                                                                                                                                     gCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgatattc          <  1:551303‑M1/251‑215 (MQ=255)
                                                                                                                                                                                                                                       ggACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgatattcgc        <  2:279087‑M1/251‑217 (MQ=255)
                                                                                                                                                                                                                                             ttttttTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgatattcgccgcttt  >  1:34279‑M1/1‑29 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
ACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTCGGAGGATGCTCCAGCGTGATGTCCCATACCGGAGGGAAA  >  NC_000913/3867333‑3867820

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