Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,488,424 G→T 10.1% 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,4240GT10.1% 337.0 / 21.8 99intergenic (+24/+137)lptG/yjgRlipopolysaccharide export ABC permease/DUF853 family protein with NTPase fold
Reads supporting (aligned to +/- strand):  ref base G (44/45);  new base T (3/7);  total (47/52)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

ACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT  >  NC_000913/4488280‑4488570
                                                                                                                                                |                                                                                                                                                  
aCTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAg                                                                                                                                                <  1:615013/149‑1 (MQ=255)
   ggACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTc                                                                                                                                             <  2:669020/149‑1 (MQ=255)
    gACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTcc                                                                                                                                            >  1:304724/1‑149 (MQ=255)
      ccAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCt                                                                                                                                          >  1:267462/1‑149 (MQ=255)
       cAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCtt                                                                                                                                         >  2:163924/1‑149 (MQ=255)
       cAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCtt                                                                                                                                         >  1:685006/1‑149 (MQ=255)
       cAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCtt                                                                                                                                         >  2:236852/1‑149 (MQ=255)
       cAGATCTGCGGCCCGCTGACGTTGGTTTATAGCAGCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGTGAAAGCTCCCtt                                                                                                                                         <  3:51211/149‑1 (MQ=255)
         gATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCTTTCGCTCCCTtta                                                                                                                                       <  4:458135/149‑1 (MQ=255)
         gATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTtta                                                                                                                                       >  2:396991/1‑149 (MQ=255)
         gATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAACTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTtta                                                                                                                                       <  1:439386/149‑1 (MQ=255)
          aTCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTttat                                                                                                                                      <  2:182447/149‑1 (MQ=255)
           tCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCTTTCGCTCCCTttatt                                                                                                                                     <  3:439382/149‑1 (MQ=255)
            cTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCTTTCGCTCCCTttattc                                                                                                                                    <  1:262520/149‑2 (MQ=255)
            cTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCTTTCGCTCCCTttattc                                                                                                                                    <  2:349737/149‑2 (MQ=255)
            cTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCTTTCGCTCCCTttattc                                                                                                                                    <  3:245670/149‑2 (MQ=255)
             ttCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTtattat                                                                                                                                   >  1:563830/1‑149 (MQ=255)
              tCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTattatt                                                                                                                                  <  1:236852/149‑1 (MQ=255)
                 gCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCTTTCGCTCCCttt                                                                                                                                        >  2:519238/1‑140 (MQ=255)
                 gCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCTATCGCTCCCttt                                                                                                                                        <  1:519238/140‑1 (MQ=255)
                  cccGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTa                                                                                                                              <  1:147370/149‑1 (MQ=255)
                  cccGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTa                                                                                                                              <  2:491202/149‑1 (MQ=255)
                           gTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAgtgt                                                                                                                     <  1:409222/149‑1 (MQ=255)
                           gTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAgtgt                                                                                                                     >  1:330650/1‑149 (MQ=255)
                            ttGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGtgtg                                                                                                                    >  3:313724/1‑149 (MQ=255)
                              ggTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCt                                                                                                                  >  3:21112/1‑149 (MQ=255)
                              ggTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCt                                                                                                                  <  4:123974/149‑1 (MQ=255)
                                tttATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTc                                                                                                                >  4:653174/1‑149 (MQ=255)
                                  tATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGt                                                                                                              <  2:679834/149‑1 (MQ=255)
                                   aTGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTGCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGtt                                                                                                             <  2:38861/149‑1 (MQ=255)
                                    tGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCTTTCGCTCCCTttattcgtt                                                                                                                                 <  3:114922/128‑5 (MQ=39)
                                    tGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCTTTCGCTCCCTttattcgtt                                                                                                                                 >  4:114922/1‑124 (MQ=38)
                                      gCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTattatt                                                                                                                                  <  3:317649/125‑1 (MQ=255)
                                      gCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTattatt                                                                                                                                  >  4:317649/1‑125 (MQ=255)
                                      gCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCAt                                                                                                          >  3:533172/1‑149 (MQ=255)
                                       cATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCAtt                                                                                                         <  3:556511/149‑1 (MQ=255)
                                        aTCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTc                                                                                                        >  1:635082/1‑149 (MQ=255)
                                         tCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTcc                                                                                                       >  1:568435/1‑149 (MQ=255)
                                         tCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTcc                                                                                                       >  3:648670/1‑149 (MQ=255)
                                             gccgATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGaa                                                                                                   >  3:274039/1‑149 (MQ=255)
                                                   catcGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAcataa                                                                                             <  2:635082/149‑1 (MQ=255)
                                                     tcGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAAtt                                                                                           >  3:194744/1‑149 (MQ=255)
                                                              tGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCaa                                                                                  <  2:176450/149‑1 (MQ=255)
                                                              tGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCaa                                                                                  <  4:153377/149‑1 (MQ=255)
                                                                ttGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAata                                                                                >  1:693869/1‑149 (MQ=255)
                                                                  gCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAatata                                                                              <  2:304724/149‑1 (MQ=255)
                                                                    cAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATAtt                                                                            >  1:524828/1‑149 (MQ=255)
                                                                      aGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGc                                                                          >  1:76626/1‑149 (MQ=255)
                                                                           cAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAgtgg                                                                     <  3:666554/149‑1 (MQ=255)
                                                                            aGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAgtggg                                                                    <  2:582492/149‑1 (MQ=255)
                                                                             gCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAgtgggg                                                                   >  3:288626/1‑149 (MQ=255)
                                                                              cttcttcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAgtggggt                                                                  <  2:479696/149‑1 (MQ=255)
                                                                               ttcttcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAgtggggtg                                                                 >  2:304477/1‑149 (MQ=255)
                                                                                tcttcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAgtggggtgg                                                                <  2:330650/149‑1 (MQ=255)
                                                                                tcttcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAgtggggtgg                                                                >  1:184389/1‑149 (MQ=255)
                                                                                  ttcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTa                                                                                                                              <  3:564939/85‑1 (MQ=255)
                                                                                  ttcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTa                                                                                                                              >  4:564939/1‑85 (MQ=255)
                                                                                  ttcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGc                                                              >  3:515219/1‑149 (MQ=255)
                                                                                  ttcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGc                                                              >  4:253000/1‑149 (MQ=255)
                                                                                  ttcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGc                                                              >  2:516112/1‑149 (MQ=255)
                                                                                    cttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCAt                                                            <  2:267462/149‑1 (MQ=255)
                                                                                     ttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATg                                                           >  1:190893/1‑149 (MQ=255)
                                                                                        aTCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATa                                                        >  2:109753/1‑149 (MQ=255)
                                                                                           aGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGa                                                     <  2:649065/149‑1 (MQ=255)
                                                                                              cTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGAcgc                                                  <  1:121996/149‑1 (MQ=255)
                                                                                                    cTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGc                                            >  4:343292/1‑149 (MQ=255)
                                                                                                     tGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGAc                                                    <  4:370961/140‑1 (MQ=255)
                                                                                                     tGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGAc                                                    >  3:370961/1‑140 (MQ=255)
                                                                                                       ttAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATAtt                                                                            >  1:525755/1‑114 (MQ=255)
                                                                                                       ttAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATAtt                                                                            <  2:525755/114‑1 (MQ=255)
                                                                                                       ttAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTc                                         <  3:219864/149‑1 (MQ=255)
                                                                                                       ttAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTc                                         <  2:76626/149‑1 (MQ=255)
                                                                                                        tAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCg                                        <  3:570210/149‑1 (MQ=255)
                                                                                                           tGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCAt                                     <  4:158857/149‑1 (MQ=255)
                                                                                                            gagaAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATc                                    >  2:416342/1‑149 (MQ=255)
                                                                                                            gagaAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATc                                    >  3:44332/1‑149 (MQ=255)
                                                                                                             agaAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCa                                   <  2:524828/149‑1 (MQ=255)
                                                                                                              gaAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAg                                  <  4:313724/149‑1 (MQ=255)
                                                                                                               aaaaTCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAgg                                 <  2:243743/149‑1 (MQ=255)
                                                                                                               aaaaTCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAg                                  <  2:58049/148‑1 (MQ=255)
                                                                                                               aaaaTCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAg                                  >  1:58049/1‑148 (MQ=255)
                                                                                                                aaaTCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGc                                <  4:160745/149‑1 (MQ=255)
                                                                                                                   tCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCaaa                             <  4:515219/149‑1 (MQ=255)
                                                                                                                        aCCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGTTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACaa                        >  4:514795/1‑149 (MQ=255)
                                                                                                                           gAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAAtc                     >  4:41619/1‑149 (MQ=255)
                                                                                                                              tGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCg                  <  2:693869/149‑1 (MQ=255)
                                                                                                                               gAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGc                 <  4:519936/149‑1 (MQ=255)
                                                                                                                               gAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGc                 >  2:381445/1‑149 (MQ=255)
                                                                                                                                aaCGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGcc                >  4:44936/1‑149 (MQ=255)
                                                                                                                                    aaTAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGAc            <  4:194744/149‑1 (MQ=255)
                                                                                                                                      tAAAGGGAGCTTTCGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACtt          >  1:62755/1‑149 (MQ=255)
                                                                                                                                       aaaGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGAc            >  2:520595/1‑146 (MQ=255)
                                                                                                                                       aaaGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGAc            <  1:520595/146‑1 (MQ=255)
                                                                                                                                         aGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTAcc       >  2:550708/1‑149 (MQ=255)
                                                                                                                                         aGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTAcc       >  1:96965/1‑149 (MQ=255)
                                                                                                                                           ggAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACctt     >  1:434281/1‑149 (MQ=255)
                                                                                                                                            gAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACcttc    <  3:69496/149‑1 (MQ=255)
                                                                                                                                              gCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACcttctt  <  3:41619/149‑1 (MQ=255)
                                                                                                                                              gCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACcttctt  <  2:96965/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
ACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTT  >  NC_000913/4488280‑4488570

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: