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.94.7% 311.5 / 10.7 75intergenic (+71/‑160)lpxC/secMUDP‑3‑O‑acyl N‑acetylglucosamine deacetylase/regulator of secA translation
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base . (36/35);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

AAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAC  >  NC_000913/107448‑107640
                                                                                                 |                                                                                               
aaaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttt                                                                                               >  1:1458451/1‑100 (MQ=255)
aaaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttt                                                                                               >  1:1458420/1‑100 (MQ=255)
aaaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttt                                                                                               >  1:944597/1‑100 (MQ=255)
 aaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCATATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttt                                                                                              >  1:1059849/1‑100 (MQ=255)
 aaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttt                                                                                              <  1:825455/100‑1 (MQ=255)
    cgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttt                                                                                           <  1:1550069/100‑1 (MQ=255)
    cgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttt                                                                                           <  2:1444870/100‑1 (MQ=255)
     gcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttttttt                                                                                          <  1:379830/100‑1 (MQ=255)
     gcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttttttt                                                                                          <  1:415819/100‑1 (MQ=255)
     gcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttttttt                                                                                          >  2:577110/1‑100 (MQ=255)
      ccTTGAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttt                                                                                           <  1:135022/98‑1 (MQ=255)
       cTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttttttt                                                                                          >  2:224435/1‑98 (MQ=255)
        ttCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttgaa                                                                                       <  2:1448590/100‑4 (MQ=255)
         tCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttgaat                                                                                      >  2:514981/1‑96 (MQ=255)
           aGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAAttt                                                                                   <  2:847207/100‑1 (MQ=255)
             cTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTAGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAAtttt                                                                                  <  2:993495/99‑1 (MQ=255)
              tGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTGTCCGGCCGGGTGAACCAGTCG‑TTTTTTTTGAATTTtat                                                                                <  1:1196359/100‑1 (MQ=255)
               gTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTtata                                                                               >  2:8814/1‑100 (MQ=255)
                 aCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAag                                                                             <  1:1007456/100‑1 (MQ=255)
                 aCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATATGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAag                                                                             <  1:387535/100‑1 (MQ=255)
                  cTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAag                                                                             >  1:1862732/1‑100 (MQ=255)
                   tGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAaga                                                                            <  2:1862744/100‑1 (MQ=255)
                   tGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAagag                                                                           <  2:508447/100‑1 (MQ=255)
                    ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGc                                                                          <  1:1449585/100‑1 (MQ=255)
                        tAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCtata                                                                      >  2:968162/1‑100 (MQ=255)
                          aCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATaaa                                                                    >  2:706287/1‑100 (MQ=255)
                          aCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATaaa                                                                    >  2:166159/1‑100 (MQ=255)
                           cGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGTATTTTATAAGAGCTATaaaa                                                                   >  1:1538779/1‑100 (MQ=255)
                            gACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATaaaaa                                                                  <  1:1098729/100‑1 (MQ=255)
                              cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTT‑AATTTTATAAGAGCTATaaaaa                                                                  >  1:1515015/1‑98 (MQ=255)
                              cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAAc                                                                <  2:1687755/100‑1 (MQ=255)
                              cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAAc                                                                <  2:944598/100‑1 (MQ=255)
                               aTTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACg                                                               <  2:150040/100‑1 (MQ=255)
                               aTTTATACTGTCGTATAAAAGTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATaaaaaagg                                                               >  1:1184642/1‑98 (MQ=255)
                                  tataCTGTCGTATAAAATTCGACTTGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGGCG‑GTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTg                                                            <  1:1783134/100‑1 (MQ=255)
                                   ataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAATACGGTGc                                                           <  2:1020652/100‑1 (MQ=255)
                                   ataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGc                                                           <  1:767588/100‑1 (MQ=255)
                                    taCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGTTGc                                                           >  1:161581/1‑99 (MQ=255)
                                    taCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCg                                                          <  2:284308/100‑1 (MQ=255)
                                      cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                        >  2:1019187/1‑100 (MQ=255)
                                       tGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAAc                                                       >  2:821351/1‑100 (MQ=255)
                                        gTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACg                                                      <  1:1677797/100‑1 (MQ=255)
                                        gTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACg                                                      <  1:879471/100‑1 (MQ=255)
                                         tCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGc                                                     <  1:568275/100‑1 (MQ=255)
                                          cGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTT‑AATTTTATAAAAGCTATAAAAAACGGTGCGAACGc                                                     >  1:1571748/1‑99 (MQ=255)
                                          cGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCt                                                    >  1:230896/1‑100 (MQ=255)
                                           gTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTg                                                   >  2:1445686/1‑100 (MQ=255)
                                           gTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTg                                                   >  1:269783/1‑100 (MQ=255)
                                             ataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGtt                                                 >  1:1443495/1‑100 (MQ=255)
                                              taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATACGCGCTATAAAAAACGGTGCGAACGCTGttt                                                >  2:1135002/1‑100 (MQ=255)
                                                 aaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCt                                             >  2:605375/1‑100 (MQ=255)
                                                 aaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCt                                             >  2:607551/1‑100 (MQ=255)
                                                  aTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                            <  1:620667/100‑1 (MQ=255)
                                                   ttCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTa                                           <  2:406395/100‑1 (MQ=255)
                                                   ttCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTa                                           <  1:1235241/100‑1 (MQ=255)
                                                   ttCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTa                                           <  1:1002636/100‑1 (MQ=255)
                                                   ttCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTa                                           <  2:1381920/100‑1 (MQ=255)
                                                    tCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTaa                                          >  1:870818/1‑100 (MQ=255)
                                                     cGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAg                                         <  1:514979/100‑1 (MQ=255)
                                                      gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGCTTTCTTAAGc                                        >  2:647191/1‑100 (MQ=255)
                                                        cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                      >  2:1780724/1‑100 (MQ=255)
                                                        cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                      <  1:51651/100‑1 (MQ=255)
                                                        cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                      >  1:993377/1‑100 (MQ=255)
                                                         tGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACt                                     >  1:1775113/1‑100 (MQ=255)
                                                          ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                    >  2:1593885/1‑100 (MQ=255)
                                                           gCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACttt                                   >  1:299645/1‑100 (MQ=255)
                                                           gCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACttt                                   >  1:365329/1‑100 (MQ=255)
                                                             aaaTCTGGCACTCTCTCCGGCCAGGTGAACCTGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTc                                 <  2:230898/100‑1 (MQ=255)
                                                              actCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTcc                                <  1:1182192/98‑1 (MQ=255)
                                                                    gCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACaa                          >  2:976874/1‑100 (MQ=255)
                                                                     cACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAAc                         <  1:1899594/100‑1 (MQ=255)
                                                                        tctctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTa                      >  2:1458900/1‑100 (MQ=255)
                                                                         ctctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                     >  2:1480446/1‑100 (MQ=255)
                                                                         ctctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                     >  1:1104722/1‑100 (MQ=255)
                                                                            tcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCtt                  <  1:284044/100‑1 (MQ=255)
                                                                             ccGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTc                 >  1:1306012/1‑100 (MQ=255)
                                                                             ccGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTc                 >  2:1476274/1‑100 (MQ=255)
                                                                              cGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCa                >  2:1610440/1‑100 (MQ=255)
                                                                                gCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGATCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCAtt              <  2:875987/100‑1 (MQ=255)
                                                                                gCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTATTAAGCACTTTTCCGCACAACTTATCTTCAtt              >  1:282679/1‑100 (MQ=255)
                                                                                 ccAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTc             <  2:531605/100‑1 (MQ=255)
                                                                                 ccAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTc             >  2:630724/1‑100 (MQ=255)
                                                                                  cAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCg            <  2:1366886/100‑1 (MQ=255)
                                                                                   aGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGt           >  2:1495228/1‑100 (MQ=255)
                                                                                     gTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGc         <  2:1538787/100‑1 (MQ=255)
                                                                                          ccAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTgg    <  2:1467558/100‑1 (MQ=255)
                                                                                           cAGTCG‑TTTTTTTTGTTTTTTATAAAAGCTATAAAAATTGGTGCGAACGTTGTTTTCTTTAGTAATTTTACGCTCGACTTATCTTTATTTTTGCt        >  1:174545/1‑95 (MQ=255)
                                                                                            aGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAc  <  2:299647/100‑1 (MQ=255)
                                                                                            aGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAc  >  2:1196013/1‑100 (MQ=255)
                                                                                                 |                                                                                               
AAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAC  >  NC_000913/107448‑107640

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