Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 107,554 G→T 13.1% intergenic (+80/‑151) 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,5540GT13.1% 528.2 / 16.1 168intergenic (+80/‑151)lpxC/secMUDP‑3‑O‑acyl N‑acetylglucosamine deacetylase/regulator of secA translation
Reads supporting (aligned to +/- strand):  ref base G (15/126);  new base T (25/2);  total (40/128)
Fisher's exact test for biased strand distribution p-value = 2.70e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 169 total aligned reads displayed.

TATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAAT  >  NC_000913/107406‑107701
                                                                                                                                                      |                                                                                                                                                     
 tggggACCTTCCTGGCCGACGCCGAACTGCGGTTGGCCTTCAACGCGCCTTCTGGTGTACGGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCGGGCACTCTCTCGGGCCAGGTGAACCTGTCG‑T‑TTTTTTTTGaa                                                                                                                                                     <  2:2415923/146‑1 (MQ=255)
     gACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TGTTTTTTTTTAAttt                                                                                                                                                  >  2:3490904/1‑150 (MQ=255)
     gACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGTAtttt                                                                                                                                                 >  2:2673870/1‑150 (MQ=255)
      aCCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTTAttttt                                                                                                                                                >  1:1024312/1‑149 (MQ=255)
        cTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑GTTTttttttatttttta                                                                                                                                              >  2:443766/1‑150 (MQ=255)
             aGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑Gc                                                                                                                                        <  1:4111870/150‑1 (MQ=255)
             aGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TGTTTTTTTTTTTTTATaattgg                                                                                                                                         >  2:445350/1‑146 (MQ=255)
                  gacgCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTT‑TTTTTTTTGAATTTTATAAGA‑GCtat‑a                                                                                                                                   <  2:1255810/150‑1 (MQ=255)
                      ctgAACTGCCGTTGGTCTTCAATGCGCGTTCAGCTGTACTTTAATATCGTGATTTATAGTGTCGTTTAATATTAGACTGGATTTTCTTTAACTCTCTCCGGCCAGGTGAACATGTCG‑T‑TTTTTTTTGAATTTTATATGA‑GCTTT‑taaaaa                                                                                                                              <  1:3760423/148‑1 (MQ=255)
                      cAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑aaaaaa                                                                                                                              <  2:3350747/150‑1 (MQ=255)
                      cAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑aaaaaa                                                                                                                              <  2:3350653/150‑1 (MQ=255)
                        gAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTGCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACg                                                                                                                            <  2:4033696/150‑1 (MQ=255)
                         aaCTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGTATTTTTTAAGA‑GGTTT‑AAAAAAAgg                                                                                                                           >  2:3346396/1‑150 (MQ=255)
                          aCTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAATTTCGACTGGCAAATCTGGCACTCTCTCCGGCCTGGTTAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGt                                                                                                                          <  1:1392633/150‑1 (MQ=255)
                           cTGCCGTTGGCCTACAACGCGCCTTCAGCTGTACTGGCATACCGAAATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCACCGGCCTGTTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTg                                                                                                                         <  1:2449900/150‑1 (MQ=255)
                            tGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑GTTTTTTTTGAGTTTATAAGATGGAAA‑AAAAAACgggg                                                                                                                         >  2:2140346/1‑148 (MQ=255)
                              ccGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑GTTTTTTTGTATTTTTTAAGA‑GCGAA‑AAAAAAACGTGCgc                                                                                                                      >  2:1697048/1‑149 (MQ=255)
                                 ttGGCCTTCAAAGCGCCTTCAGCTGTACTGGAATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACg                                                                                                                   <  2:3397651/150‑1 (MQ=255)
                                  tGGCCTTCAATGCGCCTTAAGCTGTACTGGCATAACGACTTTTATACTGTCGTATAATTTTCGACTGGCAATTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGc                                                                                                                  <  1:3183833/150‑1 (MQ=255)
                                  tGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGc                                                                                                                  <  2:487666/150‑1 (MQ=255)
                                  tGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGc                                                                                                                  <  1:3374893/150‑1 (MQ=255)
                                   ggCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCt                                                                                                                 <  1:2407497/150‑1 (MQ=255)
                                   ggCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCt                                                                                                                 <  1:1290275/150‑1 (MQ=255)
                                      cTTCAAAGCTCCTTCAGCTGTATTGGCATAACGACATTTATACTTTCGTATAAATTTCGACGGTCAAATCTGGCACTCTCTCGGGCCTGGTGAACCAGTCG‑T‑TTTTTTTTTAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGtt                                                                                                              <  1:2804205/150‑1 (MQ=255)
                                      cTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGtt                                                                                                              <  2:1404357/150‑1 (MQ=255)
                                       ttCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTTATTTTTTAATA‑TCTAT‑AAAAAACGGTTGGAAAGCGGttt                                                                                                             >  2:2742156/1‑150 (MQ=255)
                                       ttCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTTATTTTATAATA‑TCGAT‑AAAAAAAGGTGCGAACGCTGttt                                                                                                             >  2:1353187/1‑150 (MQ=255)
                                         cAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTc                                                                                                           <  2:4080930/150‑1 (MQ=255)
                                         cAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTc                                                                                                           <  2:3806858/150‑1 (MQ=255)
                                           aaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCtt                                                                                                         <  2:1241774/150‑1 (MQ=255)
                                           aaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCtt                                                                                                         <  2:3793634/150‑1 (MQ=255)
                                            aGCGCCTTCAGCTGTACTGGAATAACGACATTTATACTGTCTTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCTGTTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTa                                                                                                        <  2:3649060/150‑1 (MQ=255)
                                             gcgcCTTCTGCTGTCCTGTCATACCGACTTTTATACTGTCGTATAAAATTCGACTGGAAAATCTGGCACTCTCTCCGGCCAGGTGAACCATTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAATCGGTGCGAACGCTGTTTTCTTaa                                                                                                       <  1:2691689/150‑1 (MQ=255)
                                             gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTaa                                                                                                       <  1:1729851/150‑1 (MQ=255)
                                             gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTaa                                                                                                       <  2:51657/150‑1 (MQ=255)
                                             gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTaa                                                                                                       <  1:70014/150‑1 (MQ=255)
                                                ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCa                                                                                                    <  1:813978/150‑1 (MQ=255)
                                                     atcTGTTCTGTCATACCGACATTTTTACTGTCGTTTAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCTTGTTAAACATTCG‑T‑TTTTTTTTGAATTTTATACGA‑GCTAT‑TATAATCGGTGCGAACGCTGTTTTCTTAAGCACtttt                                                                                               <  1:1027442/148‑1 (MQ=255)
                                                     aGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTTATTTTATAATA‑TATAT‑AAAAAAAAGTGAGAAACCGGTTTTATTAAGCaaattt                                                                                               >  2:1969403/1‑145 (MQ=255)
                                                       cTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTTAATTTATAATA‑GCTAA‑AAAAAAACGTGCGAACCCGGTTTTATTAAGAACatttac                                                                                             >  2:312328/1‑148 (MQ=255)
                                                       cTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTcc                                                                                             <  2:2554977/150‑1 (MQ=255)
                                                       cTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTcc                                                                                             <  1:2196528/150‑1 (MQ=255)
                                                        tGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                                                                                            <  1:2134765/150‑1 (MQ=255)
                                                        tGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                                                                                            <  2:748948/150‑1 (MQ=255)
                                                         gTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                                                                                           <  1:901578/150‑1 (MQ=255)
                                                         gTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                                                                                           <  2:4053342/150‑1 (MQ=255)
                                                               gCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TGTTTTTTTTTAATTTT‑TAATA‑GCAAA‑AAAAAAAAGGGCGAAAGCGGTTTTGttaagcacaattaagcacaaaa                                                                                     >  2:3226175/1‑148 (MQ=255)
                                                                  tAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTTATTTTTTAAAA‑TCAAT‑AAAAAAACGTTCGGAACCGGTTTTATTAAGAACATTTGAGCACAACCTAt                                                                                  >  2:1143336/1‑150 (MQ=255)
                                                                  tAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑GTTTTTTTTTATTTTTTAAGA‑GCAAT‑AAAAAAAAGTGGGAACGCGGTTTTGTTAATCACAATTTCGCCCaaagtat                                                                                  >  2:1143422/1‑145 (MQ=255)
                                                                   aaCGACATTTATACTGTCGTATAAAATACGACTGGCAAATCTGCCACACACTCCGCCCAGGTGAACCAGTCG‑A‑TTTTTTTTGAAATTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATc                                                                                 <  1:3040187/150‑1 (MQ=255)
                                                                    aCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCt                                                                                <  1:1709456/150‑1 (MQ=255)
                                                                    aCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCt                                                                                <  1:3593384/150‑1 (MQ=255)
                                                                          tttATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTc                                                                          <  1:3517147/150‑1 (MQ=255)
                                                                            tataCTGTCGTAAAAAATTCGACTGGCAAATCTGGCACTCTCACCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGt                                                                        <  1:3568030/150‑1 (MQ=255)
                                                                             ataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTg                                                                       <  2:3188402/150‑1 (MQ=255)
                                                                                      tataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTAAATTTATAAGA‑GCTAT‑AAAAAAACGTGCGAACGCTGTTTTATTAAGCAATTTTGCGCACAACGTATAATAAATAGTTCTGTGGACt                                                              >  2:3084008/1‑150 (MQ=255)
                                                                                       aaaaaaTTCGCCTGGAAAATCTGGCACCCTCTCCGGCCTGTTGAACATTTGG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTg                                                             <  1:1978089/148‑1 (MQ=255)
                                                                                         aaaaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCa                                                           <  2:2820949/150‑1 (MQ=255)
                                                                                            aTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGc                                                        <  1:1905579/150‑1 (MQ=255)
                                                                                             ttCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCt                                                       <  2:695629/150‑1 (MQ=255)
                                                                                                 aCTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGT‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTaa                                                   <  2:4004928/150‑1 (MQ=255)
                                                                                                 aCTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTaa                                                   <  2:2631308/150‑1 (MQ=255)
                                                                                                 aCTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTaa                                                   <  1:2859933/150‑1 (MQ=255)
                                                                                                 aCTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTaa                                                   <  2:2306190/150‑1 (MQ=255)
                                                                                                  cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAAt                                                  <  1:1668751/150‑1 (MQ=255)
                                                                                                     gCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGAt                                               <  2:2520919/150‑1 (MQ=255)
                                                                                                     gCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGAt                                               <  1:3716508/150‑1 (MQ=255)
                                                                                                      cAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATa                                              <  2:956410/150‑1 (MQ=255)
                                                                                                      cAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATa                                              <  2:1046030/150‑1 (MQ=255)
                                                                                                       aaaTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTTATTTTATTAGA‑GCGATGAAAAAACGG‑GGGCAAGCGGTTTTGTTAAGCACATTTACGCACAACATATAATGAATCGTGCTGTGGCCTGCACGCTTTAATGATaa                                             >  2:235413/1‑150 (MQ=255)
                                                                                                            tGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑GTTTTTTTTTATTTTTTAAGA‑TCTAA‑AAAAAAAGGTGGGTAAGCAGTTTTATTAAGCAAATTTACCCACAACATATATTAAATAGTGCTTTGGTCAGCAGGCTTTAAAGATaaaaatt                                        >  2:3279483/1‑145 (MQ=255)
                                                                                                             ggCACTCTCTCCGGCCAGGTTAACCTGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTtg                                       <  1:443766/150‑1 (MQ=255)
                                                                                                              gCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTTTATTTTATAATA‑GCTAT‑AAAAAAAGGTGCGAAAGCTGTTTTATTAATCACATTTAAGCACAAAATATATTAATTAGTTCTTTGGTCAGCAGGCTTTAAATATAATAttttt                                      >  2:265401/1‑148 (MQ=255)
                                                                                                              gCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTtgt                                      <  1:2753748/150‑1 (MQ=255)
                                                                                                                 ctctctCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTgcg                                   <  1:1411073/150‑1 (MQ=255)
                                                                                                                 ctctctCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTgcg                                   <  2:1204547/150‑1 (MQ=255)
                                                                                                                  tctctcCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGAGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGACAAGATTTGTgcgc                                  <  2:3108664/150‑1 (MQ=255)
                                                                                                                      tcCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTaaa                              <  2:1883078/150‑1 (MQ=255)
                                                                                                                             aGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGttt                       <  1:744355/150‑1 (MQ=255)
                                                                                                                                tGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGaa                    <  1:378831/150‑1 (MQ=255)
                                                                                                                                tGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGaa                    <  1:3061364/150‑1 (MQ=255)
                                                                                                                                tGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGaa                    <  2:2943123/150‑1 (MQ=255)
                                                                                                                                ggAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACACCTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGaa                    <  2:2843611/149‑1 (MQ=255)
                                                                                                                                   aCCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAatat                 <  1:101012/150‑1 (MQ=255)
                                                                                                                                      aGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGAt              <  2:3620331/150‑1 (MQ=255)
                                                                                                                                       gTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATc             <  2:3512832/150‑1 (MQ=255)
                                                                                                                                         cG‑T‑TTTTTTTTTAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAATGCAGTCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCgg           <  1:3931451/150‑1 (MQ=255)
                                                                                                                                         cG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCgg           <  1:251923/150‑1 (MQ=255)
                                                                                                                                         cG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCgg           <  2:815555/150‑1 (MQ=255)
                                                                                                                                          g‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCggg          >  1:610869/1‑150 (MQ=255)
                                                                                                                                          g‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCggg          >  1:610843/1‑150 (MQ=255)
                                                                                                                                            t‑ttttttttGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGGCTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATAgggg         >  1:2685653/1‑149 (MQ=255)
                                                                                                                                            t‑ttttttttGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGa         <  2:1528437/150‑1 (MQ=255)
                                                                                                                                               tttttttGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATg       <  1:1790425/150‑1 (MQ=255)
                                                                                                                                                ttttttGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATgg      <  2:2058917/150‑1 (MQ=255)
                                                                                                                                                ttttttGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATgg      <  2:3545582/150‑1 (MQ=255)
                                                                                                                                                ttttttGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCAAAAATTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATgg      <  2:943570/150‑1 (MQ=255)
                                                                                                                                                 tttttGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGTATATGATCGGGATGGc     >  1:1262216/1‑150 (MQ=255)
                                                                                                                                                  ttttGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCa    <  1:394440/150‑1 (MQ=255)
                                                                                                                                                   tttGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCaa   >  1:1411350/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TATGTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑T‑TTTTTTTTGAATTTTATAAGA‑GCTAT‑AAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAAT  >  NC_000913/107406‑107701

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: ‑