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)
           aGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAAttt                                                                                   <  2:847207/100‑1 (MQ=255)
             cTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTAGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAAtttt                                                                                  <  2:993495/99‑1 (MQ=255)
              tGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTGTCCGGCCGGGTGAACCAGTCGTTTTTTTTGAATTTtat                                                                                <  1:1196359/100‑1 (MQ=255)
               gTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTtata                                                                               >  2:8814/1‑100 (MQ=255)
                 aCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAag                                                                             <  1:1007456/100‑1 (MQ=255)
                 aCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATATGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAag                                                                             <  1:387535/100‑1 (MQ=255)
                  cTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAag                                                                             >  1:1862732/1‑100 (MQ=255)
                   tGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAaga                                                                            <  2:1862744/100‑1 (MQ=255)
                   tGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAagag                                                                           <  2:508447/100‑1 (MQ=255)
                    ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGc                                                                          <  1:1449585/100‑1 (MQ=255)
                        tAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCtata                                                                      >  2:968162/1‑100 (MQ=255)
                          aCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATaaa                                                                    >  2:706287/1‑100 (MQ=255)
                          aCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATaaa                                                                    >  2:166159/1‑100 (MQ=255)
                           cGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGTATTTTATAAGAGCTATaaaa                                                                   >  1:1538779/1‑100 (MQ=255)
                            gACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATaaaaa                                                                  <  1:1098729/100‑1 (MQ=255)
                              cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTAATTTTATAAGAGCTATaaaaa                                                                  >  1:1515015/1‑98 (MQ=255)
                              cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAAc                                                                <  2:1687755/100‑1 (MQ=255)
                              cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAAc                                                                <  2:944598/100‑1 (MQ=255)
                               aTTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACg                                                               <  2:150040/100‑1 (MQ=255)
                               aTTTATACTGTCGTATAAAAGTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATaaaaaagg                                                               >  1:1184642/1‑98 (MQ=255)
                                  tataCTGTCGTATAAAATTCGACTTGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGGCGGTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTg                                                            <  1:1783134/100‑1 (MQ=255)
                                   ataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAATACGGTGc                                                           <  2:1020652/100‑1 (MQ=255)
                                   ataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGc                                                           <  1:767588/100‑1 (MQ=255)
                                    taCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGTTGc                                                           >  1:161581/1‑99 (MQ=255)
                                    taCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCg                                                          <  2:284308/100‑1 (MQ=255)
                                      cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                        >  2:1019187/1‑100 (MQ=255)
                                       tGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAAc                                                       >  2:821351/1‑100 (MQ=255)
                                        gTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACg                                                      <  1:1677797/100‑1 (MQ=255)
                                        gTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACg                                                      <  1:879471/100‑1 (MQ=255)
                                         tCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGc                                                     <  1:568275/100‑1 (MQ=255)
                                          cGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTAATTTTATAAAAGCTATAAAAAACGGTGCGAACGc                                                     >  1:1571748/1‑99 (MQ=255)
                                          cGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCt                                                    >  1:230896/1‑100 (MQ=255)
                                           gTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTg                                                   >  2:1445686/1‑100 (MQ=255)
                                           gTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTg                                                   >  1:269783/1‑100 (MQ=255)
                                             ataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGtt                                                 >  1:1443495/1‑100 (MQ=255)
                                              taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATACGCGCTATAAAAAACGGTGCGAACGCTGttt                                                >  2:1135002/1‑100 (MQ=255)
                                                 aaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCt                                             >  2:605375/1‑100 (MQ=255)
                                                 aaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCt                                             >  2:607551/1‑100 (MQ=255)
                                                  aTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                            <  1:620667/100‑1 (MQ=255)
                                                   ttCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTa                                           <  2:406395/100‑1 (MQ=255)
                                                   ttCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTa                                           <  1:1235241/100‑1 (MQ=255)
                                                   ttCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTa                                           <  1:1002636/100‑1 (MQ=255)
                                                   ttCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTa                                           <  2:1381920/100‑1 (MQ=255)
                                                    tCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTaa                                          >  1:870818/1‑100 (MQ=255)
                                                     cGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAg                                         <  1:514979/100‑1 (MQ=255)
                                                      gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGCTTTCTTAAGc                                        >  2:647191/1‑100 (MQ=255)
                                                        cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                      >  2:1780724/1‑100 (MQ=255)
                                                        cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                      <  1:51651/100‑1 (MQ=255)
                                                        cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                      >  1:993377/1‑100 (MQ=255)
                                                         tGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACt                                     >  1:1775113/1‑100 (MQ=255)
                                                          ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                    >  2:1593885/1‑100 (MQ=255)
                                                           gCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACttt                                   >  1:299645/1‑100 (MQ=255)
                                                           gCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACttt                                   >  1:365329/1‑100 (MQ=255)
                                                             aaaTCTGGCACTCTCTCCGGCCAGGTGAACCTGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTc                                 <  2:230898/100‑1 (MQ=255)
                                                              actCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTcc                                <  1:1182192/98‑1 (MQ=255)
                                                                    gCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACaa                          >  2:976874/1‑100 (MQ=255)
                                                                     cACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAAc                         <  1:1899594/100‑1 (MQ=255)
                                                                        tctctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTa                      >  2:1458900/1‑100 (MQ=255)
                                                                         ctctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                     >  2:1480446/1‑100 (MQ=255)
                                                                         ctctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                     >  1:1104722/1‑100 (MQ=255)
                                                                            tcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCtt                  <  1:284044/100‑1 (MQ=255)
                                                                             ccGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTc                 >  1:1306012/1‑100 (MQ=255)
                                                                             ccGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTc                 >  2:1476274/1‑100 (MQ=255)
                                                                              cGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCa                >  2:1610440/1‑100 (MQ=255)
                                                                                gCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGATCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCAtt              <  2:875987/100‑1 (MQ=255)
                                                                                gCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTATTAAGCACTTTTCCGCACAACTTATCTTCAtt              >  1:282679/1‑100 (MQ=255)
                                                                                 ccAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTc             <  2:531605/100‑1 (MQ=255)
                                                                                 ccAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTc             >  2:630724/1‑100 (MQ=255)
                                                                                  cAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCg            <  2:1366886/100‑1 (MQ=255)
                                                                                   aGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGt           >  2:1495228/1‑100 (MQ=255)
                                                                                     gTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGc         <  2:1538787/100‑1 (MQ=255)
                                                                                          ccAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTgg    <  2:1467558/100‑1 (MQ=255)
                                                                                           cAGTCGTTTTTTTTGTTTTTTATAAAAGCTATAAAAATTGGTGCGAACGTTGTTTTCTTTAGTAATTTTACGCTCGACTTATCTTTATTTTTGCt        >  1:174545/1‑95 (MQ=255)
                                                                                            aGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAc  <  2:299647/100‑1 (MQ=255)
                                                                                            aGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAc  >  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: