Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,745,836 T→A 18.7% intergenic (+78/+674) gltA → / ← yeiW citrate synthase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,745,8360TA18.7% 324.3 / 94.2 198intergenic (+78/+674)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (93/51);  new base A (18/18);  total (114/84)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.65e-04

Only 100 of 201 total aligned reads displayed.

GCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTG  >  NC_002947/4745691‑4745985
                                                                                                                                                 |                                                                                                                                                     
gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCATGccc                                                                                                                                                   >  1:574627/1‑150 (MQ=255)
gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCATAGCGGCCGTGgcc                                                                                                                                                   >  1:2935587/1‑150 (MQ=255)
    tCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGa                                                                                                                                               <  2:1694035/150‑1 (MQ=255)
     ccGCGCCAGCTGTCCACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGTCCGCGaa                                                                                                                                              <  2:4697508/150‑1 (MQ=255)
     ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                                                                              <  1:2887384/150‑1 (MQ=255)
     ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                              <  2:1506171/150‑1 (MQ=255)
     ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                              <  1:2779222/150‑1 (MQ=255)
     ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                              <  1:2500540/150‑1 (MQ=255)
     ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                              <  2:2356271/150‑1 (MQ=255)
      cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAag                                                                                                                                             <  1:1617942/150‑1 (MQ=255)
      cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAag                                                                                                                                             >  2:472574/1‑150 (MQ=255)
      cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  1:1035489/150‑1 (MQ=39)
       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  1:1175202/150‑1 (MQ=255)
       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  1:4067661/150‑1 (MQ=255)
       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  1:945738/150‑1 (MQ=255)
       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  1:2825216/150‑1 (MQ=255)
       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  2:4776996/150‑1 (MQ=255)
       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                            <  1:3860660/150‑1 (MQ=255)
       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                            <  1:4694013/150‑1 (MQ=255)
       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACGTCGCGGCCGTGCCCGCGAaga                                                                                                                                            >  1:3779192/1‑150 (MQ=255)
       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTATTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  1:984532/150‑1 (MQ=255)
       gcgcCAGCTGTACAACGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                            <  1:2216936/150‑1 (MQ=255)
         gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCATAGCGGCCGGGCCCCCGaatatg                                                                                                                                          >  1:2489199/1‑146 (MQ=255)
            aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGGATTGCGGCCGTGCCCCCGaacaaggcg                                                                                                                                       >  1:4228103/1‑143 (MQ=255)
             gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGGAGAGGCCgg                                                                                                                                      >  2:495394/1‑150 (MQ=255)
                 tACACCGGCGAGCATAAGCGCGACATCGTTGCCCTGAAATACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGACGGCTGCCCTCGGGCAGCCTTTTTTATTTGGCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                                                                                  <  2:5608982/150‑1 (MQ=255)
                  acacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCGTCGCGGCCGTGCCCGCGCACcggccggccccg                                                                                                                                 >  2:37853/1‑146 (MQ=255)
                     ccGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCAGAGacc                                                                                                                              >  1:4111265/1‑149 (MQ=255)
                        gCGAGCAGACGCGCCACACCGTGGCCCTCAAACACCGCACTTACGCTGCAGCAGCCGAACGCAAGACGGCGGCCCCCGGGCAGCCTTTTTTATTTGCCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCCAAGAGGCCGGCACAGAcacca                                                                                                                           <  2:3920932/150‑1 (MQ=255)
                          gAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGGCTCATAGCCGCCGTGCCCGCGAAGAGGCCGGCACCAACAAcaca                                                                                                                         >  1:661427/1‑150 (MQ=255)
                          gAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGc                                                                                                                                         >  2:3838813/1‑134 (MQ=255)
                          gAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGAcaccacc                                                                                                                         >  2:3149989/1‑149 (MQ=255)
                              aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAc                                                                                                                     >  2:5513035/1‑150 (MQ=255)
                               gAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCGCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCTTCGCGGCCTTGCCCGCGCAGAGGCCGGCCCAGACACCACACAAcc                                                                                                                    >  2:920288/1‑150 (MQ=255)
                                     cgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcgc                                                                                                                                                 >  2:5617699/1‑115 (MQ=255)
                                         aTCTTTGCCCTTAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTTCCAGCGGGCTGCATTTTTTTTTTGTCCTGTGCCGGCATCTTCGCGGCATTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtc                                                                                                          <  2:2190342/150‑1 (MQ=255)
                                           cggtGCCCTGAACGACCGCAAGTAAGCTGCATTTGCCGAACGCAAGAAGGCTGCCCTCGGGATGCCTTTTTTATTTTTCCTGTGCCGGCCTCTTCGGGGCTTTGCCGGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtctc                                                                                                        <  2:2123067/147‑1 (MQ=255)
                                            gTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCCCAAAACCTCAGTTTTTCTcc                                                                                                       <  2:878244/150‑1 (MQ=255)
                                                 ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCACAGTTTTTCTccgcccc                                                                                                  >  2:2950323/1‑149 (MQ=255)
                                                   tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCATGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCCGTTTATCAccgcccgcg                                                                                                >  2:4566644/1‑149 (MQ=255)
                                                   tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCCGCCTTGCCCCCGTAGAGGCCGGCCCAGAACCCACAAAACCACAGTTTTTATccgcccgcg                                                                                                >  2:1638344/1‑149 (MQ=255)
                                                     aaaGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTGTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACCCCCCAAAACCTCAGTTTTTGTCCGCCCGCTCa                                                                                              >  2:186060/1‑150 (MQ=255)
                                                          ccGCAAGTAAGCTGCAGTAGCCGATCGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTGTTGGCGGGCTTTGCGGGGAAGAGGCCGGCACAGGGACCACTAAAGCTCAGTTTTTGTCGGCCCGCTCTCGCAg                                                                                         >  2:5348015/1‑150 (MQ=255)
                                                           cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGc                                                                                        >  1:4627298/1‑150 (MQ=255)
                                                            gCAAGTACGCTGCATTAGCCGAACGAAAGACGGCTGCCCTCGGGCAGCATTTTTTTTTTGCCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGcc                                                                                       <  1:3230502/150‑1 (MQ=255)
                                                            gCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCTTCGCCGCCGTTCCCGCGCAGATGCCGGCACAGACACCACAAAAAAACAGATTTTATCCGCCCCCGCACGCAGcc                                                                                       >  1:2250490/1‑150 (MQ=255)
                                                                  atgctgctTTAGCCTAACGCATGATGGCTGCCCTCGGGCAGCATTTTTTTTTTGTCCTGGGCGGGCCTCTTCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                                                                 <  1:1293059/148‑1 (MQ=255)
                                                                  aaGCTGCTGTCGCCGACCGCAAGCAGGCGGCCCCCGGCCAGCCTTTTTTATTGGGCCGGGGCCGCCCCCGGCGCGGCCTTGCCCGCGACGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                                                                 <  2:4303956/150‑1 (MQ=255)
                                                                  aaGCTGCAGTAGCCTACCGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                                                                 <  2:1049363/150‑1 (MQ=255)
                                                                  aaGCGGCCGTAGCCGACCGCAAGACGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACGTCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                                                                 <  2:4650967/150‑1 (MQ=255)
                                                                   aGCGGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGGCCTGTGCCGGCCTCTTCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                                                                <  2:4918387/150‑1 (MQ=255)
                                                                    gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCGTGTGCCGGCCTCTTCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                                                               <  1:5552656/150‑1 (MQ=255)
                                                                       gCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTa                                                                            <  2:3252884/150‑1 (MQ=255)
                                                                             gCCGACCGCAAGACGGCTGCCCGCGGGCCGCCTTTTTTTTTGGTCCTGTCCGGGCCTCTACGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAAc                                                                      <  1:2089045/150‑1 (MQ=255)
                                                                               ccaccGCAAGACGGCTGCCCTCGGGCAGCCTTTTTTATTGGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACgg                                                                    <  2:4237891/146‑1 (MQ=255)
                                                                                gAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACGTCGCGGCCATTCCCGCGAAGAGGCCGGCCCACACACCACAAAACCACAGTTTTTATCCGCCCCCACACGCCCCCCCTGCACGGTATTGAACGgc                                                                   >  1:1778875/1‑150 (MQ=255)
                                                                                   ataaaaaagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCTTCGCGGCCGGGCCCGCGCAGAGGCCGGCCCCGACCCCACCAAAACTCAGTTTTTCTCCCCCCCCGCACGCAGCCCAGTCCGGGGAGTGCACGGCGCa                                                                >  1:1335055/7‑150 (MQ=255)
                                                                                    gCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCAt                                                               >  2:2268050/1‑150 (MQ=255)
                                                                                     cAAGACGGCTGCCCTCGGGATGCATTTTTTATTTGTCCTGTGCCGGCCTCTCCGCGGCCTGGCCCGCCAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc                                                              <  1:3251071/150‑1 (MQ=255)
                                                                                      aagaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCACAGTTTTTCTCCGCCCGCTCACGCAGCCCGTTCACGGTATTTAACGGCGCATcc                                                             >  1:5152610/1‑150 (MQ=255)
                                                                                       cgacgggtGCCCTGGGGCTGCCTTTTTTATTTGGCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcca                                                            <  2:159759/143‑1 (MQ=255)
                                                                                       agaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcca                                                            >  2:5559095/1‑150 (MQ=255)
                                                                                         atggcggcCCTGGGGCCGACTTTTTTTTTTGTCCTGTGCGGGCATCTGCGCGGCGTTGCCCGCGAAGCGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacc                                                          <  1:2706908/148‑1 (MQ=255)
                                                                                              tGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCCCCAAGAGGCCGGCACAGACACCAAAAAACCACAGTTTTTATACGCCCGCTCACCCACCCCGTTCAGGGTATTTAAACGCGCCACCACCAAGaa                                                     >  2:2062053/1‑150 (MQ=255)
                                                                                                cccTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCATCGCCGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTGTCCGGGTATTGAACGGCGCAGCCACCACGAACt                                                   >  2:3416277/1‑150 (MQ=255)
                                                                                                 ccTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACTTAGCGGCCTTTCCCGCGAAGAGGCCGGCACAGACACCACAAAACCGCCGTTTTTCTCCGCCCGCTCACGCAGCCCGTTCAGGGTATTTAAAGGCGCATCCACCACGAACtt                                                  >  2:522993/1‑150 (MQ=255)
                                                                                                  cTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcgc                                                                                                                                                 >  1:4368597/1‑54 (MQ=255)
                                                                                                       gCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGc                                            <  2:4096344/150‑1 (MQ=255)
                                                                                                          gatttttttttTGGGCGTGGGCCGGCCTCTGCGGGGCCTGGCCCGCGACGCGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGccac                                         <  1:198916/147‑1 (MQ=255)
                                                                                                             ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCGCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAAAGGCGCCTCCACCCCGAACGTGTTGGccaccaa                                      >  2:4280439/1‑149 (MQ=255)
                                                                                                             ttttttATTTGTCCTGTGCCGGCCTCATCCCGGCCGTGCCCGCGAAGAGGCCCGCACCGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACCCAGCCCGTTCAGGGTATTGAAAGGCGCAACCACCACGAACATGTTGGccaccaa                                      >  2:2910621/1‑149 (MQ=255)
                                                                                                                  tttttGGCCGGGGCCGGCCTCTGGGGGGCCGCCCCCGCCATGCGCCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCg                                 <  1:495394/148‑1 (MQ=255)
                                                                                                                           tgtgCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTac                        <  2:687251/150‑1 (MQ=255)
                                                                                                                                 ggCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTg                  >  2:4619400/1‑150 (MQ=255)
                                                                                                                                  gCCTCTTCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTgg                 <  2:2588723/150‑1 (MQ=255)
                                                                                                                                    ctctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGGCACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt               >  2:726464/1‑150 (MQ=39)
                                                                                                                                    ctctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGGCACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt               >  2:1356143/1‑150 (MQ=39)
                                                                                                                                    ctctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt               >  1:4622706/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCTCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTac                        <  2:2814733/140‑1 (MQ=39)
                                                                                                                                     tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc              >  1:2509188/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc              >  1:3434962/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc              >  1:3271531/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc              >  2:548807/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc              >  1:268122/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc              >  1:3678916/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc              >  2:810719/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc              >  2:3389586/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc              >  1:231215/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc              >  2:4041981/1‑150 (MQ=255)
                                                                                                                                      ctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCg             >  2:1464847/1‑150 (MQ=255)
                                                                                                                                       ttCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg            >  1:137196/1‑150 (MQ=255)
                                                                                                                                       ttCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg            >  1:5259513/1‑150 (MQ=255)
                                                                                                                                        tCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgt           <  2:3663727/150‑1 (MQ=255)
                                                                                                                                         cgcgGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCa                              <  2:3010816/130‑1 (MQ=255)
                                                                                                                                         cgcgGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCa                              >  1:3010816/1‑130 (MQ=255)
                                                                                                                                          gcggcgTTGCCCGCGAAGTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgtgt         <  1:3142462/150‑1 (MQ=255)
                                                                                                                                          gcgGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgtgt         >  2:5004252/1‑150 (MQ=255)
                                                                                                                                          gcgGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgtgt         >  1:2031995/1‑150 (MQ=255)
                                                                                                                                                ttGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCt   >  2:592431/1‑150 (MQ=255)
                                                                                                                                                 tGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTg  <  2:5127877/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTG  >  NC_002947/4745691‑4745985

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