Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,198,644 A→C F333C (TTC→TGC)  glnE ← fused deadenylyltransferase/adenylyltransferase for glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,198,6440AC98.4% 677.2 / 0.4 194F333C (TTC→TGC) glnEfused deadenylyltransferase/adenylyltransferase for glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base C (78/111);  total (80/114)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01

Only 100 of 196 total aligned reads displayed.

GTTTTCAGAAAGCAGATGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGAACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACGTCGACGCACTTCACGGGCAATCATC  >  NC_000913/3198547‑3198739
                                                                                                 |                                                                                               
 ttttCAGAAAGCAGATGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGcaca                                                                                              >  2:1249354/1‑100 (MQ=255)
  tttCAGAAAGCAGATGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCacac                                                                                             >  1:2111072/1‑100 (MQ=255)
  tttCAGAAAGCAGATGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCacac                                                                                             <  2:3507128/100‑1 (MQ=255)
   ttCAGAAAGCAGATGCAGCTCGGCAATGGCGCTTAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACAcc                                                                                            <  1:3351546/100‑1 (MQ=255)
   ttCAGAAAGCAGATGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACAcc                                                                                            <  1:1250932/100‑1 (MQ=255)
   ttCAGAAAGCAGATGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACAcc                                                                                            >  1:3814117/1‑100 (MQ=255)
   ttCAGAAAGCAGATGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACAcc                                                                                            <  1:4718123/100‑1 (MQ=255)
    tCAGAAAGCAGATGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCt                                                                                           >  2:4067117/1‑100 (MQ=255)
       gAAAGCAGATGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGaa                                                                                        <  2:102280/100‑1 (MQ=255)
        aaaGCAGATGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAAc                                                                                       <  1:2978819/100‑1 (MQ=255)
        aaaGCAGATGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAAc                                                                                       <  1:2132099/100‑1 (MQ=255)
        aaaGCAGATGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAAc                                                                                       <  2:2217063/100‑1 (MQ=255)
         aaGCAGATGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACg                                                                                      >  2:3161405/1‑100 (MQ=255)
           gCAGATGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGAt                                                                                    <  2:1835746/100‑1 (MQ=255)
           gCAGATGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGAt                                                                                    >  1:4626980/1‑100 (MQ=255)
            cAGATGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATa                                                                                   >  2:3577492/1‑100 (MQ=255)
             aGATGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATaa                                                                                  >  1:4593798/1‑100 (MQ=255)
             aGATGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATaa                                                                                  <  1:1183739/100‑1 (MQ=255)
                tGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAAtt                                                                               <  2:3081111/100‑1 (MQ=255)
                tGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAAtt                                                                               <  2:765692/100‑1 (MQ=255)
                 gCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTc                                                                              >  2:5375434/1‑100 (MQ=255)
                  cAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGTGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCa                                                                             >  2:4766370/1‑100 (MQ=255)
                  cAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGAGGGTTCGCGTCCGCCGCGTATGAGCTGGCACACCTGTACGATAAATTCa                                                                             >  1:5298912/1‑100 (MQ=255)
                  cAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCa                                                                             <  2:4502610/100‑1 (MQ=255)
                  cAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCa                                                                             >  2:4401574/1‑100 (MQ=255)
                  cAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCa                                                                             <  2:4031081/100‑1 (MQ=255)
                    gCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAAt                                                                           >  2:3991332/1‑100 (MQ=255)
                    gCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAAt                                                                           >  2:3152148/1‑100 (MQ=255)
                     cTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAAtt                                                                          <  1:5158147/100‑1 (MQ=255)
                     cTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAAtt                                                                          >  1:4693837/1‑100 (MQ=255)
                      tCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAAttt                                                                         <  2:1329442/100‑1 (MQ=255)
                       cGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTc                                                                        >  1:3190175/1‑100 (MQ=255)
                       cGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTc                                                                        >  1:2332288/1‑100 (MQ=255)
                       cGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTc                                                                        >  1:4079845/1‑100 (MQ=255)
                          cAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTcgcg                                                                     <  1:990466/100‑1 (MQ=255)
                           aaTGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGa                                                                    >  2:610353/1‑100 (MQ=255)
                           aaTGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGa                                                                    <  1:5474212/100‑1 (MQ=255)
                             tGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAAt                                                                  <  2:1513420/100‑1 (MQ=255)
                             tGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAAt                                                                  >  2:2800299/1‑100 (MQ=255)
                              ggCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATg                                                                 <  1:25235/100‑1 (MQ=255)
                              ggCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATg                                                                 <  2:3391567/100‑1 (MQ=255)
                              ggCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATg                                                                 >  1:3960329/1‑100 (MQ=255)
                              ggCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATg                                                                 <  1:712692/100‑1 (MQ=255)
                                  cTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATgccgc                                                             <  1:4816182/100‑1 (MQ=255)
                                   tGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATgccgcc                                                            <  2:2894142/100‑1 (MQ=255)
                                   tGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATgccgcc                                                            <  2:1559084/100‑1 (MQ=255)
                                        gTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTgcgc                                                       >  1:653199/1‑100 (MQ=255)
                                         ttGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGcc                                                      <  1:3416040/100‑1 (MQ=255)
                                           ggCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGa                                                    >  1:2744131/1‑100 (MQ=255)
                                            gCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAg                                                   <  1:1581124/100‑1 (MQ=255)
                                             cAGTCAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGTACACCTGAACGATAAATTCAAGTTCGCGAATGCCGCCTGCGCCGAGt                                                  <  1:4985544/100‑1 (MQ=255)
                                             cAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGt                                                  <  2:2679585/100‑1 (MQ=255)
                                             cAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGt                                                  <  1:3233270/100‑1 (MQ=255)
                                              aGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGtt                                                 >  2:936007/1‑100 (MQ=255)
                                              aGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGtt                                                 <  1:4344334/100‑1 (MQ=255)
                                               gTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGttt                                                <  1:5356907/100‑1 (MQ=255)
                                               gTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGttt                                                >  1:2876327/1‑100 (MQ=255)
                                               gTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGttt                                                <  1:4394331/100‑1 (MQ=255)
                                                tAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTg                                               <  1:371247/100‑1 (MQ=255)
                                                 aaaGAGCGCGATTGCAGCGACGGTTCGCGTCCGCGGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGa                                              <  2:2082032/100‑1 (MQ=255)
                                                 aaaGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGa                                              <  1:3570202/100‑1 (MQ=255)
                                                 aaaGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGa                                              <  1:1997404/100‑1 (MQ=255)
                                                 aaaGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGa                                              >  1:1009970/1‑100 (MQ=255)
                                                 aaaGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGa                                              <  2:5092309/100‑1 (MQ=255)
                                                 aaaGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGa                                              >  2:5127686/1‑100 (MQ=255)
                                                  aaGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGat                                             <  1:598234/100‑1 (MQ=255)
                                                      gcgcGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTg                                         <  1:1292409/100‑1 (MQ=255)
                                                       cgcgATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGt                                        <  2:3750311/100‑1 (MQ=255)
                                                        gcgATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTc                                       <  2:485927/100‑1 (MQ=255)
                                                        gcgATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTc                                       <  2:2717850/100‑1 (MQ=255)
                                                        gcgATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTc                                       <  1:3028288/100‑1 (MQ=255)
                                                          gATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCgg                                     <  1:3651789/100‑1 (MQ=255)
                                                             tGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATGAATTCAAGTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCa                                  <  1:5274376/100‑1 (MQ=255)
                                                             tGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCa                                  <  2:2960203/100‑1 (MQ=255)
                                                              gCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCaa                                 >  1:4315922/1‑100 (MQ=255)
                                                              gCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCAAACCTGAACGATAAATTCCATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTaaa                                 >  2:4940939/1‑100 (MQ=255)
                                                               cAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCaaa                                <  2:1724334/100‑1 (MQ=255)
                                                                aGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAAc                               >  1:4295592/1‑100 (MQ=255)
                                                                 gCGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACAGCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAAcc                              <  2:4107709/100‑1 (MQ=255)
                                                                  cGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCa                             <  1:2618028/100‑1 (MQ=255)
                                                                  cGACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCa                             <  1:3468059/100‑1 (MQ=255)
                                                                   gACGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCAc                            <  1:3540/100‑1 (MQ=255)
                                                                    aCGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACg                           <  2:185743/100‑1 (MQ=255)
                                                                     cGGTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACGt                          <  2:2722129/100‑1 (MQ=255)
                                                                      ggTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACGTc                         >  1:3175974/1‑100 (MQ=255)
                                                                      ggTTCGCGTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACGTc                         >  1:2491903/1‑100 (MQ=255)
                                                                            cgTCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACGTCGACGCa                   <  2:3829356/100‑1 (MQ=255)
                                                                              tCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACGTCGACGCACt                 <  1:4269143/100‑1 (MQ=255)
                                                                              tCCGCCGCGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACGTCGACGCACt                 <  1:3064296/100‑1 (MQ=255)
                                                                                cgccgcGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACGTCGACGCACTTc               <  1:3610815/100‑1 (MQ=255)
                                                                                  ccgcGAATGAGCTGGCACACCTGAACTATCAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACGTCGACGCACTTCAc             <  2:5226966/100‑1 (MQ=255)
                                                                                  ccgcGAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACGTCGACGCACTTCAc             <  1:2175746/100‑1 (MQ=255)
                                                                                    gcgAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACGTCGACGCACTTCACgg           <  2:5340510/100‑1 (MQ=255)
                                                                                      gAATGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACGTCGACGCACTTCACGGGc         <  1:1867347/100‑1 (MQ=255)
                                                                                       aaTGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACGTCGACGCACTTCACGGGCa        <  1:2278945/100‑1 (MQ=255)
                                                                                        aTGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACGTCGACGCACTTCACGGGCaa       <  1:820624/100‑1 (MQ=255)
                                                                                        aTGAGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACGTCGACGCACTTCACGGGCaa       >  1:4496138/1‑100 (MQ=255)
                                                                                           aGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACGTCGACGCACTTCACGGGCCatca    >  2:4711258/1‑100 (MQ=255)
                                                                                           aGCTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACGTCGACGCACTTCACGGGCAatca    >  1:1947519/1‑100 (MQ=255)
                                                                                             cTGGCACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACGTCGACGCACTTCACGGGCAatcatc  <  1:5085215/100‑1 (MQ=255)
                                                                                                 |                                                                                               
GTTTTCAGAAAGCAGATGCAGCTCGGCAATGGCGCTGAGCGTTGGCAGTAAAGAGCGCGATTGCAGCGACGGTTCGCGTCCGCCGCGAATGAGCTGGAACACCTGAACGATAAATTCAATTTCGCGAATGCCGCCTGCGCCGAGTTTGATATTGTCGGTCAAACCACGTCGACGCACTTCACGGGCAATCATC  >  NC_000913/3198547‑3198739

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