Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 107,578 C→A 13.2% intergenic (+104/‑127) lpxC → / → secM UDP‑3‑O‑acyl N‑acetylglucosamine deacetylase/regulator of secA translation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913107,5780CA13.2% 711.5 / 40.3 223intergenic (+104/‑127)lpxC/secMUDP‑3‑O‑acyl N‑acetylglucosamine deacetylase/regulator of secA translation
Reads supporting (aligned to +/- strand):  ref base C (49/139);  new base A (32/2);  total (82/141)
Fisher's exact test for biased strand distribution p-value = 2.59e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 223 total aligned reads displayed.

AACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCT  >  NC_000913/107431‑107726
                                                                                                                                                      |                                                                                                                                                    
aaCTGCCTTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACgg                                                                                                                                                    <  2:440913/150‑1 (MQ=255)
 aCTGCCGTTGGCCTTCAACGCGCCTTCAGCTGTACTGGCATAACGAAATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCT‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGt                                                                                                                                                   <  1:2048359/150‑1 (MQ=255)
        ttGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTAATTTTTTTAAATCAA‑TAAAAAAAAGTGCGaaac                                                                                                                                            >  2:1945445/1‑148 (MQ=255)
          ggCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTTATTTTTTAAAATGTA‑TAAAAAAAAGTTGGAACGcg                                                                                                                                          >  1:3691623/1‑149 (MQ=255)
           gCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTTATTTTTTAAAAACGA‑TAAAAAAAAGTTGGAACGcgg                                                                                                                                         >  1:4415026/1‑148 (MQ=255)
             cTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGtt                                                                                                                                       >  2:3851062/1‑150 (MQ=255)
              ttCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTT‑TTTTTTTTGAATTTT‑TAAGAGCGA‑TAAAAAAAAGTGCGAACGCGGttt                                                                                                                                      >  2:2496770/1‑150 (MQ=255)
                cAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTc                                                                                                                                    <  1:1059990/150‑1 (MQ=255)
                cAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTc                                                                                                                                    <  2:4840444/150‑1 (MQ=255)
                 aaaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCt                                                                                                                                   <  2:3938941/150‑1 (MQ=255)
                 aaaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCt                                                                                                                                   <  1:3170000/150‑1 (MQ=255)
                 aaaGCGCCTTCAGCGTTTCTGGCAATACCCAATTTATACGTGCGTAAAAAATTCGCCGGGCAAAAGCGGGCCCCCCACCGGCCAGGTTAACCCGTCGTT‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTc                                                                                                                                    <  1:2772054/150‑1 (MQ=255)
                  aaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCtt                                                                                                                                  <  2:645052/150‑1 (MQ=255)
                  aaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCtt                                                                                                                                  <  1:1477352/150‑1 (MQ=255)
                   aGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTa                                                                                                                                 <  1:3182057/150‑1 (MQ=255)
                    gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTaa                                                                                                                                <  1:179870/150‑1 (MQ=255)
                    gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTaa                                                                                                                                <  1:515063/150‑1 (MQ=255)
                    gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTaa                                                                                                                                <  1:3051225/150‑1 (MQ=255)
                    gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTaa                                                                                                                                <  1:3773144/150‑1 (MQ=255)
                     cgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAg                                                                                                                               <  2:39499/150‑1 (MQ=255)
                       ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCa                                                                                                                             <  2:162476/150‑1 (MQ=255)
                             gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TGTTTTTTTTTAATTTT‑TAAGATCTATTAAAAAACGG‑GCGCACCCCGGTTTATTTAGCACatttg                                                                                                                       >  2:3689955/1‑149 (MQ=255)
                                         tAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                                                                                                           <  2:902785/150‑1 (MQ=255)
                                          aaCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATc                                                                                                          <  2:3454007/150‑1 (MQ=255)
                                           aCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCt                                                                                                         <  1:1810833/150‑1 (MQ=255)
                                            cGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCtt                                                                                                        <  1:4122208/150‑1 (MQ=255)
                                             gACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑GTTTTTTTGGATTTTATAAAAGCGA‑TAAAAAAACGTGCGACCCGAGGTTTGTTTAGCACGTTTAAGCACAACATATAttg                                                                                                       >  2:5038327/1‑149 (MQ=255)
                                               cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACTCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCAt                                                                                                     <  2:3610244/150‑1 (MQ=255)
                                                 tttATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGTATTTTTTAACAGCTA‑TAAAAAAAGGTGCGGACGCGGTTTTGTTTAGCACTTTTTCGCACAACTTTTATTccata                                                                                                   >  2:4253512/1‑146 (MQ=255)
                                                      aCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCt                                                                                              <  2:1679467/150‑1 (MQ=255)
                                                        tGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgt                                                                                            <  2:912400/150‑1 (MQ=255)
                                                        tGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgt                                                                                            <  2:912187/150‑1 (MQ=255)
                                                        tGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgt                                                                                            <  1:1775851/150‑1 (MQ=255)
                                                        tGTAGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgt                                                                                            <  2:462720/150‑1 (MQ=255)
                                                              ataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTg                                                                                      <  1:3143353/150‑1 (MQ=255)
                                                              ataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAAAGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTg                                                                                      <  2:1507563/150‑1 (MQ=255)
                                                               taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGTATTTTATAAAATCTA‑TAAAAAAAGGTGCGAACGCTGTTTTCTTAATCACATTTCCGCACAACTTATCATAAATAGTGCTGTGGACttc                                                                                     >  2:4295858/1‑148 (MQ=255)
                                                               taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑GTTTTTTTTAAATTTATAAGATCTA‑TAAAAAACCGTGCGAACCCAGGTTTCGTAAGCACATTTACGCACAACATATCTTCAATCGTGCTGTgggcagc                                                                                     >  2:2585020/1‑145 (MQ=255)
                                                               taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑GTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGGACGCGGTTTTCTTAAGCACGTTTACGCGCAACGTATATTCATTAGTGCTggggtcagc                                                                                     >  2:5097717/1‑145 (MQ=255)
                                                                  aagtCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAgg                                                                                  <  1:2248096/147‑1 (MQ=255)
                                                                     tCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCtt                                                                               <  2:1029593/150‑1 (MQ=255)
                                                                         cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAAt                                                                           <  1:4648796/150‑1 (MQ=255)
                                                                         cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAAt                                                                           <  2:3062403/150‑1 (MQ=255)
                                                                          tGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTTATTTTATAATAGCTA‑TAAAAAAAGGTGCGGACGCGGTTTTGTTAAGCACATTTAAGGACAAATTATGATGATTGGTGGTGTGGACGGCGGGCGTTAAtt                                                                          >  2:2105394/1‑149 (MQ=255)
                                                                          tGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTAAATTTATAAAAGCTA‑AAAAAAAAGGTGCGCACGCGGGTTTGTTAAGCACATTTGCGCACAACGTATCGTAATTAGGGGGGTGGGCTGCCGGCGTTAATg                                                                          >  2:297733/1‑150 (MQ=255)
                                                                          tGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATg                                                                          <  1:4591693/150‑1 (MQ=255)
                                                                           ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGa                                                                         <  1:2824274/150‑1 (MQ=255)
                                                                            gCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTTATTTTTTAAGAGCGA‑TAAAAAAAAGTGCGAAAGCGGTTTTGTTAAGAAAATTTAAGCACAAAAAAACTTCAGAAGTTCGTTGGTCTGCAGGCGTTAAttgt                                                                        >  2:2676155/1‑147 (MQ=255)
                                                                            gCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTTATTTTATAATAGCTA‑TAAAAAAAGGTGCGAAAGCGGTTTTGTTAATCACTTTTGAGCACAAATTATATTCATTAGTGCTGTGGACTGGAGGCTTTTATGAt                                                                        >  2:3870083/1‑150 (MQ=255)
                                                                             cAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATa                                                                       <  2:3373973/150‑1 (MQ=255)
                                                                                aTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGa                                                                    <  2:2203923/150‑1 (MQ=255)
                                                                                 tCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGAt                                                                   <  2:2792661/150‑1 (MQ=255)
                                                                                 tCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGAt                                                                   <  2:4202008/150‑1 (MQ=255)
                                                                                 tCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑GTTTTTTTGTATTTTATAAGAGCTA‑TAAAAAAAGGGGCGAAACCAGTTTTCTTAAGCACATTTCCCCCCAACGTAACATCATTCGTGCTGTGGACTGCAGGCGGTAATGATAAGag                                                                   >  2:3752530/1‑149 (MQ=255)
                                                                                   tGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAAATTTATAATATCTA‑TAAAAAAACCTGCGAACGCAGTTTTATTTAACACTATTACGCACAAAAAAACTTAAATAGTGGAGTGGTCTGCAGGCTTTAATGATAAAAttt                                                                 >  2:2051164/1‑150 (MQ=255)
                                                                                       aCTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTgc                                                             <  2:4740191/150‑1 (MQ=255)
                                                                                       aCTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTgc                                                             <  1:4212089/150‑1 (MQ=255)
                                                                                         tctctcCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTTATTTTATAAGAGCTA‑AAAAAAAAGGGGCGGACCCGGGTTTCTTTAGCACATTTACCCCCAACATATCTTACATAGTGCGGTGGGCGGCAGGCTTTAAAGGTAAGATTTTTgcgc                                                           >  2:2920576/1‑150 (MQ=255)
                                                                                           tctcCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCGGTTTTCTTAAGCACTTTTGCGCACAAAGTATAGTCATTGGTGGGGTGGGCTGCGGGCGTTTATTATAAGATTTTTGCGCTa                                                         >  2:5115669/1‑150 (MQ=255)
                                                                                             tcCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTaaa                                                       <  1:3894823/150‑1 (MQ=255)
                                                                                               cGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTAATTTTTTAAGAGCTA‑TAAAAAAAGGTGCGAACCCGGTTTTCTTAATCAAATTTACGCACAAAGTATCATAAATAGTTCAGTGGGCAGCACGCTTTAATTATAAGAATTTTGCGCTAAATa                                                     >  2:3347968/1‑150 (MQ=255)
                                                                                                 gCCAGGTGAACCAGTCG‑T‑TTTTTTTTTTATTTTATAATAGCTA‑AAAAAAAAGGTGCGAACCCGTTTTTGTTAATCACATTTACGCACAAAAAAAATTAAATAGTGCGGTTGACTGCAGGCTATAAAGATAAAATTTTTTCGCAAAATaaa                                                   >  2:3226680/1‑148 (MQ=255)
                                                                                                    aGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGttt                                                <  1:2155622/150‑1 (MQ=255)
                                                                                                    aGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGttt                                                <  1:2454353/150‑1 (MQ=255)
                                                                                                    aGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTCTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGttt                                                <  1:4315792/150‑1 (MQ=255)
                                                                                                    aGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAAAGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGGTGTGTACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGttt                                                <  1:4741053/150‑1 (MQ=255)
                                                                                                       tGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGaa                                             <  1:2621834/150‑1 (MQ=255)
                                                                                                       tGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGAGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGaa                                             <  1:3549668/150‑1 (MQ=255)
                                                                                                         aaCCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAata                                           <  1:1956752/150‑1 (MQ=255)
                                                                                                          aCCAGTCGTT‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAata                                           <  1:4233996/150‑1 (MQ=255)
                                                                                                               tCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCg                                     <  1:1402737/150‑1 (MQ=255)
                                                                                                               tCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCg                                     <  1:687766/150‑1 (MQ=255)
                                                                                                               tCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCg                                     <  1:2345039/150‑1 (MQ=255)
                                                                                                                 g‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCggg                                   <  1:1741254/150‑1 (MQ=255)
                                                                                                                   t‑ttttttttGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGa                                  <  1:1177815/150‑1 (MQ=255)
                                                                                                                   t‑ttttttttGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGa                                  <  2:1152350/150‑1 (MQ=255)
                                                                                                                   t‑ttttttttGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGa                                  >  2:2297354/1‑150 (MQ=255)
                                                                                                                   t‑ttttttttGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGAACGGGa                                  <  2:345802/150‑1 (MQ=255)
                                                                                                                      tttttttGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATg                                >  1:739274/1‑150 (MQ=255)
                                                                                                                      tttttttGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATg                                >  1:1527441/1‑150 (MQ=255)
                                                                                                                      tttttttGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATg                                >  1:2680328/1‑150 (MQ=255)
                                                                                                                        tttttGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGc                              <  1:3953494/150‑1 (MQ=255)
                                                                                                                          tttGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCaa                            >  1:782914/1‑150 (MQ=255)
                                                                                                                                ttttATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGt                      <  1:3669502/150‑1 (MQ=255)
                                                                                                                                 tttATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTg                     >  2:1290102/1‑150 (MQ=255)
                                                                                                                                      aaGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTgg                >  1:101067/1‑150 (MQ=255)
                                                                                                                                      aaGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTgg                >  2:259627/1‑150 (MQ=255)
                                                                                                                                       agagCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGGTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGa               >  2:3294796/1‑150 (MQ=255)
                                                                                                                                       agagCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGa               >  1:5099198/1‑150 (MQ=255)
                                                                                                                                       agagCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGa               <  1:440036/150‑1 (MQ=255)
                                                                                                                                        gagCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGaa              >  2:4164082/1‑150 (MQ=255)
                                                                                                                                        gagCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGaa              >  2:1806622/1‑150 (MQ=255)
                                                                                                                                           cTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATAc           >  2:1917892/1‑150 (MQ=255)
                                                                                                                                            ta‑taAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACt          >  2:2310630/1‑150 (MQ=255)
                                                                                                                                               taAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGa        <  2:2410048/150‑1 (MQ=255)
                                                                                                                                                aaaaaaCGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGAc       <  2:779202/150‑1 (MQ=255)
                                                                                                                                                 aaaaaCGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGAcg      >  2:4813573/1‑150 (MQ=255)
                                                                                                                                                  aaaaCGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATTGGGATGGCAATAACg                       >  1:2580108/1‑132 (MQ=255)
                                                                                                                                                   aaaCGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGAcgcg    >  1:3290507/1‑150 (MQ=255)
                                                                                                                                                     aCGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCt  <  1:1567174/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
AACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGAGCTA‑TAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAACGTGAGTGGAATACTGACGCGCT  >  NC_000913/107431‑107726

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑