Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 107,553 (T)9→8 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.93.8% 454.6 / 25.0 112intergenic (+71/‑160)lpxC/secMUDP‑3‑O‑acyl N‑acetylglucosamine deacetylase/regulator of secA translation
Reads supporting (aligned to +/- strand):  ref base T (7/0);  new base . (48/57);  total (55/57)
Fisher's exact test for biased strand distribution p-value = 5.60e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

AAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAC  >  NC_000913/107448‑107640
                                                                                                 |                                                                                               
aaaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttt                                                                                               >  2:1852913/1‑100 (MQ=255)
 aaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttt                                                                                              >  2:991314/1‑100 (MQ=255)
 aaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttt                                                                                              >  1:348660/1‑100 (MQ=255)
 aaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttt                                                                                              >  1:364322/1‑100 (MQ=255)
 aaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtt                                                                                                >  2:2480888/1‑98 (MQ=255)
  aGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttt                                                                                             >  2:1189719/1‑100 (MQ=255)
   gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttttt                                                                                            >  1:1989500/1‑100 (MQ=255)
   gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttttt                                                                                            >  1:485552/1‑100 (MQ=255)
    cgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttt                                                                                           <  2:1367725/100‑1 (MQ=255)
    cgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttt                                                                                           <  1:2441292/100‑1 (MQ=255)
    cgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttn                                                                                           <  2:2379960/100‑2 (MQ=255)
      ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGTCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttg                                                                                         >  1:781933/1‑99 (MQ=255)
      ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttt                                                                                         <  2:2542525/100‑1 (MQ=255)
      ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttg                                                                                         >  2:37206/1‑99 (MQ=255)
      ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttttttt                                                                                          >  2:1618599/1‑99 (MQ=255)
           aGCTGTACTGGCATATCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAAttt                                                                                   <  1:37207/100‑1 (MQ=255)
           aGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAAttt                                                                                   <  1:1489961/100‑1 (MQ=255)
            gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTT‑AAtttt                                                                                  >  2:1380173/1‑100 (MQ=255)
             cTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTta                                                                                 <  2:1033754/100‑1 (MQ=255)
              tGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTtat                                                                                <  1:239515/100‑1 (MQ=255)
              tGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTtat                                                                                >  2:1351938/1‑100 (MQ=255)
              tGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTtat                                                                                <  2:527582/100‑1 (MQ=255)
              ngTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTT‑AATTTtat                                                                                >  2:2379910/2‑100 (MQ=255)
               gTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTtata                                                                               <  2:1552747/100‑1 (MQ=255)
                 aCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAag                                                                             >  2:351255/1‑100 (MQ=255)
                   tGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAagag                                                                           <  2:2540869/100‑1 (MQ=255)
                   tGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGCGAACCAGTCG‑TTTTTTTTGAATTTTATAagag                                                                           <  1:1364986/100‑1 (MQ=255)
                    ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGc                                                                          >  1:917962/1‑100 (MQ=255)
                    ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGc                                                                          >  2:2255013/1‑100 (MQ=255)
                      cATAACGACATTTATACTGTCGAATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCta                                                                        <  1:1368640/100‑1 (MQ=255)
                      cATAACGACATTTATACCGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCta                                                                        <  1:829411/100‑1 (MQ=255)
                       aTAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTTAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCtat                                                                       <  2:826449/100‑1 (MQ=255)
                       aTAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCtat                                                                       >  1:934308/1‑100 (MQ=255)
                       aTAACGACATTTATACTGTCGTATAAAATTCGACTGGAAAATCTTGCATCATCTTCGGGCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCta                                                                        <  2:760305/99‑1 (MQ=255)
                          aCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATaaa                                                                    >  2:2424155/1‑100 (MQ=255)
                          aCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATaaa                                                                    >  1:1268033/1‑100 (MQ=255)
                            gACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCAAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATaaaaa                                                                  <  2:1595973/100‑1 (MQ=255)
                             aCATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATaaaaa                                                                  <  2:1989483/99‑1 (MQ=255)
                                tttATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACgg                                                              >  2:1470684/1‑100 (MQ=255)
                                tttATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACgg                                                              >  1:2279550/1‑100 (MQ=255)
                                 ttATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGt                                                             <  2:348656/100‑1 (MQ=255)
                                  tataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGACTTTTATAAGAGATATAAAAAACGGTg                                                            >  2:1273580/1‑100 (MQ=255)
                                  tataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTg                                                            >  1:2286330/1‑100 (MQ=255)
                                  tataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTg                                                            <  1:239948/100‑1 (MQ=255)
                                    taCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCg                                                          >  1:2543314/1‑100 (MQ=255)
                                      cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGTATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                        >  2:2703490/1‑100 (MQ=255)
                                        gTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACg                                                      >  2:427163/1‑100 (MQ=255)
                                        gTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACg                                                      >  1:1590784/1‑100 (MQ=255)
                                          tgTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCt                                                    <  1:2669630/99‑1 (MQ=255)
                                           gTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCta                                                   >  2:2669610/1‑99 (MQ=255)
                                            tataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCt                                                    <  2:1283399/98‑1 (MQ=255)
                                            tataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTg                                                   <  1:2463416/99‑1 (MQ=255)
                                            tataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGt                                                  <  2:1115127/100‑1 (MQ=255)
                                             ataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGtt                                                 >  2:831555/1‑100 (MQ=255)
                                              taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                                >  1:2590098/1‑100 (MQ=255)
                                              taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                                >  2:701531/1‑100 (MQ=255)
                                                 aaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCt                                             <  1:2301314/100‑1 (MQ=255)
                                                 aaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCt                                             <  1:1234786/100‑1 (MQ=255)
                                                  aTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                            <  2:2124621/100‑1 (MQ=255)
                                                   ttCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTa                                           <  2:75014/100‑1 (MQ=255)
                                                    tCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTT‑TTTTTTTTAAGAGATATAAAAAACGGTGCGAACGCTGTTTTCtt                                            >  2:2725746/1‑98 (MQ=255)
                                                      gacttactaaTCTCGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGc                                        <  2:118839/92‑1 (MQ=255)
                                                      gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGc                                        <  1:991321/100‑1 (MQ=255)
                                                      gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGc                                        <  2:1435113/100‑1 (MQ=255)
                                                      gACTGGCAAATCTCGCACTCTCACCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGc                                        <  1:2480907/100‑1 (MQ=255)
                                                         tGGCATATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACt                                     <  2:980710/100‑1 (MQ=255)
                                                         tGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACt                                     >  1:2538603/1‑100 (MQ=255)
                                                         tGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACt                                     >  1:2403483/1‑100 (MQ=255)
                                                          ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                    >  2:2495323/1‑100 (MQ=255)
                                                          ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                    >  2:2649770/1‑100 (MQ=255)
                                                          ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                    >  2:2134748/1‑100 (MQ=255)
                                                             ataTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAAAAAAAAAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTc                                 <  2:1672304/98‑1 (MQ=255)
                                                               aTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGATTTTTATAAGAGCTATACAAAACGGTTCGAACGCTGTTTTCTTAAGCACTTTTCCg                               >  2:2223276/1‑100 (MQ=255)
                                                               aTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                               >  1:214044/1‑100 (MQ=255)
                                                               aTCTGGCACTCTCTCCGGCAAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                               <  1:855750/100‑1 (MQ=255)
                                                                ttccgatctctctCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                              <  1:1924977/93‑1 (MQ=255)
                                                                ttccgatctctctCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                              <  1:1924864/93‑1 (MQ=255)
                                                                tCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                              <  1:2116646/100‑1 (MQ=255)
                                                                    gCACTCTCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAAc                         >  1:939248/1‑100 (MQ=255)
                                                                     cACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAAc                         <  1:1591771/100‑1 (MQ=255)
                                                                      aCTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAGCt                        >  2:2061475/1‑100 (MQ=255)
                                                                       ctctctCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACtt                       <  2:934302/100‑1 (MQ=255)
                                                                         ctctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcagatcggaa                     >  2:1924961/1‑92 (MQ=255)
                                                                         ctctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcagatcggaa                     >  2:1924848/1‑92 (MQ=255)
                                                                         ctctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                     >  2:1038711/1‑100 (MQ=255)
                                                                           ctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCt                   <  1:709910/100‑1 (MQ=255)
                                                                                  cAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCg            >  1:381576/1‑100 (MQ=255)
                                                                                  cAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCACACAACTTATCTTCATTCg            <  1:1528493/100‑1 (MQ=255)
                                                                                    ggTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTg          <  1:105261/100‑1 (MQ=255)
                                                                                    ggTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTg          <  2:2590079/100‑1 (MQ=255)
                                                                                     gTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGc         <  2:478500/100‑1 (MQ=255)
                                                                                     gTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGc         <  1:260124/100‑1 (MQ=255)
                                                                                     gTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGc         >  1:422924/1‑100 (MQ=255)
                                                                                      tGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCt        >  1:2076253/1‑100 (MQ=255)
                                                                                       gAACCAGTCGTTTTTTTTT‑TTTTTTTTAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCt        >  2:289010/1‑99 (MQ=255)
                                                                                         aCCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgtg     >  2:2291891/1‑100 (MQ=255)
                                                                                          ccAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAATGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTgg    >  1:276863/1‑100 (MQ=255)
                                                                                          ccAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTgg    <  1:1038718/100‑1 (MQ=255)
                                                                                           cAGTCGTTTTTTTTT‑TTTTTTTTAAGAGCTATAAAAACCGGTGCGAACTCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgt      >  2:190977/1‑97 (MQ=255)
                                                                                           cAGTCGTTTTTTTTT‑TTTTTTTTAAGAGCTATAAAAAACGGTGCGAACTCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCttt      >  2:890585/1‑95 (MQ=255)
                                                                                                 |                                                                                               
AAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAC  >  NC_000913/107448‑107640

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