Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,099,195 G→T 6.5% intergenic (+216/+76) PP_mr66 → / ← acsA‑I P6/acetyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,099,1950GT 6.5% 287.7 / 6.4 91intergenic (+216/+76)PP_mr66/acsA‑IP6/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base G (45/38);  new base T (3/3);  total (49/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGT  >  NC_002947/5099047‑5099342
                                                                                                                                                    |                                                                                                                                                    
cTTGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAAGCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGGATGCGTGAGGt                                                                                                                                                     >  1:1986405/1‑150 (MQ=255)
cTTGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGAAAGCGTGAGGt                                                                                                                                                     >  1:2959883/1‑150 (MQ=255)
  tGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCa                                                                                                                                                   <  1:3324769/150‑1 (MQ=255)
   gATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAg                                                                                                                                                  <  2:4456484/150‑1 (MQ=255)
   gATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAg                                                                                                                                                  <  2:2588930/150‑1 (MQ=255)
    aTGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAgg                                                                                                                                                 <  2:817586/150‑1 (MQ=255)
    aTGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAgg                                                                                                                                                 <  1:4819408/150‑1 (MQ=255)
    aTGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAgg                                                                                                                                                 <  2:3870145/150‑1 (MQ=255)
     tGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAggg                                                                                                                                                <  2:3405505/150‑1 (MQ=255)
      gCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTACCAATGCAAAAGCGGGCGTGCCCGCGCATGCGGGAGGGCCgggg                                                                                                                                               >  1:4774323/1‑150 (MQ=255)
      gCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGAATGCGGGAGGTCAgggg                                                                                                                                               >  1:4774392/1‑150 (MQ=255)
      gCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGAATGCGGGAGGTCAgggg                                                                                                                                               >  1:4774453/1‑150 (MQ=255)
        acaAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCTAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCCACAAAGCGTGAAGTCAGGGGAc                                                                                                                                             >  1:2240327/1‑150 (MQ=255)
        acaAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTc                                                                                                                                                    <  2:2263987/143‑1 (MQ=255)
        acaAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTc                                                                                                                                                    >  1:2263987/1‑143 (MQ=255)
        acaAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGAc                                                                                                                                             >  1:2523510/1‑150 (MQ=255)
        acaAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCCCGAATGCGGGAGGTCAggggcc                                                                                                                                             >  1:2240218/1‑148 (MQ=255)
         caAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACa                                                                                                                                            <  2:4694474/150‑1 (MQ=255)
         caAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACa                                                                                                                                            <  2:1044032/150‑1 (MQ=255)
          aaaaTGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTTGCCGCGAAAGCGTGAGGTCAGGGGACAg                                                                                                                                           >  1:934443/1‑150 (MQ=255)
          aaaaTGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAg                                                                                                                                           >  2:2875024/1‑150 (MQ=255)
          aaaaTGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAAAGCGTGAGGTCAGGGGACAg                                                                                                                                           >  1:934444/1‑150 (MQ=255)
          aaaaTGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTGTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAGGAGCGGGCTTGCCCGCGAATGCGTCAGTTCAggctaagc                                                                                                                                           <  2:2829238/150‑7 (MQ=38)
              tGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTa                                                                                                                                      <  1:820100/150‑1 (MQ=255)
              tGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTa                                                                                                                                      <  1:820101/150‑1 (MQ=255)
              tGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTa                                                                                                                                      <  2:2363853/150‑1 (MQ=255)
                cTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGc                                                                                                                                    <  2:2059257/150‑1 (MQ=255)
                cTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGc                                                                                                                                    <  2:1554396/150‑1 (MQ=255)
                cTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGc                                                                                                                                    <  2:4315600/150‑1 (MQ=255)
                  gACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCt                                                                                                                                  <  2:4287039/150‑1 (MQ=255)
                     tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGTCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGaa                                                                                                                               <  1:1741833/150‑1 (MQ=255)
                     tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAATGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCGTAGccggga                                                                                                                               >  2:55485/1‑146 (MQ=255)
                     tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATTCGGGAGGTCAGGGGACAGCGTAGCCTGaa                                                                                                                               >  1:3670081/1‑150 (MQ=255)
                     tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGc                                                                                                                                    >  2:1434487/1‑145 (MQ=255)
                     tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGc                                                                                                                                    <  1:1434487/145‑1 (MQ=255)
                     tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGaa                                                                                                                               <  1:3004917/150‑1 (MQ=255)
                     tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGaa                                                                                                                               <  2:5815/150‑1 (MQ=255)
                     tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGCATGCGTGAGGTGAGGGGCCAGCTTTGccggta                                                                                                                               >  1:1811391/1‑146 (MQ=255)
                     gccGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGaa                                                                                                                               <  1:1044937/149‑1 (MQ=255)
                         aTGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAATAGCGGGGGTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGTACTGa                                                                                                                           >  2:3179121/1‑150 (MQ=255)
                         aTGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGa                                                                                                                           >  2:3179033/1‑150 (MQ=255)
                          tGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGAc                                                                                                                          <  1:3661234/150‑1 (MQ=255)
                          atcACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCACTCCGGGGCCGGTGACAGAACAGTAACCGATGATTCTCGCGTAACAATGCAAAACCGGGCTTGCCCGCGAATGCGAGAGGTCAGGGGACAGCTTAGCCTGAACTGAc                                                                                                                          <  2:4601999/148‑1 (MQ=255)
                                   aCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTcgcg                                                                                                                 <  1:2152263/150‑1 (MQ=255)
                                   aCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTcgcg                                                                                                                 <  1:4468984/150‑1 (MQ=255)
                                    cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCgg                                                                                                                <  1:1484398/150‑1 (MQ=255)
                                    cGGGCCGGAATTGTGTCCTGGCGCCCCGCACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCgg                                                                                                                <  2:3290311/150‑1 (MQ=255)
                                       gCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGGATGCGTGAGGTCAGGGGACAGCTTAGCCTGGACGGACGCATTCGCgggga                                                                                                             >  1:951641/1‑148 (MQ=255)
                                        ccGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTcgcg                                                                                                                 >  2:4668791/1‑145 (MQ=255)
                                        ccGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTcgcg                                                                                                                 <  1:4668791/145‑1 (MQ=255)
                                          ggAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGc                                                                                                          <  2:4632985/150‑1 (MQ=255)
                                          ggAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGGGAGGTCAGGGGCCAGCTTAGCCTGAACGGACGCATTCGCGGGCAAGc                                                                                                          >  2:2411816/1‑150 (MQ=255)
                                          ggAATTGTTTCCTGGCGCCCCGAACCGGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGc                                                                                                          <  2:4633198/150‑1 (MQ=255)
                                                  ttCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCTTGCCCGcgaat                                                                                                  <  2:2098642/150‑5 (MQ=255)
                                                      tGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTGAGGGGGCAGCGTAGCCGGGACTGTCGCCGTCGCGGGCCAGCCCGCGCCTGCAg                                                                                              >  1:947913/1‑150 (MQ=255)
                                                        gcgcCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGAAAGCGTGAGGTGAGGGGTCAGCGGGTGCTGAACGGTCGCCTTGGCGGGCCAGCCCGCGccaggcggg                                                                                            >  1:2188816/1‑143 (MQ=255)
                                                          gcCCCGCACAGGTTCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTTCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAggggg                                                                                          >  2:1231821/1‑150 (MQ=39)
                                                          gcCCCGCACAGGTTCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGAATGCGTGAGGTCCGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCGGCAggggg                                                                                          >  2:1231822/1‑150 (MQ=255)
                                                                    tgtgCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGGATGCGTGAGGTCAGGGGACAGCGTAGCCTGAACAGACGCATTCGCGGGCAAGCCCCCTCCAGCAGGGGGTATTGCCggg                                                                                >  1:1672941/1‑150 (MQ=255)
                                                                      gggggggCCGGTCAAACAACTTAACCCCATGTTTGTGGGTTAACATTGAAAAGGGGGGGTGGCCCGCGATTGGTGGCGGGCCGGGGACTCCTTAGCGTTAACTGACCATTTCGGGGCAAAGCCCGCTCCTGCAGGGGGTATTGCCGGGaa                                                                              <  2:2644353/147‑1 (MQ=255)
                                                                             ccGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGGATGCGGGTGGTGCGGGGGCAGCGTAGCCTGGACTGTCGCCTTCGCGGGCCAGCCCGCGCCCGCCGGGGGGGTTGCCGGGGGGGGGTTa                                                                       >  1:4768419/1‑150 (MQ=255)
                                                                                 tGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGCATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGGAGGGTTTACGcc                                                                   >  2:1476189/1‑150 (MQ=255)
                                                                                     aGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCCAAATCGTGAGGTGAGGGGGCAGCTTATCCGGAAATGACGCATTAGCGGGCAAGCCCGCGCCAGCAGGGGGTATTGCCGGGGAGGGTTTAAGCCGAAg                                                               >  1:4972232/1‑150 (MQ=255)
                                                                                       aaCAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGcc                                                             >  2:3348736/1‑150 (MQ=255)
                                                                                            tAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGc                                                                                                          >  1:3010066/1‑100 (MQ=255)
                                                                                            tAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGc                                                                                                          <  2:3010066/100‑1 (MQ=255)
                                                                                                cGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGaa                                                                                                                               <  1:5116517/75‑1 (MQ=255)
                                                                                                cGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGaa                                                                                                                               >  2:5116517/1‑75 (MQ=255)
                                                                                                         tCGGGTAACAATGCAAAAGCGGGCTTGCCCCCCAAAGCGTGAGGTCAGGGGACAGCTTAGCCAGCACGGACCCATTCCCGGGCAAGCCCCCACCAGCAGGGGGTAGTGCCGGGAAGGGTTTACGCCCCAGCCCCGGTCATccccaagggg                                           >  2:2913896/1‑144 (MQ=255)
                                                                                                            ggTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCCGGTTCATCGCCTTGTGGGTa                                        >  2:2012848/1‑150 (MQ=255)
                                                                                                                       aaaaGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCgg                                                                                                                <  1:4781243/67‑1 (MQ=255)
                                                                                                                       aaaaGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCgg                                                                                                                >  2:4781243/1‑67 (MQ=255)
                                                                                                                          cgcgggcggGCCCGCGAAGGCGTGCGTTCTGGGGACCGCTTACCCTGAACTGCCGCATTCGGGGGCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTg                          <  2:3260008/149‑1 (MQ=255)
                                                                                                                           gCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGt                         >  2:2345217/1‑150 (MQ=255)
                                                                                                                                 ttGCCCGCGAATGCGTGCGGTCTGGGGACAGCTTCGCCTGAACTGCGGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACcaca                   <  2:2523510/150‑1 (MQ=255)
                                                                                                                                     ccGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCg               >  1:782352/1‑150 (MQ=255)
                                                                                                                                     ccGCGAATGCGTCAGTTCAGGCTAAGCTGTCCCCTGACCTCTCGCATTCGCGGGCAAGCCCGGTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCg               >  2:4726017/1‑150 (MQ=21)
                                                                                                                                      cgcgACTGGGTGCGGTCAGGGGACCGCTTAGCCCTACCTGACGCATTCGCGGGAAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCgg              <  2:1649924/150‑1 (MQ=255)
                                                                                                                                      cgcgAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCgg              >  1:3239937/1‑150 (MQ=255)
                                                                                                                                      cgcgAATGCGTCAGCTCAGGCTAAGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCGGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGAt                                   >  2:988047/1‑129 (MQ=18)
                                                                                                                                      cgcgAATGCGTCAGCTCAGGCTAAGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCGGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGAt                                   <  1:988047/129‑1 (MQ=18)
                                                                                                                                       gcgAATGCGTCAGTTCAGGCTAAGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCCGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCggg             >  2:3934435/1‑150 (MQ=21)
                                                                                                                                       gcgAATGCGTCAGTTCAGGCTAAGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCggg             >  1:1123130/1‑150 (MQ=21)
                                                                                                                                               gTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTAc                       >  2:1234456/1‑132 (MQ=255)
                                                                                                                                               gTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTAc                       <  1:1234456/132‑1 (MQ=255)
                                                                                                                                               gTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCaa     >  2:376503/1‑150 (MQ=255)
                                                                                                                                               gTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCaa     >  1:5019730/1‑150 (MQ=255)
                                                                                                                                               gTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCaa     >  1:308074/1‑150 (MQ=255)
                                                                                                                                               gTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAg    >  1:4103387/1‑150 (MQ=255)
                                                                                                                                                 gAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAgg   >  2:2350516/1‑150 (MQ=255)
                                                                                                                                                  aGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGt  >  1:1622696/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTTGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGT  >  NC_002947/5099047‑5099342

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