Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 107,554 G→T 10.7% 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,5540GT10.7% 827.7 / 21.9 253intergenic (+80/‑151)lpxC/secMUDP‑3‑O‑acyl N‑acetylglucosamine deacetylase/regulator of secA translation
Reads supporting (aligned to +/- strand):  ref base G (21/197);  new base T (31/2);  total (54/199)
Fisher's exact test for biased strand distribution p-value = 1.93e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 259 total aligned reads displayed.

GTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAA  >  NC_000913/107409‑107703
                                                                                                                                                 |                                                                                                                                                     
 tgACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTTAAttt                                                                                                                                                  >  2:761834/1‑150 (MQ=255)
   aCCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAACGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTta                                                                                                                                                <  1:2716455/150‑1 (MQ=255)
      ttCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGAATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTtatta                                                                                                                                             <  1:2073682/150‑3 (MQ=255)
      ttCCAGGACGACGCAGAACTGCCCTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCACTCGTTTTTTTTTGAATTTTATaa                                                                                                                                             <  2:4549039/150‑1 (MQ=255)
       tCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTGTTTTTTTTAATATTATaaa                                                                                                                                            >  2:2150511/1‑149 (MQ=255)
        ccAGGCCGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAaga                                                                                                                                           <  1:5593500/150‑1 (MQ=255)
        ccAGGACGACGCAGAACTGCCGTTGGTCTTCAAACCGCCTTCAGCTGTACAGGCATAACCACAATTATACTGACGAATAAAATTCGACTGCCAAATCTGCCACTCTCTCCGGCCAGGTGAACCAGTCGTATTTTTTTGAATTTTATAaga                                                                                                                                           <  2:1466585/150‑1 (MQ=255)
          aGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGTATTTTTTAAGGGc                                                                                                                                         >  2:873231/1‑150 (MQ=255)
          aGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTGTTTTTTtttattttattatatg                                                                                                                                         >  2:1403584/1‑144 (MQ=255)
           ggACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCt                                                                                                                                        <  1:1851479/150‑1 (MQ=255)
           ggACGACGCAGAACTGCCGTTGGCCTTCAAAGCCCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCt                                                                                                                                        <  1:4547178/150‑1 (MQ=255)
             acgacgCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGTATTTTTTaaaatcaat                                                                                                                                      >  1:4440304/1‑143 (MQ=255)
                  gCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATaaaaa                                                                                                                                 >  2:329720/1‑150 (MQ=255)
                    aGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAAc                                                                                                                               <  1:4812696/150‑1 (MQ=255)
                     gAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACg                                                                                                                              <  2:3742154/150‑1 (MQ=255)
                     gAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACAGGCAAATATGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTTTTATTTATAAAATATATaaaaaaaa                                                                                                                              >  1:1705533/1‑148 (MQ=255)
                       aCTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGt                                                                                                                            <  1:4627026/150‑1 (MQ=255)
                       aCTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGt                                                                                                                            <  2:362164/150‑1 (MQ=255)
                        cTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTg                                                                                                                           <  2:2070700/150‑1 (MQ=255)
                        cTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTg                                                                                                                           <  2:5171123/150‑1 (MQ=255)
                         tGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTGTTTTTTTTGTTTTATAAAATCTTAaaaaaaaagtgg                                                                                                                          >  1:3679695/1‑144 (MQ=255)
                           ccGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTGTTTTTTTAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                                                                                       >  2:3998070/1‑150 (MQ=255)
                             gTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAATATCTATAAAAAAAGGTGCGAAc                                                                                                                      >  2:4428923/1‑150 (MQ=255)
                                   cTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGtt                                                                                                                <  1:4270579/150‑1 (MQ=255)
                                   cTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGtt                                                                                                                <  2:3392686/150‑1 (MQ=255)
                                   cTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTGTTTTTTTTTATTTTTTAATATCTATAAAAAAAAGTGCGAAAGCAGtt                                                                                                                >  2:605843/1‑150 (MQ=255)
                                    ttCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                                                                                               <  1:3058650/150‑1 (MQ=255)
                                    ttCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTGTTTTTTTGTATTTTTTAAGTGCTTTAAAAAAAGGGGGGGACGcgggtt                                                                                                               >  2:1711807/1‑145 (MQ=255)
                                     tCAACGCGCCTTCAGCTGTACTGGCATAACGACATTTTTACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGtttt                                                                                                              <  2:2527037/150‑1 (MQ=255)
                                      cAACGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTc                                                                                                             <  2:3442834/150‑1 (MQ=255)
                                      cAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGTAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTc                                                                                                             <  2:5292287/150‑1 (MQ=255)
                                       acaGCGCCTTCATCTGTTCTGGCTTAACGACATTTTTACTTTCGTATAAACTTCGACTGGCAAATGTGGCACTCTCTCCGGCCAGGTTAACGAGTCTTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCt                                                                                                            <  1:4959285/148‑1 (MQ=255)
                                       aaaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCt                                                                                                            <  2:5781184/150‑1 (MQ=255)
                                       aaaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCt                                                                                                            <  2:3506210/150‑1 (MQ=255)
                                          gcgcCTTCCGCTGTACTGGCATAACGACATTTATACTTTGGTATAAAATTCGCCTGGCAAATCTGGCACTCTCTCCGGCCTGTTGAACCAGTGGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTaa                                                                                                         <  2:2918832/150‑1 (MQ=255)
                                          gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTTTATTTTATAATATAAATAAAAAAAAGTGCGAAAGCGGTTTTGTTaa                                                                                                         >  2:1347207/1‑150 (MQ=255)
                                          gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTaa                                                                                                         <  1:728053/150‑1 (MQ=255)
                                          gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTGTTTTTTTTTATTTTATAATAGCTATAAAAAACGGTGCGCACGCGGGTTTCTTaa                                                                                                         >  2:2423955/1‑150 (MQ=255)
                                                 cAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACttt                                                                                                  <  2:3025217/150‑1 (MQ=255)
                                                  aGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtttt                                                                                                 <  1:25704/150‑1 (MQ=255)
                                                    cTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTcc                                                                                               <  1:1934265/150‑1 (MQ=255)
                                                    cTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTcc                                                                                               <  1:2318921/150‑1 (MQ=255)
                                                    cTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTcc                                                                                               <  1:4897102/150‑1 (MQ=255)
                                                     tGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                                                                                              <  1:1748730/150‑1 (MQ=255)
                                                     tGTACGGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                                                                                              <  1:4287696/150‑1 (MQ=255)
                                                       aactcccATAACCACAATAATACTCTCGAAAAAAAATCGACTGCCAAATCTCCCACACACACCCGCCAGCTGAACCAGTCGATATTTTTTGAATATTATAAGAGCTATAAAAAACGGTGCGAACGCTGATTTCTTAAGCACATTTCCGca                                                                                            <  2:17928/144‑1 (MQ=255)
                                                        aCTGGCATCACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcac                                                                                           <  2:4030698/150‑1 (MQ=255)
                                                        aCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcac                                                                                           <  2:685237/150‑1 (MQ=255)
                                                          tGGAATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTTAACCTTTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACaa                                                                                         <  1:4874409/150‑1 (MQ=255)
                                                                aaCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATc                                                                                   <  2:4370185/150‑1 (MQ=255)
                                                                aaCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATc                                                                                   <  2:3734118/150‑1 (MQ=255)
                                                                     cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCAt                                                                              <  1:118520/150‑1 (MQ=255)
                                                                     aaTTTATACTGTCTTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCAt                                                                              <  1:2288115/149‑1 (MQ=255)
                                                                       tttATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTc                                                                            <  1:4172792/150‑1 (MQ=255)
                                                                              tGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgt                                                                     <  2:631373/150‑1 (MQ=255)
                                                                                   tataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACt                                                                <  1:2176042/150‑1 (MQ=255)
                                                                                   tataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACt                                                                <  1:5223445/150‑1 (MQ=255)
                                                                                    ataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTTAAATTTATAAAAGCTATAAAAAAACGTGCGAAAGCGGTTTTGTTAAGCAAAATTGAGCACAAAAAATATTAATTAGTTCTGTGGACTg                                                               >  2:5300489/1‑150 (MQ=255)
                                                                                    ataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTg                                                               <  1:1795985/150‑1 (MQ=255)
                                                                                    ataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTg                                                               <  2:314626/150‑1 (MQ=255)
                                                                                       aaaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAg                                                            <  1:3878899/150‑1 (MQ=255)
                                                                                        aaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAgg                                                           <  2:736675/150‑1 (MQ=255)
                                                                                          ttCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTTAATTTTATAATAGCTATAAAAAAAGGTGGGAAAGCGGTTTTCTTAATCACATTTTCGCACAACGTATATTCATTCGTTCTGTGGTCTGCAGGCt                                                         >  2:4498084/1‑150 (MQ=255)
                                                                                              aCTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTaa                                                     <  2:3163664/150‑1 (MQ=255)
                                                                                              aCTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTaa                                                     <  1:3131362/150‑1 (MQ=255)
                                                                                              aCTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTaa                                                     <  2:1906389/150‑1 (MQ=255)
                                                                                              aCTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTaa                                                     <  1:2984775/150‑1 (MQ=255)
                                                                                              aCTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTaa                                                     <  1:3456819/150‑1 (MQ=255)
                                                                                              aCTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTaa                                                     <  1:4702402/150‑1 (MQ=255)
                                                                                               cTGGCAATTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAAt                                                    <  1:3571203/150‑1 (MQ=255)
                                                                                               cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAAt                                                    <  2:964704/150‑1 (MQ=255)
                                                                                                tGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTGTTTTTTTTTATTTTTTAATAGCGATTAAAAAAGGTGCGAACGCGGTTTTTTTAATAAAATTTTCGCACAAAATATATTGATTAGTGGTGTGGGCGGCAGGCttttaat                                                   >  2:895453/1‑146 (MQ=255)
                                                                                                    aaagcgggCACTCTCTCCGGCCTGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATTAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATaa                                               <  1:1090507/144‑1 (MQ=255)
                                                                                                    aaaTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATaa                                               <  1:5460022/150‑1 (MQ=255)
                                                                                                       tCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGAt                                            <  1:5407122/150‑1 (MQ=255)
                                                                                                           gccctcgctCCGGCCAGGGGAACCCGTCGTTTTTTTTTGACTTTTATAAGCGCGATAAAAAACGGTGCGTACGCTGTTTTCTTCAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTtgt                                        <  2:3069785/147‑1 (MQ=255)
                                                                                                               tctctcCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTgcgc                                    <  1:5233743/150‑1 (MQ=255)
                                                                                                               tctctcCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTgcgc                                    <  2:1468376/150‑1 (MQ=255)
                                                                                                                      ggCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATAc                             <  2:1831632/150‑1 (MQ=255)
                                                                                                                        ccAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGt                           <  1:5806731/150‑1 (MQ=255)
                                                                                                                         cAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGtt                          <  2:1035674/150‑1 (MQ=255)
                                                                                                                          aGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGttt                         <  2:5037927/150‑1 (MQ=255)
                                                                                                                          aGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGttt                         <  2:5220227/150‑1 (MQ=255)
                                                                                                                             tGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGaa                      <  1:659591/150‑1 (MQ=255)
                                                                                                                             tGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGaa                      <  1:5786094/150‑1 (MQ=255)
                                                                                                                             tGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACATTTGaa                      <  2:4781372/150‑1 (MQ=255)
                                                                                                                                aCCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAatat                   <  1:1272757/150‑1 (MQ=255)
                                                                                                                                aCCAGTAGTTTTTTTTTGAATTTTATAGGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAatat                   <  2:3861847/150‑1 (MQ=255)
                                                                                                                                    ttCGTTTATTTTTGAATTTTATAAGAGCTATAAAATACGGTGCGAACGCTGTTTTCTTATGCACTTTTCCTCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATTTGATc               <  1:4457559/149‑1 (MQ=255)
                                                                                                                                    gTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATc               <  1:4624003/150‑1 (MQ=255)
                                                                                                                                      cGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCTCAACTTATCTTCTTTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCgg             <  2:231528/150‑1 (MQ=255)
                                                                                                                                      cGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCgg             <  1:146346/150‑1 (MQ=255)
                                                                                                                                      cGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCgg             <  2:764012/150‑1 (MQ=255)
                                                                                                                                      agTTTTTTTTTTAATTTTATAAGAGCTATAAAAAACGGTGCGAACGATGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCgg             <  2:5618188/149‑1 (MQ=255)
                                                                                                                                          tttttttGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATg         <  2:1705533/150‑1 (MQ=255)
                                                                                                                                          tttttttGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATg         >  1:3818947/1‑150 (MQ=255)
                                                                                                                                          tttttttGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATg         <  2:4726576/150‑1 (MQ=255)
                                                                                                                                           ttttttGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATgg        <  2:3639094/150‑1 (MQ=255)
                                                                                                                                            tttttGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGc       <  2:2510093/150‑1 (MQ=255)
                                                                                                                                                tGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATa   >  2:1560084/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GTGACCTTCCAGGACGACGCAGAACTGCCGTTGGCCTTCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACTGCAGGCTTTAATGATAAGATTTGTGCGCTAAATACGTTTGAATATGATCGGGATGGCAATAA  >  NC_000913/107409‑107703

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: