Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,316,757 G→T 6.6% A257E (GCG→GAG)  nusA ← transcription termination/antitermination L factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,316,7570GT 6.6% 279.6 / 2.9 90A257E (GCG→GAG) nusAtranscription termination/antitermination L factor
Reads supporting (aligned to +/- strand):  ref base G (42/41);  new base T (3/4);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACG  >  NC_000913/3316610‑3316899
                                                                                                                                                   |                                                                                                                                              
tCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACAcgcgc                                                                                                                                               <  1:411932/149‑1 (MQ=255)
  cATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACgcgcgc                                                                                                                                             >  2:1396461/1‑149 (MQ=255)
  cATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGACACCCAGTTCAGTAGACACCGCCTGAACACgcgcgc                                                                                                                                             >  2:1008266/1‑149 (MQ=255)
   aTGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATAATCCCACAGGTCGATAACCATACGATCGCCACCCAGTTCAGTAGACACCGCCTCAACACGCGCGcc                                                                                                                                            >  2:310298/1‑149 (MQ=255)
    tGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCa                                                                                                                                           >  1:534267/1‑149 (MQ=255)
      gtgtgtTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACg                                                                                                                                         <  1:149536/149‑1 (MQ=255)
       tgtgtTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGc                                                                                                                                        <  2:1483656/149‑1 (MQ=255)
           tttATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATAc                                                                                                                                    <  2:591213/149‑1 (MQ=255)
           tttATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATAc                                                                                                                                    >  2:1112761/1‑149 (MQ=255)
           tttATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGAGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATAc                                                                                                                                    >  1:977686/1‑149 (MQ=255)
              aTCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTa                                                                                                                                 >  1:921066/1‑149 (MQ=255)
                cTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACg                                                                                                                               >  1:518107/1‑149 (MQ=255)
                  tCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCa                                                                                                                             >  2:632401/1‑149 (MQ=255)
                  tCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCa                                                                                                                             >  1:1112638/1‑149 (MQ=255)
                    aTCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAg                                                                                                                           <  1:772803/149‑1 (MQ=255)
                      ccaccaCGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCa                                                                                                                         <  1:1577883/149‑1 (MQ=255)
                        nccacGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCAcc                                                                                                                       >  2:66496/2‑149 (MQ=255)
                         ccacGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGGGCCGGGTTATCATCCCACAGGTCGATATCGATACGCTCGCCACCGAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCt                                                                                                                      <  1:856852/149‑1 (MQ=255)
                            cGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTAcc                                                                                                                   <  2:1573744/149‑1 (MQ=255)
                               tAGAAGCAACGTCTGCAGGTGCCATTGCGTTAATCACGAACTGAGCAGGGTTATCATCCCACAGGACGATATCGATACGCTCACAACCAACTGCAGTAGAAACCTCCTGAAAACGCTCGCCACGCATACCTCCGCATGAACCAACAGGa                                                                                                                >  1:56317/1‑149 (MQ=255)
                                aGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCTCGCCACGCATACCTACGCAAGCACCTACCGGAt                                                                                                               <  2:619156/149‑1 (MQ=255)
                                 gAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCCTACCTACGCAAGCACCTACCGGATc                                                                                                              >  2:753059/1‑149 (MQ=255)
                                  aaGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCg                                                                                                             >  2:1510515/1‑149 (MQ=255)
                                   aGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCTCGCCACGCATACCTACGCAAGCACCTACCGGATCGa                                                                                                            >  2:451143/1‑149 (MQ=255)
                                    gCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGAt                                                                                                           >  1:1507359/1‑149 (MQ=255)
                                     cAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCTCGCCACGCATACCTACGCAAGCACCTACCGGATCGATa                                                                                                          <  1:162256/149‑1 (MQ=255)
                                     cAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATa                                                                                                          <  1:1394033/149‑1 (MQ=255)
                                      aaCGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATAc                                                                                                         >  1:743885/1‑149 (MQ=255)
                                        cGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGt                                                                                                       <  2:1242053/149‑1 (MQ=255)
                                        cGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGt                                                                                                       <  2:1248198/149‑1 (MQ=255)
                                         gTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGTATCGATACGtt                                                                                                      <  2:977686/149‑1 (MQ=255)
                                           cTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTa                                                                                                    >  1:959763/1‑149 (MQ=255)
                                                 gTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTCCGCAAGCACCTACCGGATCGATACGTTTATCGTTg                                                                                              >  1:1029861/1‑149 (MQ=255)
                                                 gTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTg                                                                                              <  1:1249105/149‑1 (MQ=255)
                                                    ccATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGtt                                                                                           <  1:485707/149‑1 (MQ=255)
                                                      aTTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGtttt                                                                                         <  1:1179080/149‑1 (MQ=255)
                                                        tGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCa                                                                                       <  2:959763/149‑1 (MQ=255)
                                                              aaTCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGa                                                                                 >  2:1179315/1‑149 (MQ=255)
                                                              aaTCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGa                                                                                 >  2:226098/1‑149 (MQ=255)
                                                              aaTCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGa                                                                                 >  2:1550724/1‑149 (MQ=255)
                                                              aaTCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGa                                                                                 <  2:1314438/149‑1 (MQ=255)
                                                              aaTCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGa                                                                                 <  2:572405/149‑1 (MQ=255)
                                                               aTCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCCCGCATACCTACGCAAGCACCTCCCGGATCGATACGTGTATCGTTGGTTTTAACCGCGAt                                                                                >  1:1314913/1‑149 (MQ=255)
                                                               aTCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCGCCGCGAt                                                                                <  2:534267/149‑1 (MQ=255)
                                                                tCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGCCACCGCCTGAACACGCGCGCAACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAtt                                                                               >  2:810866/1‑149 (MQ=255)
                                                                tCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCCCCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAtt                                                                               >  2:1124426/1‑149 (MQ=255)
                                                                    gtaCTGCGCCTGGTTAGCATCCCACAGGACGATATCGATACGCTCGCAAACCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCg                                                                           <  2:278799/147‑1 (MQ=255)
                                                                    gAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCg                                                                           <  2:556698/149‑1 (MQ=255)
                                                                       cTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCAc                                                                        >  1:816570/1‑149 (MQ=255)
                                                                               ggTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAAccc                                                                <  2:518107/149‑1 (MQ=255)
                                                                                gTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACAGCCTGAACACTCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAAccca                                                               <  2:921066/149‑2 (MQ=255)
                                                                                 ttATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCgg                                                              >  1:251840/1‑149 (MQ=255)
                                                                                    tcatgaaacaGGACGATATCGATACGCTCTCCAACCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTCCCGGATCGATACGGTTATCGTTGGTGTTCACCGCGATTTTCGCACGAGAACCCGGATc                                                           <  2:151906/142‑1 (MQ=255)
                                                                                      atcCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCTACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATcgc                                                         >  2:1138673/1‑149 (MQ=255)
                                                                                       ccccACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATcgcg                                                        <  2:1393120/148‑1 (MQ=255)
                                                                                        cccACAGGTCGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGa                                                       <  1:310298/149‑1 (MQ=255)
                                                                                         ccACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAg                                                      <  2:1484755/149‑1 (MQ=255)
                                                                                             aGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTCCGCAAGCACCTACCGGATCGATACGTTTTTCGTGGGTTTTCCCCGCGATTTTCGCACGAGACCCCGGATCGCGAGCCgg                                                  >  2:320687/1‑148 (MQ=255)
                                                                                               gACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCgctg                                                >  1:906694/1‑149 (MQ=255)
                                                                                                aCGATATCGCTACGCTCGCCACCCAGTTCAGTAGACACCGCATGATCACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGCTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCgctgc                                               <  1:1296551/149‑1 (MQ=255)
                                                                                                   atatCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCttt                                            >  1:522055/1‑149 (MQ=255)
                                                                                                     atCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCAACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTaa                                          <  1:660878/149‑1 (MQ=255)
                                                                                                      tCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAAt                                         <  1:66496/149‑1 (MQ=255)
                                                                                                        gATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAAttt                                       >  2:1186055/1‑149 (MQ=255)
                                                                                                         aTACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCTCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTc                                      >  2:1401311/1‑149 (MQ=255)
                                                                                                              cTCGCCTCCCAGTTCAGTAGACACCGCCTGAACACGCTCGCCACGCATACCTACCCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTCTCGCTCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCa                                 <  2:138746/149‑1 (MQ=255)
                                                                                                              cTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCa                                 >  1:1305092/1‑149 (MQ=255)
                                                                                                               tCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTGGGGTTTCACCGCGATTTTCGAACGAGAACCCGGCTCGCGAGCCGCTGCTTTAATTTCAATCAc                                >  1:242217/1‑149 (MQ=255)
                                                                                                                 gCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCActt                              >  1:1302562/1‑149 (MQ=255)
                                                                                                                  ccACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCActtc                             <  1:632401/149‑1 (MQ=255)
                                                                                                                     cccAGTTCAGTAGACACCGCCTGCACACGCGCGCCACGCATACCTACGCAAGCGCCTACCGGATCGCTACGTTTATCGTTGGTTTTCCCCGCGATTTTCGCGCGAGACCCCGGATCGCGAGCCGCTGCTGTAATTTCAATCACttcttc                          >  1:151129/1‑149 (MQ=255)
                                                                                                                     cccAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACttcttc                          >  2:751849/1‑149 (MQ=255)
                                                                                                                       cAGTTCAGTAGACACCGCCTGAACACGCTCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGc                        <  1:451143/149‑1 (MQ=255)
                                                                                                                         gTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCg                      >  1:570101/1‑149 (MQ=255)
                                                                                                                            cAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGAtt                   >  2:663804/1‑149 (MQ=255)
                                                                                                                             aGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTATATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGCTCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGAttt                  >  1:246686/1‑149 (MQ=255)
                                                                                                                              gTAGACACCGCCTGAACACGCGCGCCACTCATACCTACGCAAGCACCTACCGGATCGATTCGTTTATCGTTGGTTTTCACCGATATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTc                 <  1:1008266/149‑1 (MQ=255)
                                                                                                                              gTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTc                 <  1:1396461/149‑1 (MQ=255)
                                                                                                                              gTAGAAAACGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTAAATAACTTCTTAGCCGATTTc                 >  2:1097670/1‑149 (MQ=255)
                                                                                                                               tcgACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATGGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCt                <  1:1080930/147‑1 (MQ=255)
                                                                                                                                aGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTg               <  1:1112761/149‑1 (MQ=255)
                                                                                                                                 gACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTgg              <  1:1051151/149‑1 (MQ=255)
                                                                                                                                  acacCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGc             >  2:983651/1‑149 (MQ=255)
                                                                                                                                    acCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAc           >  1:16428/1‑149 (MQ=255)
                                                                                                                                    acCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAc           >  2:904332/1‑149 (MQ=255)
                                                                                                                                      cGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACtt         <  2:251840/149‑1 (MQ=255)
                                                                                                                                      cGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACtt         <  2:1112638/149‑1 (MQ=255)
                                                                                                                                      cGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACtt         >  1:1379690/1‑149 (MQ=255)
                                                                                                                                         cTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGAGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCaa      <  1:1510515/149‑1 (MQ=255)
                                                                                                                                           gAACACGCGCGCCACGCATACCTCCGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCGTAAATTTCAATCACTTCTTCGCCGATTTCTGGCCCTTCAATa    >  2:773494/1‑149 (MQ=255)
                                                                                                                                             acacGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACg  <  2:1305092/149‑1 (MQ=255)
                                                                                                                                             acacGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACg  <  2:921354/149‑1 (MQ=255)
                                                                                                                                             acacGCGCGCCACGCAAACCTACCCACCCACCTCCCGTACCGATACGTTAATCGTTGGTTTTCACCGCGATTTTCGCACCACACCCCGCACCGCCCGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTAAATACg  <  2:753800/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACG  >  NC_000913/3316610‑3316899

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