Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,362,529 T→C 12.4% intergenic (‑176/+22) yjdQ ← / ← pheU pseudogene, P4‑like integrase remnant/tRNA‑Phe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,362,5290TC12.4% 209.7 / 21.0 81intergenic (‑176/+22)yjdQ/pheUpseudogene, P4‑like integrase remnant/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base T (32/39);  new base C (4/6);  total (36/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.61e-01

GTGCTATAGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAG  >  NC_000913/4362381‑4362676
                                                                                                                                                    |                                                                                                                                                   
gTGCTATAGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGt                                                                                                                                                     >  4:719491/1‑149 (MQ=255)
 tGCTATAGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTg                                                                                                                                                    <  3:640916/149‑1 (MQ=255)
 tGCTATAGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTg                                                                                                                                                    <  3:44156/149‑1 (MQ=255)
   cTATAGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGa                                                                                                                                                  <  1:641650/149‑1 (MQ=255)
    tataGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAg                                                                                                                                                 <  1:204091/149‑1 (MQ=255)
    tataGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAg                                                                                                                                                 <  3:612326/149‑1 (MQ=255)
    tataGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGCCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAg                                                                                                                                                 <  2:119115/149‑1 (MQ=255)
          tCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCACGGAGGTCCGt                                                                                                                                           <  4:121194/149‑1 (MQ=255)
           cAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCACGGAGGTCCGtt                                                                                                                                          <  2:232933/149‑1 (MQ=255)
           cAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTAATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGtt                                                                                                                                          <  1:452614/149‑1 (MQ=255)
            aGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCACGGAGGTCCGttt                                                                                                                                         <  2:231992/149‑1 (MQ=255)
               aaCAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTtata                                                                                                                                      <  3:131700/149‑1 (MQ=255)
               aaCAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTtata                                                                                                                                      >  4:239815/1‑149 (MQ=255)
                aCAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTatat                                                                                                                                     >  2:714323/1‑149 (MQ=255)
                 cAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCACGGAGGTCCGttt                                                                                                                                         >  3:746442/1‑144 (MQ=255)
                 cAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCACGGAGGTCCGttt                                                                                                                                         >  1:242167/1‑144 (MQ=255)
                 cAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCACGGAGGTCCGttt                                                                                                                                         <  4:746442/144‑1 (MQ=255)
                 cAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCACGGAGGTCCGttt                                                                                                                                         <  2:242167/144‑1 (MQ=255)
                  aaaTTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGa                                                                                                                                   >  3:432177/1‑149 (MQ=255)
                    aTTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAAt                                                                                                                                 >  1:682447/1‑149 (MQ=255)
                         cccGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTCTATATGAATTTGGt                                                                                                                            <  1:109163/149‑1 (MQ=255)
                                       tGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAACCGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGaa                                                                                                              <  1:343237/149‑1 (MQ=255)
                                        gCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAAt                                                                                                             <  2:764475/149‑1 (MQ=255)
                                         cccGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGt                                                                                                                                           >  4:718669/1‑118 (MQ=255)
                                         cccGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGt                                                                                                                                           <  3:718669/118‑1 (MQ=255)
                                          ccGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCg                                                                                                           >  2:192991/1‑149 (MQ=255)
                                          ccGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCg                                                                                                           >  4:40088/1‑149 (MQ=255)
                                          ccGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCg                                                                                                           >  4:3269/1‑149 (MQ=255)
                                             cAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAAc                                                                                                        >  3:207088/1‑149 (MQ=255)
                                                aTTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGTCAGAATGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGCCCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCaa                                                                                                     <  2:551333/149‑1 (MQ=255)
                                                aTTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGTCCGGACTCGGAATCGAACCaa                                                                                                     <  4:419887/149‑1 (MQ=255)
                                                aTTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCaa                                                                                                     <  3:702393/149‑1 (MQ=255)
                                                aTTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCaa                                                                                                     <  2:675877/149‑1 (MQ=255)
                                                   gCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGa                                                                                                  <  1:478987/149‑1 (MQ=255)
                                                               gggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGAtttt                                                                                      >  1:417429/1‑149 (MQ=255)
                                                                  agaATCTGGTTGAGTTCGGTAGAATTGATTTCGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATCTTCaa                                                                                   <  1:220272/149‑1 (MQ=255)
                                                                   gaATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGAttgaaaat                                                                                  >  2:249780/1‑143 (MQ=255)
                                                                   gaATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAAt                                                                                  <  1:567265/149‑1 (MQ=255)
                                                                      tCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATccc                                                                               <  2:417429/149‑1 (MQ=255)
                                                                      tCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTGAAAATccc                                                                               <  3:705880/149‑1 (MQ=255)
                                                                        tGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCt                                                                             >  3:232059/1‑149 (MQ=255)
                                                                        tGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCt                                                                             <  3:712346/149‑1 (MQ=255)
                                                                          gTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGtt                                                                                                                                          >  4:526102/1‑86 (MQ=255)
                                                                          gTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGtt                                                                                                                                          <  3:526102/86‑1 (MQ=255)
                                                                                  cGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGa                                                                   <  3:40088/149‑1 (MQ=255)
                                                                                  cGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGa                                                                   >  2:568265/1‑149 (MQ=255)
                                                                                    gTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACt                                                                 <  3:239815/149‑1 (MQ=255)
                                                                                      aGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGa                                                               <  2:86316/149‑1 (MQ=255)
                                                                                      aGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGa                                                               <  4:432177/149‑1 (MQ=255)
                                                                                      aGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGa                                                               <  4:160318/149‑1 (MQ=255)
                                                                                         aTTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTAc                                                                      <  4:443190/139‑1 (MQ=255)
                                                                                         aTTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTAc                                                                      >  3:443190/1‑139 (MQ=255)
                                                                                         aTTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCt                                                            >  4:445499/1‑149 (MQ=255)
                                                                                          ttGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTa                                                           >  1:387701/1‑149 (MQ=255)
                                                                                            gATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCTGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATc                                                         >  1:493994/1‑149 (MQ=255)
                                                                                             aTTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATcc                                                        >  3:540751/1‑149 (MQ=255)
                                                                                             aTTTGGAGGCAGAACGCTTAAATCGTGGCGTCCCGAAACGAAAAACGGACCTCAGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATcc                                                        >  4:262091/1‑149 (MQ=255)
                                                                                                     gCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGa                                                               <  4:265349/134‑1 (MQ=255)
                                                                                                     gCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGa                                                               >  3:265349/1‑134 (MQ=255)
                                                                                                      cAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACgg                                               <  2:387701/149‑1 (MQ=255)
                                                                                                        gAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAAc                                                 <  4:599782/145‑1 (MQ=255)
                                                                                                        gAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAAc                                                 >  3:599782/1‑145 (MQ=255)
                                                                                                            gCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCa                                         >  3:726458/1‑149 (MQ=255)
                                                                                                            gCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCa                                         >  2:593558/1‑149 (MQ=255)
                                                                                                             cTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCAt                                        >  1:326796/1‑149 (MQ=255)
                                                                                                              ttAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCAtt                                       >  3:50237/1‑149 (MQ=255)
                                                                                                              ttAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCAtt                                       >  2:725996/1‑149 (MQ=255)
                                                                                                                aaaTCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTaa                                     >  3:210995/1‑149 (MQ=255)
                                                                                                                           tCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGAc                                                                  >  3:501932/1‑109 (MQ=255)
                                                                                                                           tCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGAc                                                                  <  4:501932/109‑1 (MQ=255)
                                                                                                                             cTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGca                        >  4:83718/1‑149 (MQ=255)
                                                                                                                                    gAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCg                 >  4:433362/1‑149 (MQ=255)
                                                                                                                                        aaCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCaa             <  2:326796/149‑1 (MQ=255)
                                                                                                                                         aCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAAc            <  2:698708/149‑1 (MQ=255)
                                                                                                                                         aCGGACCTCCACGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAAc            >  1:638988/1‑149 (MQ=255)
                                                                                                                                           ggACCTCCACGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGa                                                               <  4:354709/96‑1 (MQ=255)
                                                                                                                                           ggACCTCCACGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGa                                                               >  3:354709/1‑96 (MQ=255)
                                                                                                                                             aCCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCAGAATCGAACCAAGGACACGGGGATTTTCCATCCCGTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTaa        <  3:433362/149‑1 (MQ=255)
                                                                                                                                              ccTCCGTGGAGGTCCGTTTATATGAGTTTGGTGCCCGGAGTCGGAATCGAACCAAGGACCCGGGGATTTTCAATCCCCTGCTCTACCGACTGCTCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAAt       >  4:611826/1‑149 (MQ=21)
                                                                                                                                                  cGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCa   <  3:440745/149‑1 (MQ=255)
                                                                                                                                                  cGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCa   <  1:410038/149‑1 (MQ=255)
                                                                                                                                                   gTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCGACCTAATTTCAg  <  2:170364/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTGCTATAGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAG  >  NC_000913/4362381‑4362676

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