Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,560,380 T→A 8.3% intergenic (+27/+47) sbcB → / ← PP_1366 exodeoxyribonuclease I/transcriptional regulator MvaT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,560,3800TA 8.3% 988.4 / 52.2 320intergenic (+27/+47)sbcB/PP_1366exodeoxyribonuclease I/transcriptional regulator MvaT
Reads supporting (aligned to +/- strand):  ref base T (172/121);  new base A (10/17);  total (182/138)
Fisher's exact test for biased strand distribution p-value = 4.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 321 total aligned reads displayed.

GGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCAC  >  NC_002947/1560233‑1560525
                                                                                                                                                   |                                                                                                                                                 
gggCGAATTCTAAATGGCGCGCCAGCACGCGTGGGTGGGCGCGGACGAGGCTGGAAGGCGGGTGTTGGAGGCCTGGAAACTGAATGCTCGAAATTTGCTGTAACGTTTTTCAATCGGCTGTTGGCAGATAAACAAAAACGCCGGCAAtgt                                                                                                                                                 <  1:5911344/150‑1 (MQ=255)
  gCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGTGGGCGCCGACGAGGCTGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCATATAAACAAAAACGCCGGCAATGTGc                                                                                                                                               <  2:3120279/150‑1 (MQ=255)
    gACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCg                                                                                                                                             >  1:986396/1‑150 (MQ=255)
       ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCg                                                                                                                                          >  2:3703119/1‑150 (MQ=255)
       ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCg                                                                                                                                          >  2:5989866/1‑150 (MQ=255)
       ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCg                                                                                                                                          >  1:1965330/1‑150 (MQ=255)
       ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGAGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAATCGCCGGCAATGTGCCGGCg                                                                                                                                          >  1:3463181/1‑150 (MQ=255)
        tCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGt                                                                                                                                         <  2:6375242/150‑1 (MQ=255)
        tCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGGTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACATTGCCGGCGt                                                                                                                                         <  2:5334615/150‑1 (MQ=255)
         ctaaatGGCGCGCCAGCACGCCTGGGCGGGCGCCGACGAGGCAGGCAGGGGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATGGGCTGATGGCAGATAAAAAAACACGCCGGCAATGTGCCGGCGtt                                                                                                                                        <  1:1028518/148‑1 (MQ=255)
         ctaaaaGGCGCGCCCGCAAGCCTGGGAGGGCGCCGAGGAGGCAGGCAGGCGCGTGTTGTTGGCCTGGAAAGTGCATGCTCGACAATTGCAGTAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGtt                                                                                                                                        <  1:5374226/148‑1 (MQ=255)
         cGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGTTGGCTTATAAACAAAAACGCCGGCACATTGCCGGCGttt                                                                                                                                       <  1:1317537/150‑1 (MQ=255)
          gAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttt                                                                                                                                       >  1:5638354/1‑150 (MQ=255)
          gAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttt                                                                                                                                       >  1:6804097/1‑150 (MQ=255)
           aaaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGtttt                                                                                                                                      >  1:1805452/1‑150 (MQ=255)
           aaaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCCACAATTGCAGTAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACATTGCCGGCGtttt                                                                                                                                      <  2:5772270/150‑1 (MQ=255)
            aacggCGCGCCAGCAAGCCTGGGCGGGCGCCGACGAGGCTGGCTGGCGCGTGTTGGAGGCCTGGCAACTGAATGCTCGAAAATTGCCGGAACGCTATTCAATCGGCTGATGGAATATAAACAAAAACGCCGGCAATGTGCCGGCGttttt                                                                                                                                     <  1:4586958/147‑1 (MQ=255)
            aaaGGCGCGCCCGCACGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTGGGAGGCCTGGCAACGGCATGCTCGACAATGGCAGGAACGCTATTCAATCGGCTGATGGCAGATAACCAAAAACGCCGGCAATGTGCCGGCGttttt                                                                                                                                     <  2:7539634/150‑1 (MQ=255)
            aaaGGCGCGCCAGCACGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttttt                                                                                                                                     <  2:22257/150‑1 (MQ=255)
             aaGGCGCGCCAGCAAGCCTGGGAGGGCGCCTACGTGGCAGGCAGGCGCGTGTTGGAGGCCGGGAACCTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTttg                                                                                                                                    <  2:6719502/150‑1 (MQ=255)
             aaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACATTGCCGGCGTTTttg                                                                                                                                    <  2:8046787/150‑1 (MQ=255)
              aGGCGCCCCAGCAAGCCTGGGCGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACATTGCCGGCGTTTttgt                                                                                                                                   <  1:667678/150‑1 (MQ=255)
               ggCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTGGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACATTGCCGGCGTTTttgtt                                                                                                                                  <  2:5943477/150‑1 (MQ=255)
                gcgcgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACTGTGCCGGCGTTTttgttg                                                                                                                                 >  2:3275130/1‑150 (MQ=255)
                      ccgCACGCCTGGGCGGGCGCCGACGAGGCAGGCAGGCGCGGGTGGGCGGCCTGGCAACTGCAGGCTCGACACTTGCAGGACCGCTATTCAATCGGCGGAGGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGaa                                                                                                                           <  1:7711426/148‑1 (MQ=255)
                      cAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGaa                                                                                                                           >  1:1632373/1‑150 (MQ=255)
                      cAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGaa                                                                                                                           >  2:7227259/1‑150 (MQ=255)
                      cAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGaa                                                                                                                           <  2:6581955/150‑1 (MQ=255)
                            gCCTGGGAGGGCGCCGACGAGGCAGGCTGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATCCAATCGGCTGATGGCATATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAaga                                                                                                                     <  1:7557020/150‑1 (MQ=255)
                                gggAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATGGCAGGAACGCTATCCAATCGGCTGAGGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGc                                                                                                                 <  1:6188911/150‑1 (MQ=255)
                                   aGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACGGCATGCACGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTa                                                                                                              <  1:1887383/150‑1 (MQ=255)
                                     ggCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATg                                                                                                            >  1:8510249/1‑150 (MQ=255)
                                        gcCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCgg                                                                                                         <  1:3275130/150‑1 (MQ=255)
                                             cgaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAAt                                                                                                    <  2:2724178/150‑1 (MQ=255)
                                             cgaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAAt                                                                                                    <  2:2261657/150‑1 (MQ=255)
                                                  caggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTc                                                                                               <  1:6121641/150‑1 (MQ=255)
                                                    ggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCca                                                                                             <  2:2851703/150‑1 (MQ=255)
                                                        ggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCcagcag                                                                                         >  2:418175/1‑150 (MQ=255)
                                                         gcgcGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAgg                                                                                        >  1:577113/1‑150 (MQ=255)
                                                         gcgcGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAgg                                                                                        >  2:1134231/1‑150 (MQ=255)
                                                                   gAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGACGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGc                                                                              >  2:1343306/1‑150 (MQ=255)
                                                                    tggCCGGGCAACTGCATGCTCGAAATTTGCTGTAACGATATTCACTCGGCTGTTTCCACATAAACTACCACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCt                                                                             <  1:743768/149‑1 (MQ=255)
                                                                     ggCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCtt                                                                            >  1:2239536/1‑150 (MQ=255)
                                                                        cTGGCCACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACTGTGCCGGCGTGTTTGTTGTGCGAAGCAAGAGGGCGGATGCGGATAATCCGTCCAGCCGGGGTGCCCCGCTGTCa                                                                         >  2:6269892/1‑150 (MQ=255)
                                                                          ggCAACTGCATGCTCGACATTTGCAGGAACGCTATGCAATCGGCGGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAAc                                                                       <  1:3951350/150‑1 (MQ=255)
                                                                           gCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAAcc                                                                      >  1:1193824/1‑150 (MQ=255)
                                                                             aaCTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCAc                                                                    <  1:7014137/150‑1 (MQ=255)
                                                                             aaCTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAATCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGAAGGGTTGCCCAGCTTTCAACCaa                                                                    >  1:8055591/1‑149 (MQ=255)
                                                                               cTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAATCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGTTAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGt                                                                  >  1:8415457/1‑150 (MQ=255)
                                                                                tGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTc                                                                 <  2:5150397/150‑1 (MQ=255)
                                                                                  cATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTcgc                                                               >  1:6524278/1‑150 (MQ=255)
                                                                                     gCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTcgcctc                                                            >  2:1609758/1‑148 (MQ=255)
                                                                                     gCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTcgccgc                                                            >  1:8756723/1‑150 (MQ=255)
                                                                                     gCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTcgccgc                                                            <  1:7762770/150‑1 (MQ=255)
                                                                                       tCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGccc                                                          >  2:4630425/1‑150 (MQ=255)
                                                                                       tCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGccc                                                          >  2:1679942/1‑150 (MQ=255)
                                                                                       tCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTATTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATCATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGccc                                                          >  2:3758390/1‑150 (MQ=255)
                                                                                        cGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGcccc                                                         >  2:1435000/1‑150 (MQ=255)
                                                                                                gCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCt                                                 >  2:5515846/1‑150 (MQ=255)
                                                                                                 cAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCtt                                                <  2:1172408/150‑1 (MQ=255)
                                                                                                   ggAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTc                                              <  2:1946533/150‑1 (MQ=255)
                                                                                                   ggAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTc                                              <  2:8376712/150‑1 (MQ=255)
                                                                                                   ggAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTc                                              >  2:481166/1‑150 (MQ=255)
                                                                                                    gAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTcc                                             >  2:8258138/1‑150 (MQ=255)
                                                                                                    gAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTcc                                             >  1:4752172/1‑150 (MQ=255)
                                                                                                     aaCGCTATTCAATCGGCTGATGGCAGATAACCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCa                                            <  1:8260105/150‑1 (MQ=255)
                                                                                                        gCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCActc                                         >  2:852057/1‑150 (MQ=255)
                                                                                                          tATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCtt                                       >  1:4314667/1‑150 (MQ=255)
                                                                                                          tATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTAAACCACGTCGCCGCCCCAAGTGGATTTACACTCtt                                       >  1:2053544/1‑150 (MQ=255)
                                                                                                            ttCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTc                                     >  2:8355906/1‑150 (MQ=255)
                                                                                                              cAATCGGCTGATGGCAGAGAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAg                                   <  2:7735369/150‑1 (MQ=255)
                                                                                                              aaaTCGGCTGATGGAATATAAACAAAAACGCCGGCAATGTGCCGGCGATTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAg                                   <  1:3872302/149‑1 (MQ=255)
                                                                                                               aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTTTGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg                                  >  2:3508450/1‑150 (MQ=255)
                                                                                                                aTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAggg                                 >  2:7375614/1‑150 (MQ=255)
                                                                                                                    gCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTGTGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCtt                             >  2:794451/1‑150 (MQ=255)
                                                                                                                     cTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGATTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTtg                            <  2:6902540/150‑1 (MQ=255)
                                                                                                                      tGATGGCATATAAACAAAAACGCCGGCAATGTGCCGGCGATTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTtgt                           <  1:1553086/150‑1 (MQ=255)
                                                                                                                      tGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTtgt                           >  2:244957/1‑150 (MQ=255)
                                                                                                                       gATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTtgtg                          <  2:2464919/150‑1 (MQ=255)
                                                                                                                             aGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgcc                    >  2:4461746/1‑150 (MQ=255)
                                                                                                                              gATAAACAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgccg                   <  1:5593752/150‑1 (MQ=255)
                                                                                                                                 aaaCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTACACTCTTTGAGGGTCTTGTGGTTGCCGCCt                >  1:355546/1‑150 (MQ=255)
                                                                                                                                 aaaCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTGTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCt                >  2:4903709/1‑150 (MQ=255)
                                                                                                                                  aaCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCtt               <  2:5638289/150‑1 (MQ=255)
                                                                                                                                   aCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCttt              >  1:437091/1‑150 (MQ=255)
                                                                                                                                   aCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCttt              <  1:2788463/150‑1 (MQ=255)
                                                                                                                                    cAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTg             <  1:7202565/150‑1 (MQ=255)
                                                                                                                                     aaaaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTgg            >  1:3400867/1‑150 (MQ=255)
                                                                                                                                     aaaaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTACACACTTTCAGGGGCTTGTGGTTGCCGCCTTTgg            >  1:5999724/1‑150 (MQ=255)
                                                                                                                                       aaaCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGtt          >  1:7971559/1‑150 (MQ=255)
                                                                                                                                         aCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTTGTTTc        >  1:522500/1‑150 (MQ=255)
                                                                                                                                           gCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCACGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCGTTCCACTCTTTCAGGGTCGTGTGGGTGCCGCCTGTGGTTTCGa      >  2:395308/1‑150 (MQ=255)
                                                                                                                                            ccGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGAt     <  1:2619309/150‑1 (MQ=255)
                                                                                                                                            ccGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGAt     >  2:2964035/1‑150 (MQ=255)
                                                                                                                                              ggCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCa   <  2:7056245/150‑1 (MQ=255)
                                                                                                                                              ggCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCa   >  2:4367707/1‑150 (MQ=255)
                                                                                                                                              ggCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCa   <  1:5313639/150‑1 (MQ=255)
                                                                                                                                               gCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATTCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAAACACCTCGCCGCCCCAAGTGGCTTTACACTATTTAAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATAAc  >  1:28967/1‑150 (MQ=255)
                                                                                                                                               gCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCAAGTCGCCGCCCCAAGTGGCTTTCCACACTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCAc  >  1:2680847/1‑150 (MQ=255)
                                                                                                                                               gCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCATTCAACCACGTCCCCGCCCCAAGTGGCTTTCCACTCATTAAGGGTCTTTTGGTTGCCGCCTTTTGTTTAGATAAc  >  1:89133/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCAC  >  NC_002947/1560233‑1560525

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: