Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 107,553 (T)9→8 100% intergenic (+79/‑152) 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,5450T.85.5% 260.3 / 48.5 83intergenic (+71/‑160)lpxC/secMUDP‑3‑O‑acyl N‑acetylglucosamine deacetylase/regulator of secA translation
Reads supporting (aligned to +/- strand):  ref base T (8/4);  new base . (34/37);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

TCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAC  >  NC_000913/107446‑107640
                                                                                                    |                                                                                               
tCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑t                                                                                                 >  2:1429707/1‑100 (MQ=255)
 cAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑tt                                                                                                <  1:1009465/100‑1 (MQ=255)
 cAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑tt                                                                                                >  2:1116847/1‑100 (MQ=255)
 cAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑tt                                                                                                >  1:1345562/1‑100 (MQ=255)
   aaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑tttt                                                                                              >  2:2171302/1‑100 (MQ=255)
    aGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑ttttt                                                                                             <  1:572634/100‑1 (MQ=255)
    aGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑ttttt                                                                                             >  1:113426/1‑100 (MQ=255)
     gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑tttttt                                                                                            <  1:1894346/100‑1 (MQ=255)
       gcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑tttttttt                                                                                          >  2:933592/1‑100 (MQ=255)
        ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGTACTCTCTCCGGCCAGGTGAACCAGTCG‑ttttttttg                                                                                         >  1:1629728/1‑99 (MQ=255)
         cTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑ttttttttga                                                                                        >  1:901426/1‑98 (MQ=255)
          ttCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑ttttttttgaa                                                                                       <  2:929486/100‑4 (MQ=255)
          ttCAGCTGTACTGGCATAACGACATTTATACTGTCGCATAAAATTCTACTGTCAAATCTTGCACTCTCACCGGCCAGGTGAACCAGTCG‑ttttttttgaa                                                                                       <  2:1310651/100‑4 (MQ=255)
           tCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTTGaa                                                                                      <  1:251577/100‑1 (MQ=255)
              gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTTGAAttt                                                                                   >  1:2175239/1‑100 (MQ=255)
              gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAAtttt                                                                                  >  1:1267086/1‑100 (MQ=255)
              gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAAtttt                                                                                  <  2:928171/100‑1 (MQ=255)
               cTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTta                                                                                 <  2:1555690/100‑1 (MQ=255)
               cTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTta                                                                                 >  1:766359/1‑100 (MQ=255)
                tGTACTGGCATAACGACATTTATACTTTCGTATAAAATTCGACTTGCATTTCTTTCTTTCTCTCTTGTCAGTTGATCCAGTCG‑‑TTTTTTTTGAATTTtat                                                                                <  2:1377467/100‑1 (MQ=255)
                tGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTtat                                                                                >  1:264052/1‑100 (MQ=255)
                   ttccgatctaCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAag                                                                             <  2:1655688/91‑1 (MQ=255)
                    cTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAaga                                                                            >  2:1628002/1‑100 (MQ=255)
                     tGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAagag                                                                           <  1:1282788/100‑1 (MQ=255)
                     tGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAagag                                                                           <  1:50181/100‑1 (MQ=255)
                      ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGc                                                                          >  2:1662277/1‑100 (MQ=255)
                       gCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTTGAATTTTATAagag                                                                           <  2:1898866/100‑1 (MQ=255)
                         aTAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTT‑TTTTTTATAACAGCtat                                                                       >  2:2032486/1‑100 (MQ=255)
                           aaCGACATTTTTCCGGTCGTATAAAATTCGACTGGCAAATTTGGCATTCTCTCCGGCGAGGTGAACCGGTCG‑‑TTTTTTTTGAATTTTATAAGAGCta                                                                        <  2:615617/97‑1 (MQ=255)
                           aaCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCtata                                                                      <  1:2084376/99‑1 (MQ=255)
                            aCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAagagatcggaa                                                                    >  1:1655699/1‑93 (MQ=255)
                            aCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATaaa                                                                    >  1:359420/1‑100 (MQ=255)
                             cGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAATATCTATaaaa                                                                   >  1:1993632/1‑100 (MQ=255)
                                cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAAc                                                                <  1:218331/100‑1 (MQ=255)
                                cATTTATACTGTAGTATAAAATTCGACTGGCAAATCTTGCACTCTCTCCGGCAAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAAc                                                                <  2:2021100/100‑1 (MQ=255)
                                 aTTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTT‑TTTTTTTTAAAAGCTaaaaaa                                                                   >  2:558940/1‑96 (MQ=255)
                                 aTTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAAc                                                                <  1:507951/99‑1 (MQ=255)
                                 aTTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACg                                                               >  2:1507789/1‑100 (MQ=255)
                                  tttATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACgg                                                              >  1:278246/1‑100 (MQ=255)
                                  tttATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTAAATTTTATAAGAGCTATAAAAAACgg                                                              >  1:1832388/1‑100 (MQ=255)
                                   ttatacgacagTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGt                                                             >  2:733529/11‑100 (MQ=255)
                                   ttATACTGTCGGATAAAATTCGATTGGTAATTCTAGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGt                                                             <  2:163453/100‑1 (MQ=255)
                                    tataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTg                                                            <  1:682013/100‑1 (MQ=255)
                                     ataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGc                                                           <  1:829747/100‑1 (MQ=255)
                                     ataCTGTCATATAAAATTCGACTGGCAAATCTGGAACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGc                                                           <  1:829953/100‑1 (MQ=255)
                                        cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                        >  1:320608/1‑100 (MQ=255)
                                          gTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACg                                                      >  1:72632/1‑100 (MQ=255)
                                             gTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTg                                                   >  1:425189/1‑100 (MQ=255)
                                              tataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGt                                                  >  1:74382/1‑100 (MQ=255)
                                               ataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGtt                                                 >  1:1391814/1‑100 (MQ=255)
                                                taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                                <  2:467358/100‑1 (MQ=255)
                                                taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                                >  1:1113516/1‑100 (MQ=255)
                                                   aaTTCGACTGGCAAATCTGGTACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCt                                             <  2:1582367/100‑1 (MQ=255)
                                                   aaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTT‑TATTTTATAAGAGATATAAAAAACGGTGCGAACGCTGTTTTCt                                             >  2:1948244/1‑100 (MQ=255)
                                                     ttCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTa                                           <  2:1267080/100‑1 (MQ=255)
                                                      tCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTaa                                          <  1:733529/100‑1 (MQ=255)
                                                       cGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTT‑AATTTTATAATTTCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAg                                         >  1:1269852/1‑100 (MQ=255)
                                                       cGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTaa                                          >  1:585250/1‑99 (MQ=255)
                                                        gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGc                                        >  1:17563/1‑100 (MQ=255)
                                                         aCTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCa                                       <  2:864786/100‑1 (MQ=255)
                                                          cTGTCAAATCTTGCACTCTCTCCTGCCAGTTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAATAGCTATAAAAAACGGTTCGAACGCTGTTTTCTTAAGCAc                                      <  2:1129945/100‑1 (MQ=255)
                                                          cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGc                                        <  1:975913/99‑1 (MQ=255)
                                                          cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                      <  2:602011/100‑1 (MQ=255)
                                                            ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                    >  2:1175759/1‑100 (MQ=255)
                                                            ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                    >  1:300825/1‑100 (MQ=255)
                                                            ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                    >  2:2167781/1‑100 (MQ=255)
                                                              cAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtttt                                  >  1:1225762/1‑100 (MQ=255)
                                                               aaaTCTGGAACTCTCTCAGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTc                                 <  2:1013544/100‑1 (MQ=255)
                                                                  tCTGGCTCTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                              <  1:1793471/100‑1 (MQ=255)
                                                                  tCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                              <  2:711410/100‑1 (MQ=255)
                                                                  tCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                              <  1:556794/100‑1 (MQ=255)
                                                                     ggCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcaca                           >  2:1396088/1‑100 (MQ=255)
                                                                      gCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcac                            <  1:226189/99‑1 (MQ=255)
                                                                      gCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCAGCACaa                          >  1:79086/1‑100 (MQ=255)
                                                                      gCACTCTCTCCGGCCAAGTGAAACAATCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGGGAGAACGCATATTTCTTAAGCACTTTTCCACACaa                          >  2:980497/1‑100 (MQ=255)
                                                                       cACTCTCTTCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAAc                         <  1:1219042/100‑1 (MQ=255)
                                                                       cACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAAc                         >  1:915161/1‑100 (MQ=255)
                                                                         ctctctCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACtt                       >  1:611392/1‑100 (MQ=255)
                                                                          tctctcCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTa                      <  2:1840547/100‑1 (MQ=255)
                                                                           ctctcCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                     >  2:656369/1‑100 (MQ=255)
                                                                           ctctcCGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                     >  1:236353/1‑100 (MQ=255)
                                                                              tcCGGCCAGGTGAACCAGTCG‑TTTTTTTTT‑TTTTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCtt                  >  2:151348/1‑100 (MQ=255)
                                                                               ccGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTc                 <  1:1436722/100‑1 (MQ=255)
                                                                                cGGCCAGGTGAACCAGTCG‑TTTTTTTTT‑AATTTTATAAGAGCTATAAAAAACGGTACTAACGCTGTTTTCTTAAGCAATTTTCCGCACAACTTATCTTc                 >  2:1578417/1‑99 (MQ=255)
                                                                                cGGCCAGGTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCa                <  1:626678/100‑1 (MQ=255)
                                                                                      ggTGAACCAGTCG‑‑TTTTTTTTGTATTTTTTAAAAGCTATAAAAAACCGTGCGAACGATGTTTTATTAAGCACTTTTCAGCACAACTTCTATTCATTCGta          >  2:1776230/1‑99 (MQ=255)
                                                                                      ggTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTg          >  1:2226247/1‑100 (MQ=255)
                                                                                      ggTGAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTg          <  2:1754721/100‑1 (MQ=255)
                                                                                         gAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtg       <  1:301337/100‑1 (MQ=255)
                                                                                         gAACCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtg       >  1:709911/1‑100 (MQ=255)
                                                                                         gAACCAGTCG‑‑TTAGTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGc         <  1:1388426/98‑1 (MQ=255)
                                                                                           aCCAGTCG‑TTTTTTTTT‑TTTTTTATAAAAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgt      >  2:1203865/1‑99 (MQ=255)
                                                                                           aCCAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgtg     <  2:1345556/100‑1 (MQ=255)
                                                                                             cAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgtg     <  2:1957561/98‑1 (MQ=255)
                                                                                             cAGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGa   <  2:1052776/100‑1 (MQ=255)
                                                                                              aGTCG‑‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAc  <  1:1628013/100‑1 (MQ=255)
                                                                                                    |                                                                                               
TCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAC  >  NC_000913/107446‑107640

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