Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,725,949 G→A 5.1% M21I (ATG→ATA wecH → O‑acetyltransferase for enterobacterial common antigen (ECA)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,725,9490GA 5.1% 258.5 / 2.2 79M21I (ATG→ATAwecHO‑acetyltransferase for enterobacterial common antigen (ECA)
Reads supporting (aligned to +/- strand):  ref base G (34/40);  new base A (2/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

TAGTGATAAAACTCAGATTATTCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGTTCTGAATTCTGCCTCTCGTGTCAGCGTGCCGCTATTTTTCATGATTTCCGGCTATCTCTTTTTTGGCGAACG  >  NC_000913/3725804‑3726092
                                                                                                                                                 |                                                                                                                                               
tAGTGATAAAACTCAGATTATTCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTgatgat                                                                                                                                               >  3:240144/1‑148 (MQ=255)
  gTGATAAAACTCAGATTATTCCAGTGCGCCCACCAAATGATTTATCGGTCCGACTTTTTTTAACCAGATACAGAGCCAAAGTTGCAGCCCATAATTTACTGGATTGATAACCTGCGAGGGATAGAGTGTTTAATGGTTGTTATGAtt                                                                                                                                              <  6:297256/147‑1 (MQ=255)
  gTGATAAAACTCAGATTATTCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGAtt                                                                                                                                              <  2:370783/147‑1 (MQ=255)
  gTGATAAAACTCAGATTATTCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGAtt                                                                                                                                              <  4:192655/147‑1 (MQ=255)
   tGATAAAACTCAGATTATTCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTca                                                                                                                                            >  6:267818/1‑148 (MQ=255)
    gATAAAACTCAGATTATTCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTca                                                                                                                                            >  1:249160/1‑147 (MQ=255)
    gATAAAACTCAGATTATTCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTca                                                                                                                                            >  5:310636/1‑147 (MQ=255)
     aTAAAACTCAGATTATTCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCCAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTcaca                                                                                                                                          <  2:307610/148‑1 (MQ=255)
     aTAAAACTCAGATTATTCCAGTGCGCCAACCAAATGATTTATCGGGCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTcaca                                                                                                                                          <  7:129603/148‑1 (MQ=255)
      tAAAACTCAGATTATTCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCacac                                                                                                                                         <  5:267818/148‑1 (MQ=255)
      tAAAACTCAGATTATTCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCacac                                                                                                                                         >  7:187720/1‑148 (MQ=255)
        aaaCTCAGATTATTCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCAcacc                                                                                                                                        <  4:240144/147‑1 (MQ=255)
        aaaCTCAGATTATTCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCAcacc                                                                                                                                        <  4:120952/147‑1 (MQ=255)
        aaaCTCAGATTATTCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCAcacc                                                                                                                                        <  2:249160/147‑1 (MQ=255)
             cAGATTATTCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTAcc                                                                                                                                  >  1:358242/1‑148 (MQ=255)
              aGATTATTCCAGTGCTCCAACATAATGTTTTATCGTTAATACTTTTTTTAACCAGATACAGAGACACAGATGCAGAACAAAATTTACTGGATTGATAACCTGCGAGTGATAGCGTGTTTAATGGTGGTGATGATTaaca                                                                                                                                          <  4:92658/139‑1 (MQ=255)
              aGATTATTCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACAGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTAcc                                                                                                                                  <  8:12897/147‑1 (MQ=255)
              aGATTATTCCAGTGCGCCAAACAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCt                                                                                                                                 <  6:56214/148‑1 (MQ=255)
               gATTATTCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGACACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCAcaccagatcgt                                                                                                                                 >  3:92658/1‑141 (MQ=255)
                cttattCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGAAGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGGGTTTAATGGTGGTGATGATTCACACCACTACCTgg                                                                                                                               <  1:253460/147‑1 (MQ=255)
                attattCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCAcacc                                                                                                                                        <  4:116460/139‑1 (MQ=255)
                attattCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCAcacc                                                                                                                                        >  3:116460/1‑139 (MQ=255)
                      ccAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAACCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTAtgt                                                                                                                          <  3:114632/147‑1 (MQ=255)
                      ccAGTGCGCCAACCAAATGATTTATCCGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTAATGGATTGTTAACCTGCGAGGGATAGCCTGTTTAAAGGAGCTGATCATACACAACACTACCTGGTAtgt                                                                                                                          >  5:29382/1‑147 (MQ=255)
                       cAGTGCTCCAACCAAATTATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGACCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATTTGa                                                                                                                        <  6:64798/148‑1 (MQ=255)
                         gTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGAAACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGAcc                                                                                                                      <  2:378448/148‑1 (MQ=255)
                              ccaaccaaATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGAGACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGc                                                                                                                 <  1:46560/148‑1 (MQ=255)
                                   caaATGATTTATCTGTCAGACTTTTTTTAAAAAGATACAGAGACACAGATGCAGCCCAAAATTTACAGGATTGATAACCTTAGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACGTTGTATGTGAACAATGCTCa                                                                                                              >  2:331856/1‑146 (MQ=255)
                                   caaATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCa                                                                                                              >  5:327861/1‑146 (MQ=255)
                                   caaATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCa                                                                                                              >  8:2797/1‑146 (MQ=255)
                                     aaTGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTcaca                                                                                                                                          >  1:221451/1‑116 (MQ=255)
                                     aaTGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTcaca                                                                                                                                          <  2:221451/116‑1 (MQ=255)
                                      aTGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTcaca                                                                                                                                          <  6:167276/115‑1 (MQ=255)
                                      aTGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTcaca                                                                                                                                          >  5:167276/1‑115 (MQ=255)
                                      aTGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACAACCTGGTATGTGCCCAATGCTCATAgtg                                                                                                         >  1:128514/1‑148 (MQ=255)
                                        gATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCAt                                                                                                             <  3:370958/142‑1 (MQ=255)
                                        gATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCAt                                                                                                             >  4:370958/1‑142 (MQ=255)
                                          tttATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAg                                                                                                     >  6:144702/1‑148 (MQ=255)
                                          tttATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAg                                                                                                     >  2:389663/1‑148 (MQ=255)
                                            tATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGc                                                                                                    <  8:187720/147‑1 (MQ=255)
                                              tCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGccc                                                                                                  <  7:176074/147‑1 (MQ=255)
                                                  tCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCAcacc                                                                                                                                        >  1:162032/1‑105 (MQ=255)
                                                  tCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCAcacc                                                                                                                                        <  2:162032/105‑1 (MQ=255)
                                                  tCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACAAATGCTCATAGTGTTAGCCCCGTca                                                                                             >  4:46764/1‑148 (MQ=255)
                                                   cAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTa                                                                                                                                    >  1:322825/1‑108 (MQ=255)
                                                   cAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTa                                                                                                                                    <  2:322825/108‑1 (MQ=255)
                                                      aCTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATg                                                                                         <  5:277077/148‑1 (MQ=255)
                                                      aCTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATg                                                                                         >  2:13872/1‑148 (MQ=255)
                                                         ttttttAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTTCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGa                                                                                      >  5:31433/1‑148 (MQ=255)
                                                         ttttttAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTTCGAGGGATAGAGTGTTTAATGTTGGTGATGATTAACACCACTACCTGGTATGTTACCAATGCTCATAGTGTTATCCCCTTCACATGGGa                                                                                      >  7:110876/1‑148 (MQ=255)
                                                         ttttttAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGCTGATTCACACCACTACCTGGTATGTGACCAATGCTCCTAGTTTTAGCCCCGTCACAGGGGa                                                                                      >  2:72652/1‑148 (MQ=255)
                                                         ttttttAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGa                                                                                      >  2:98462/1‑148 (MQ=255)
                                                                 acAGATACAGCGACACAGATGCAGCCCAAAATTTACTTGATTTATTACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCa                                                                              <  2:128514/147‑1 (MQ=255)
                                                                    gATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTca                                                                                             <  7:413911/130‑1 (MQ=255)
                                                                    gATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTca                                                                                             >  8:413911/1‑130 (MQ=255)
                                                                     aTACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGt                                                                          <  6:310636/148‑1 (MQ=255)
                                                                     aTACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATAATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGt                                                                          >  3:218439/1‑148 (MQ=255)
                                                                     aTACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATAATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGt                                                                          >  3:213346/1‑148 (MQ=255)
                                                                        cAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGTTc                                                                        <  2:358242/147‑1 (MQ=255)
                                                                         agagACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGTTCTg                                                                      <  7:2797/148‑1 (MQ=255)
                                                                          gagaCACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGTTCTGa                                                                     <  1:389663/148‑1 (MQ=255)
                                                                            gaCACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGTTCTGAAt                                                                   <  7:254572/148‑1 (MQ=255)
                                                                                caGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAg                                                                                                           >  7:261503/1‑104 (MQ=255)
                                                                                caGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAg                                                                                                           <  8:261503/104‑1 (MQ=255)
                                                                                caGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGTTCTGAATTCTg                                                               <  6:62462/148‑1 (MQ=255)
                                                                                          ccAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATAATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGTTCTGAATTCTGCCTCTCgtgt                                                     <  4:213346/148‑1 (MQ=255)
                                                                                          ccAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATAATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGTTCTGAATTCTGCCTCTCgtgt                                                     <  4:218439/148‑1 (MQ=255)
                                                                                           cAAAAATAACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGAATCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTAACATGGGATATCTCCAATGCACTGAATTCTGCCGCTCgtgt                                                     >  2:88124/1‑147 (MQ=255)
                                                                                            aaaaTTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCACGTCACATGGGATATCGCCAATGTTCTGAATTCTGCCTCTCgtgt                                                     >  7:177826/1‑146 (MQ=255)
                                                                                                        aTTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGTTCTGAATTCTGCCTCTCGTGTCAGCGTGCCGCTAt                                       <  6:273551/148‑1 (MQ=255)
                                                                                                            aTAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGTTCTGAATTCTGCCTCTCGTGTCAGCGTGCCGCTAttttt                                   <  5:144702/148‑1 (MQ=255)
                                                                                                             tAACCTGCGAGGGATAGGGTTTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGTTCTGAATTCTGCCTCTCGTGTCAGCGTGCCGCTATTTTTc                                  <  6:29382/148‑1 (MQ=255)
                                                                                                                    cGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGTTCTGAATTCTGCCTCTCGTGTCAGCGTGCCGCTATTTTTCATGAtt                            <  1:13872/147‑1 (MQ=255)
                                                                                                                     gAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGTTCTGAATTCTGCCTCTCGTGTCAGCGTGCCGCTATTTTTCATGATTTc                          >  8:68078/1‑148 (MQ=255)
                                                                                                                          aTAGCGTGTTGAATGGTGGTGATGATTCACACCACTACCTGGTATGTGCCCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCAAATGTTCTGAATTCTGCCTATCGTGTCAGCGTGACGATATTTTTCATGAGTTCCGGc                      >  4:264/1‑147 (MQ=255)
                                                                                                                                   ttAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGTTCTGAATTCTGCCTCTCGTGTCAGCGTGCCGCTATTTTTCATGATTTCCGGCTATCTCtttt            >  2:90495/1‑148 (MQ=255)
                                                                                                                                          tggtgATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGTTCTGAATTCTGCCTCTCGTGTCAGCGTGCCGCTATTTTTCATGATTTCCGGCTATCTCTTTTTTGGCg      >  7:280687/1‑147 (MQ=255)
                                                                                                                                           ggtgATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGTTCTGAATTCTGCctct                                                          >  6:302477/1‑94 (MQ=255)
                                                                                                                                           ggtgATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGTTCTGAATTCTGCctct                                                          <  5:302477/94‑1 (MQ=255)
                                                                                                                                           ggtgATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATAGCCAATGTTCTGAATTCTGCCTCTCGTGTCAGCGTGCCGCTATTTTTCATGATTTCCGGCTATCTCTTTTTTGGCGaa    >  7:79405/1‑148 (MQ=255)
                                                                                                                                               atgatTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGTTCTGAATTCTGCCTCTCGTGTCAGCGTGCCGCTATTTTTCATGATTTCCGGCTATCTCTTTTTTGGCGAAca  >  1:269379/1‑145 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TAGTGATAAAACTCAGATTATTCCAGTGCGCCAACCAAATGATTTATCGGTCAGACTTTTTTTAACCAGATACAGAGACACAGATGCAGCCCAAAATTTACTGGATTGATAACCTGCGAGGGATAGCGTGTTTAATGGTGGTGATGATTCACACCACTACCTGGTATGTGACCAATGCTCATAGTGTTAGCCCCGTCACATGGGATATCGCCAATGTTCTGAATTCTGCCTCTCGTGTCAGCGTGCCGCTATTTTTCATGATTTCCGGCTATCTCTTTTTTGGCGAACG  >  NC_000913/3725804‑3726092

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