Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,504,197 G→A 5.2% Q97* (CAG→TAG)  fryA ← putative PTS enzyme: Hpr, enzyme I and IIA components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,504,1970GA 5.2% 524.2 / 12.4 155Q97* (CAG→TAG) fryAputative PTS enzyme: Hpr, enzyme I and IIA components
Reads supporting (aligned to +/- strand):  ref base G (75/72);  new base A (4/4);  total (79/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

Only 100 of 162 total aligned reads displayed.

TTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCT  >  NC_000913/2504053‑2504344
                                                                                                                                                |                                                                                                                                                   
   cccGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGc                                                                                                                                              >  3:377042/1‑149 (MQ=255)
    ccGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAACTCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCt                                                                                                                                             >  1:1110784/1‑149 (MQ=255)
    ccGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCAGCCTACGCAGGca                                                                                                                                             >  4:256664/1‑148 (MQ=255)
       aTAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGaa                                                                                                                                          <  1:631209/149‑1 (MQ=255)
       aTAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCAGCCTACGCAGGCAGaa                                                                                                                                          <  3:256664/149‑1 (MQ=255)
        tAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAc                                                                                                                                         <  4:1168875/149‑1 (MQ=255)
         aaCTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCAGCCTACGCAGGCAGAAct                                                                                                                                        >  1:153612/1‑148 (MQ=255)
          aCTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCAGCCTACGCAGGCAGAActg                                                                                                                                       <  3:1352040/149‑3 (MQ=255)
            tGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCAGCCTACGCAGGCAGAActgac                                                                                                                                     <  2:153612/149‑5 (MQ=255)
             gTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAcgctgt                                                                                                                                    >  3:823391/1‑149 (MQ=255)
              tCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAcgctgtc                                                                                                                                   >  4:192383/1‑149 (MQ=255)
              tCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAcgctgtc                                                                                                                                   >  1:250166/1‑149 (MQ=255)
                   gCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAcgctgtcgctgt                                                                                                                              <  3:1284716/149‑1 (MQ=255)
                    cTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACgctgtcgctgtc                                                                                                                             <  1:887845/149‑1 (MQ=255)
                       tCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAcgctgtc                                                                                                                                   >  1:315964/1‑140 (MQ=255)
                       tCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAcgctgtc                                                                                                                                   <  2:315964/140‑1 (MQ=255)
                       tCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGAt                                                                                                                          >  4:476592/1‑149 (MQ=255)
                        cGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATa                                                                                                                         >  1:64311/1‑149 (MQ=255)
                          cTTTGTAACAGGGTCAGCGTACCCACGCCGAGGCCGCTTGCCAGCACATGGCCGTAAAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTATGCTTAACGCTGTCGCTGTCGATaaa                                                                                                                       <  4:1216107/149‑1 (MQ=255)
                           tttGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAg                                                                                                                      >  4:426050/1‑149 (MQ=255)
                           tttGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAg                                                                                                                      >  4:402052/1‑149 (MQ=255)
                           tttGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAg                                                                                                                      >  4:353353/1‑149 (MQ=255)
                           tttGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCCCCGTCAGTTCTGCCTGCGTAGGCTGAACGCAGTCGCTGTCGATAAAg                                                                                                                      >  3:1244405/1‑149 (MQ=255)
                              gTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAgcgc                                                                                                                   <  1:174624/149‑1 (MQ=255)
                                aaCAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCAc                                                                                                                 >  1:745300/1‑149 (MQ=255)
                                 aCAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCAcc                                                                                                                >  4:513344/1‑149 (MQ=255)
                                 aCAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCAcc                                                                                                                >  4:266036/1‑149 (MQ=255)
                                  cAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCt                                                                                                               >  2:833270/1‑149 (MQ=255)
                                    gggTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTgg                                                                                                             <  3:1414099/149‑1 (MQ=255)
                                      gTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGAt                                                                                                           <  2:932727/149‑1 (MQ=255)
                                      gTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGAt                                                                                                           <  2:250166/149‑1 (MQ=255)
                                      gTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGAt                                                                                                           >  3:809691/1‑149 (MQ=255)
                                         aGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTa                                                                                                        >  3:794732/1‑149 (MQ=255)
                                          gCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTAc                                                                                                       <  3:266036/149‑1 (MQ=255)
                                            gTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTActc                                                                                                     <  2:1030668/149‑1 (MQ=255)
                                             tACCCACGCCGACGCCGCTTGCCAGCCCATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTActct                                                                                                    >  1:841490/1‑149 (MQ=255)
                                              aCCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATTAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCtt                                                                                                   <  2:323142/149‑1 (MQ=255)
                                               cccACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTc                                                                                                  <  2:67673/149‑1 (MQ=255)
                                                ccACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGAt                                                                                                           >  3:956015/1‑139 (MQ=255)
                                                ccACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGAt                                                                                                           <  4:956015/139‑1 (MQ=255)
                                                ccACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTcc                                                                                                 >  3:669266/1‑149 (MQ=255)
                                                ccACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTcc                                                                                                 >  2:459877/1‑149 (MQ=255)
                                                 cACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCa                                                                                                >  3:902229/1‑149 (MQ=255)
                                                 cACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCa                                                                                                <  1:1307534/149‑1 (MQ=255)
                                                  acgccgacgccgCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGTGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAg                                                                                               >  1:822110/1‑149 (MQ=255)
                                                  acgccgacgccgCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCAGAACGCTGTCGCTGTCGATAAAGCGCa                                                                                                                  >  1:370896/1‑130 (MQ=255)
                                                   cgccgacgccgCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGc                                                                                              >  4:650913/1‑149 (MQ=255)
                                                      cgacgctgCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCAcg                                                                                           >  2:1047395/1‑149 (MQ=255)
                                                        acgccgCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCAcgcg                                                                                         >  1:822330/1‑149 (MQ=255)
                                                               ttGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGAt                                                                                                           >  3:480086/1‑124 (MQ=255)
                                                               ttGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACgcgccgcgcc                                                                                  >  4:704374/1‑149 (MQ=255)
                                                                   cAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGct                                                                              >  2:301396/1‑149 (MQ=255)
                                                                   cAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGct                                                                              >  1:1294308/1‑149 (MQ=255)
                                                                   cAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGACTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGct                                                                              >  2:776222/1‑149 (MQ=255)
                                                                     gCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGctct                                                                            >  2:727800/1‑149 (MQ=255)
                                                                       acaTTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTc                                                                          <  3:179486/149‑1 (MQ=255)
                                                                         aTTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCAt                                                                        <  3:704374/149‑1 (MQ=255)
                                                                              cGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTg                                                                   <  1:1411329/149‑1 (MQ=255)
                                                                               gTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCt                                                                                                               >  3:628339/1‑104 (MQ=255)
                                                                               gTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCt                                                                                                               <  4:628339/104‑1 (MQ=255)
                                                                               gTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGc                                                                  <  2:442573/149‑1 (MQ=255)
                                                                                tACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGa                                                                                                            >  3:1434568/1‑106 (MQ=255)
                                                                                tACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGa                                                                                                            <  4:1434568/106‑1 (MQ=255)
                                                                                 aCAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCg                                                                >  3:83780/1‑149 (MQ=255)
                                                                                    gTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACgctgtcgctgtc                                                                                                                             <  3:1345725/85‑1 (MQ=255)
                                                                                    gTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACgctgtcgctgtc                                                                                                                             >  4:1345725/1‑85 (MQ=255)
                                                                                     tAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTaa                                                            >  4:536109/1‑149 (MQ=255)
                                                                                     tAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTaa                                                            >  2:1296382/1‑149 (MQ=255)
                                                                                     tAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTaa                                                            >  3:1220967/1‑149 (MQ=255)
                                                                                     tAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTaa                                                            <  4:1394168/149‑1 (MQ=255)
                                                                                        atgcGGGTTAGGGCGGCTTAATGAACGCGGCAGCGGATGCGCTGTCAGTTCTGCCTGCGTAGGCTTAACGCTGTCGCTGTCGATAGATCGCACGTGGATGTACTGTTCCAGCACGCGCGGCGCCTGCTCTTCATCGGTGCCGCTAATGt                                                         <  4:276500/146‑1 (MQ=255)
                                                                                         tCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGtt                                                        >  1:904263/1‑149 (MQ=255)
                                                                                         tCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGtt                                                        >  3:1379436/1‑149 (MQ=255)
                                                                                         tCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTACCGCTAATGtt                                                        <  3:1144687/149‑1 (MQ=255)
                                                                                           cGGGTTAAGGTGGCTTTATGAACGCGTCAGCGGATGCGCCGTCAGTTCTGCCTGCGTATGCTGAACGCTGTCGCTGTCTATAAAGCGCACCTAGATGTCCACTTCCAGCACGCGCCGCGCCTTCTCGGCATCGCTGCCGCTAACGTTCa                                                      <  2:803442/149‑1 (MQ=255)
                                                                                            gggTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAAt                                                           <  3:1027288/143‑1 (MQ=255)
                                                                                                   ggcggcTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCaa                                              >  2:1100822/1‑149 (MQ=255)
                                                                                                   ggcggcTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCaa                                              >  1:1181361/1‑149 (MQ=255)
                                                                                                      ggcTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTg                                           <  4:823391/149‑1 (MQ=255)
                                                                                                         ttAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCa                                        >  1:144874/1‑149 (MQ=255)
                                                                                                          tAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCAt                                       <  2:495507/149‑1 (MQ=255)
                                                                                                          tAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCAt                                       <  4:667071/149‑1 (MQ=255)
                                                                                                          tAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCAt                                       >  1:1062448/1‑149 (MQ=255)
                                                                                                          tAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCAt                                       <  3:353353/149‑1 (MQ=255)
                                                                                                              gAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTaaa                                   <  2:841490/149‑1 (MQ=255)
                                                                                                              gAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTaaa                                   <  2:491936/149‑1 (MQ=255)
                                                                                                                  gcgGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAgg                               >  4:440646/1‑149 (MQ=255)
                                                                                                                   cgGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAggg                              <  3:139077/149‑1 (MQ=255)
                                                                                                                        gCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGggtgccg                         <  3:426050/149‑1 (MQ=255)
                                                                                                                          ggATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGggtgccggt                       <  1:1067116/149‑1 (MQ=255)
                                                                                                                              gcgcCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCa                   >  3:592307/1‑149 (MQ=255)
                                                                                                                               cgcCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCaa                  <  4:794732/149‑1 (MQ=255)
                                                                                                                                gcCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAg                                >  1:228110/1‑134 (MQ=255)
                                                                                                                                gcCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAg                                <  2:228110/134‑1 (MQ=255)
                                                                                                                                gcCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAAt                 >  3:955586/1‑149 (MQ=255)
                                                                                                                                gcCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCAGTAAATAGGGTGCCGGTGCCAAt                 >  3:894911/1‑149 (MQ=255)
                                                                                                                                       gtcctgcctgcGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACATGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGcc          >  4:691533/4‑149 (MQ=255)
                                                                                                                                        ttCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCa         <  4:809691/149‑1 (MQ=255)
                                                                                                                                              ctgcGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGc   >  1:134455/1‑149 (MQ=255)
                                                                                                                                               tgcGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCt  >  2:1290710/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCT  >  NC_000913/2504053‑2504344

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: