Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,488,427 A→C 7.8% intergenic (+27/+134) 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,4270AC 7.8% 320.2 / 16.2 116intergenic (+27/+134)lptG/yjgRlipopolysaccharide export ABC permease/DUF853 family protein with NTPase fold
Reads supporting (aligned to +/- strand):  ref base A (51/56);  new base C (5/4);  total (56/60)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.26e-01

Only 100 of 117 total aligned reads displayed.

CTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCC  >  NC_000913/4488281‑4488575
                                                                                                                                                  |                                                                                                                                                    
cTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCt                                                                                                                                                   >  2:925201/1‑150 (MQ=255)
 tGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATAAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTc                                                                                                                                                  <  2:1556575/150‑1 (MQ=255)
  ggACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTcc                                                                                                                                                 >  1:585900/1‑150 (MQ=255)
    aCCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATTAACGAATAAAGGCAGCGAAAGCTCCCt                                                                                                                                               >  2:17381/1‑150 (MQ=255)
     ccAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCtt                                                                                                                                              >  2:1137899/1‑150 (MQ=255)
      cAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCTTTCGCTCCCttt                                                                                                                                             <  2:50491/150‑1 (MQ=255)
      cAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCTTTCGCTCCCttt                                                                                                                                             <  2:571439/150‑1 (MQ=255)
       aGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGATAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTtta                                                                                                                                            >  1:1008798/1‑150 (MQ=255)
         aTCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCTTTCGCTCCCTttatt                                                                                                                                          <  1:1647420/150‑1 (MQ=255)
         aTCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCTTTCGCTCCCTttatt                                                                                                                                          <  1:1541486/150‑1 (MQ=255)
          tCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTttatta                                                                                                                                         <  1:540447/150‑1 (MQ=255)
             tCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTg                                                                                                                                      <  2:1370137/150‑1 (MQ=255)
              cGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGt                                                                                                                                     >  1:232024/1‑150 (MQ=255)
               ggCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTTTTGtt                                                                                                                                    >  1:845823/1‑150 (MQ=255)
               ggCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGtt                                                                                                                                    <  2:1597981/150‑1 (MQ=255)
               ggCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTAAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGtt                                                                                                                                    <  2:190143/150‑1 (MQ=255)
                gCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGGGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTa                                                                                                                                   >  1:98451/1‑150 (MQ=255)
                      tGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGAACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAg                                                                                                                             <  1:197418/150‑1 (MQ=255)
                       gACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAgt                                                                                                                            <  1:1079552/150‑1 (MQ=255)
                       gACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAgt                                                                                                                            <  1:1680421/150‑1 (MQ=255)
                           ttGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGc                                                                                                                        <  2:1616893/150‑1 (MQ=255)
                                 tATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCtt                                                                                                                      <  2:1001737/146‑1 (MQ=255)
                                 tATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCtt                                                                                                                      >  1:1001737/1‑146 (MQ=255)
                                 tATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGtt                                                                                                                  <  1:220914/150‑1 (MQ=255)
                                   tGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCa                                                                                                                <  2:1209229/150‑1 (MQ=255)
                                         cccGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTg                                                                                                          <  2:1061242/150‑1 (MQ=255)
                                         cccGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTg                                                                                                          <  2:1529518/150‑1 (MQ=255)
                                             ccgATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGaaaa                                                                                                      >  1:391942/1‑150 (MQ=255)
                                                  gatcGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAaataat                                                                                                 <  1:1722585/149‑1 (MQ=255)
                                                  catcGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAACGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAaataat                                                                                                 >  2:931466/1‑150 (MQ=255)
                                                      ggCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTaaa                                                                                             >  1:860178/1‑150 (MQ=255)
                                                       gcgcACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGTACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTaaaa                                                                                            <  2:731746/150‑1 (MQ=255)
                                                         gcACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTTTGCTTCTTTCATTCCTGAAAAATAATTAAAAtt                                                                                          >  2:1540940/1‑150 (MQ=255)
                                                           aCTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCa                                                                                        <  1:1042126/150‑1 (MQ=255)
                                                            cTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCaa                                                                                       <  1:1583168/150‑1 (MQ=255)
                                                            atGTTGCCAAGCGCCATCTTCTTCTTAATCAGCCTGTGTCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCaa                                                                                       <  1:1661622/149‑1 (MQ=255)
                                                             tGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAat                                                                                      <  1:384988/150‑1 (MQ=255)
                                                             tGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAat                                                                                      <  2:437428/150‑1 (MQ=255)
                                                              gTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAata                                                                                     >  2:268058/1‑150 (MQ=255)
                                                                 gCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAatatat                                                                                  >  2:357153/1‑150 (MQ=255)
                                                                    aaGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAATCTCCCTTTATTATTGTTAGCAAATTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGc                                                                               <  1:1227858/150‑1 (MQ=255)
                                                                    aaGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGc                                                                               <  1:1476825/150‑1 (MQ=255)
                                                                    aaGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGc                                                                               <  1:995492/150‑1 (MQ=255)
                                                                        gcCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAgtg                                                                           <  1:465090/150‑1 (MQ=255)
                                                                              ttcttcttAATCAGCCTGTGGATGTTAATGAGAAAATCGTAACCGCATGAACGAATCAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAgtggggtgg                                                                     <  1:1089519/150‑1 (MQ=255)
                                                                               tcttcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAgtggggtggg                                                                    >  2:576905/1‑150 (MQ=255)
                                                                               tcttcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAgtggggtggg                                                                    <  2:434768/150‑1 (MQ=255)
                                                                                cttcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTAAATATATTGCAGTGGGGTGGGc                                                                   >  2:1542485/1‑150 (MQ=255)
                                                                                   cttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATg                                                                >  1:748531/1‑150 (MQ=255)
                                                                                    ttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGa                                                               <  2:768480/150‑1 (MQ=255)
                                                                                    ttAATCAGCCTGGTGCTGTTAATGAGAAAATCTTATCCGAATGAACGAATAAAGGTATCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCAATCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTTTGCATGa                                                               <  2:1234095/150‑1 (MQ=255)
                                                                                     tAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTACTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGTGCATGAt                                                              >  1:698622/1‑150 (MQ=255)
                                                                                     tAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCAATCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATgtt                                                              >  2:647902/1‑148 (MQ=255)
                                                                                      aaTCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATa                                                             <  1:1308760/150‑1 (MQ=255)
                                                                                      aaTCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATa                                                             <  1:1498658/150‑1 (MQ=255)
                                                                                      aaTCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATa                                                             >  1:1466246/1‑150 (MQ=255)
                                                                                       aTCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCAGTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATaa                                                            >  1:684396/1‑150 (MQ=255)
                                                                                        tCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATaaa                                                           >  1:622002/1‑149 (MQ=255)
                                                                                           gCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGAcg                                                        >  1:1551786/1‑150 (MQ=255)
                                                                                            ccTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAATAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGAcgc                                                       >  2:234262/1‑150 (MQ=255)
                                                                                            ccTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGAcgc                                                       <  1:1282119/150‑1 (MQ=255)
                                                                                              tgtgGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACgcgc                                                     <  2:1375622/150‑1 (MQ=255)
                                                                                                   cTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCg                                                >  2:1359829/1‑150 (MQ=255)
                                                                                                      ttAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCg                                             >  2:795620/1‑150 (MQ=255)
                                                                                                      ttAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCg                                             >  2:788782/1‑150 (MQ=255)
                                                                                                        aaTGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCa                                           <  1:1492685/150‑1 (MQ=255)
                                                                                                        aaTGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCa                                           >  2:1562483/1‑150 (MQ=255)
                                                                                                         aTGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCAt                                          >  2:130350/1‑150 (MQ=255)
                                                                                                          tGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATc                                         >  1:1259468/1‑150 (MQ=255)
                                                                                                          ggagaTAATCGTAACCGAATGAACGAATAAAGGGAGCGATAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAATATAATTAAAATTCAATATATTTCAGTGGGGTGGGCATGATAAGACGAGACAGCGTCGCATc                                         >  1:1259470/2‑150 (MQ=255)
                                                                                                             taaaaTCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAgg                                      <  1:1539699/149‑1 (MQ=255)
                                                                                                             gaAACTCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGGGCTGCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGGGGGCATGATAAGACGCGCCAGCGTCGCATCAgg                                      <  1:2978/150‑1 (MQ=255)
                                                                                                             gaAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAgg                                      >  2:385266/1‑150 (MQ=255)
                                                                                                              aaaaTCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGc                                     <  1:998603/150‑1 (MQ=255)
                                                                                                               aaaTCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCa                                    >  2:983905/1‑150 (MQ=255)
                                                                                                               aaaTCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCa                                    <  1:255142/150‑1 (MQ=255)
                                                                                                                  tCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAg                                 <  2:585900/150‑1 (MQ=255)
                                                                                                                    gTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGAc                               >  1:1245415/1‑150 (MQ=255)
                                                                                                                     tAACCGAATGAACGAATAAAGGGAGCGAAAGCTACCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACa                              <  2:1375874/150‑1 (MQ=255)
                                                                                                                      aaCCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTTAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCATGAAAAGACaa                             >  2:542307/1‑150 (MQ=255)
                                                                                                                      aaCCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACaa                             >  2:540029/1‑150 (MQ=255)
                                                                                                                        ccGAATGAACGAATAAAGGGAGCGAAAGCGCCCTTTATTATTGTTAGCAAAGTGGGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGGGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAAt                           <  2:1078311/150‑1 (MQ=255)
                                                                                                                            aTGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCg                       <  2:860178/150‑1 (MQ=255)
                                                                                                                              gAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGcc                     >  1:524278/1‑150 (MQ=255)
                                                                                                                                aataaTAAAGGGAGCTTTCGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTg                   >  2:252868/4‑150 (MQ=255)
                                                                                                                                aataaTAAAGGGAGCTTTCGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTACTGAAAAATAATTTAAATTCAATATATTGCAGTGGGGTGGGGATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTg                   >  2:1127907/4‑150 (MQ=255)
                                                                                                                                aCGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTg                   >  2:591782/1‑150 (MQ=255)
                                                                                                                                    aTAAAGTGATCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGTGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACtt               <  1:919276/150‑1 (MQ=255)
                                                                                                                                    aTAAAGGGAGCTTTCGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACtt               >  2:1491670/1‑150 (MQ=255)
                                                                                                                                     tAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGAATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTa              >  1:1595555/1‑150 (MQ=255)
                                                                                                                                      aaaGGGAGCTTTCGCTCCCTTTATTATTGTTAGCAAAGTGTGCGTCATTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGTCAAAGCCAAATCTCGCCTGACttcc             >  2:140089/1‑148 (MQ=255)
                                                                                                                                        aGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAAATAAAAATCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACct           >  1:934874/1‑150 (MQ=255)
                                                                                                                                         gggAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACctt          >  2:631569/1‑150 (MQ=255)
                                                                                                                                          ggAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACcttc         >  2:982195/1‑150 (MQ=255)
                                                                                                                                             gCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACATATCTCGCCTGACTTAGCttcttc      >  1:1120518/1‑150 (MQ=255)
                                                                                                                                             gCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCttcttc      <  1:788782/150‑1 (MQ=255)
                                                                                                                                             gCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCttcttc      <  1:1571227/150‑1 (MQ=255)
                                                                                                                                             gCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCttcttc      <  1:795620/150‑1 (MQ=255)
                                                                                                                                             gCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCttcttc      <  2:1131213/150‑1 (MQ=255)
                                                                                                                                                 aaGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTccc  <  1:5234/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCC  >  NC_000913/4488281‑4488575

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: