Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,488,424 G→T 6.9% intergenic (+24/+137) lptG → / ← yjgR lipopolysaccharide export ABC permease/DUF853 family protein with NTPase fold

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,488,4240GT 6.9% 235.0 / 4.6 72intergenic (+24/+137)lptG/yjgRlipopolysaccharide export ABC permease/DUF853 family protein with NTPase fold
Reads supporting (aligned to +/- strand):  ref base G (30/37);  new base T (2/3);  total (32/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.14e-01

CGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCT  >  NC_000913/4488277‑4488572
                                                                                                                                                   |                                                                                                                                                    
cGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCg                                                                                                                                                      >  5:129366/1‑148 (MQ=255)
     tGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCt                                                                                                                                                >  3:179333/1‑149 (MQ=255)
       gACCAGATCTTCGGCCCGCTGACGTTGGATTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTcc                                                                                                                                              <  2:101127/149‑1 (MQ=255)
        aCCATAGCTTCGGCCCTCTGACGTTGGATTATGGCATCCCGGCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGGTTATGAGAAAAGCGTAACGGAATTAACGAATAAAGGGAGCGAAAGCTccc                                                                                                                                             <  1:81614/149‑1 (MQ=255)
        aCCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTccc                                                                                                                                             >  5:77510/1‑149 (MQ=255)
          cAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCTTTCGCTCCCtt                                                                                                                                           <  5:36873/149‑1 (MQ=255)
          cAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCTTTCGCTCCCtt                                                                                                                                           <  7:15042/149‑1 (MQ=255)
          cAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCtt                                                                                                                                           >  4:71675/1‑149 (MQ=255)
          cAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCtt                                                                                                                                           >  3:66482/1‑149 (MQ=255)
            gATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTtta                                                                                                                                         >  2:15407/1‑149 (MQ=255)
            gATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTtta                                                                                                                                         <  3:71675/149‑1 (MQ=255)
            gATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTtta                                                                                                                                         >  3:174163/1‑149 (MQ=255)
             aTCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTTCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTtta                                                                                                                                         <  5:187929/148‑1 (MQ=255)
             aTCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTtta                                                                                                                                         >  6:187929/1‑148 (MQ=255)
                  cGGCCCGCTGTCGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTattatt                                                                                                                                    <  5:147197/148‑1 (MQ=255)
                   ggCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGt                                                                                                                                  <  3:99100/149‑1 (MQ=255)
                            aCGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAgt                                                                                                                         <  4:4131/149‑1 (MQ=255)
                              gTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGt                                                                                                                                  >  8:15933/1‑138 (MQ=255)
                              gTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGt                                                                                                                                  <  7:15933/138‑1 (MQ=255)
                              gTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAgtgt                                                                                                                       >  3:94230/1‑149 (MQ=255)
                               ttGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGtgtg                                                                                                                      >  2:178326/1‑149 (MQ=255)
                                  gTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGTGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCt                                                                                                                    <  8:232458/148‑1 (MQ=255)
                                     tATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCg                                                                                                                 <  7:8394/148‑1 (MQ=255)
                                      aTGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGtt                                                                                                               >  1:29731/1‑149 (MQ=255)
                                       tGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGtt                                                                                                               >  7:164509/1‑148 (MQ=255)
                                         gCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCa                                                                                                             <  4:174163/148‑1 (MQ=255)
                                           aTCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTc                                                                                                          <  8:15218/149‑1 (MQ=255)
                                            tCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTcc                                                                                                         >  3:164752/1‑149 (MQ=255)
                                            tCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTcc                                                                                                         >  7:185425/1‑149 (MQ=255)
                                                     tcatcGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGAGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAaat                                                                                                 >  4:62154/1‑148 (MQ=255)
                                                       atcGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAaataat                                                                                              <  5:104236/149‑1 (MQ=255)
                                                       atcGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAaataa                                                                                               <  8:164509/148‑1 (MQ=255)
                                                       atcGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAaataa                                                                                               <  6:77510/148‑1 (MQ=255)
                                                                   ttGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCaa                                                                                    <  7:11842/147‑1 (MQ=255)
                                                                   ttGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAat                                                                                   <  7:175950/148‑1 (MQ=255)
                                                                      ccAATCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAatata                                                                                <  7:114857/148‑1 (MQ=255)
                                                                       cAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATAtt                                                                              <  5:205190/149‑1 (MQ=255)
                                                                        aaGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTg                                                                             >  4:101484/1‑149 (MQ=255)
                                                                        aaGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTg                                                                             <  8:185425/149‑1 (MQ=255)
                                                                        aaGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTg                                                                             <  1:178799/149‑1 (MQ=255)
                                                                         aGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGc                                                                            <  2:29731/149‑1 (MQ=255)
                                                                         aGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTAATGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTg                                                                             >  5:197683/1‑148 (MQ=255)
                                                                          gcgcCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCa                                                                           <  4:164752/149‑1 (MQ=255)
                                                                                gCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAgtgggg                                                                     >  6:79168/1‑149 (MQ=255)
                                                                                    cttcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAgtggggtgg                                                                  <  5:79168/148‑1 (MQ=255)
                                                                                          aaTCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGa                                                            >  2:617/1‑148 (MQ=255)
                                                                                          aaTCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGa                                                            >  4:29014/1‑148 (MQ=255)
                                                                                                 ctttgGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGGTAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGAcgc                                                    <  5:2723/146‑1 (MQ=255)
                                                                                                        tGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCg                                             >  1:173721/1‑149 (MQ=255)
                                                                                                           tAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTc                                           <  1:178326/148‑1 (MQ=255)
                                                                                                            aaggagaAAATCGTAACCGAATCAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATAGTTAGCAAAGTGTGCTTCGTGCATAGCTGAAAAATAATTAATATTCCATATATTGCAGTGGGGTGGGCATGCTAAGACGCTCCAGGGTCg                                          >  6:183339/4‑148 (MQ=255)
                                                                                                            aaTGAGAAATTCGTAACCGTATGAACGAATAAAGGGAGCGAAAGCTCCCTTTattatt                                                                                                                                    <  3:10536/58‑1 (MQ=255)
                                                                                                            aaTGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTattatt                                                                                                                                    >  4:10536/1‑58 (MQ=255)
                                                                                                            aaTGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAATACGCGCCAGCGTCg                                          >  8:206573/1‑148 (MQ=255)
                                                                                                              tGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCAt                                       >  8:127641/1‑149 (MQ=255)
                                                                                                                  aaaaTCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAgg                                   <  4:66482/149‑1 (MQ=255)
                                                                                                                     aTCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCaa                                <  3:29014/149‑1 (MQ=255)
                                                                                                                      tCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCaaa                               <  6:129366/149‑1 (MQ=255)
                                                                                                                        gTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGa                             <  1:53395/149‑1 (MQ=255)
                                                                                                                           aCCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACa                           >  2:95825/1‑148 (MQ=255)
                                                                                                                           aCCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACa                           >  6:58428/1‑148 (MQ=255)
                                                                                                                              gAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAACACAAAtc                       >  2:102555/1‑149 (MQ=255)
                                                                                                                                   aaCGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGc                   <  4:94230/148‑1 (MQ=255)
                                                                                                                                    aCGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCt                 <  5:137757/149‑1 (MQ=255)
                                                                                                                                      gAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGa               <  5:10221/149‑1 (MQ=255)
                                                                                                                                       aaTAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGAc              <  4:179333/149‑1 (MQ=255)
                                                                                                                                       aaTAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGAc              <  2:173721/149‑1 (MQ=255)
                                                                                                                                         tAAAGGGAGCTTTCGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGAc                                                      >  7:175479/1‑107 (MQ=255)
                                                                                                                                         tAAAGGGAGCTTTCGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGAc                                                      <  8:175479/107‑1 (MQ=255)
                                                                                                                                         tAAAGGGAGCTTTCGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCTAAGACCAATCTCGCCTGACt             >  8:201419/1‑148 (MQ=255)
                                                                                                                                         tAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAAGTAAAATTCAATATATTGCAGTGGGGTGGGGATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAATTCGCGCCTGAcg             >  1:154491/1‑147 (MQ=255)
                                                                                                                                         tAAAGGGAGCGAAAGCTCCCTATAATATTGTTAGCAAAGGGTGCTTCGTTCATTAATTAAAAATAATTAAAATTCAATATATTGCAGAGGGGTGGGCATGATAAGACGCGCCAGCATCGCATCAGCCAAAGACAACTCTCGCCTGACt             >  8:79685/1‑148 (MQ=255)
                                                                                                                                                 gCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTGCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACcttctt    <  1:75057/149‑1 (MQ=255)
                                                                                                                                                   gAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTtcttct  <  8:41679/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCT  >  NC_000913/4488277‑4488572

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