Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,740,884 C→G 15.0% intergenic (+75/+677) gltA → / ← yeiW citrate synthase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,740,8840CG15.0% 242.0 / 26.4 92intergenic (+75/+677)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (38/40);  new base G (7/7);  total (45/47)
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

AAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGC  >  NC_002947/4740735‑4741031
                                                                                                                                                     |                                                                                                                                                   
aaGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGc                                                                                                                                                     >  1:2112164/1‑150 (MQ=255)
aaGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGc                                                                                                                                                     >  1:1310776/1‑150 (MQ=255)
  gATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCt                                                                                                                                                   >  2:592893/1‑150 (MQ=255)
   aTCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCtt                                                                                                                                                  <  2:183656/150‑1 (MQ=255)
   aTCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCtt                                                                                                                                                  >  1:2318075/1‑150 (MQ=255)
            ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                         <  1:592893/150‑1 (MQ=255)
            ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCCCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                         <  1:61113/150‑1 (MQ=255)
              gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGATCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                       <  2:200285/150‑1 (MQ=255)
              gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                       <  1:1306489/150‑1 (MQ=255)
              gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                       <  1:214083/150‑1 (MQ=255)
                 ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGc                                                                                                                                    <  1:1437791/150‑1 (MQ=255)
                   aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGcc                                                                                                                                               >  1:1383625/1‑137 (MQ=255)
                   aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGcc                                                                                                                                               <  2:1383625/137‑1 (MQ=255)
                   aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                                                                  <  1:159823/150‑1 (MQ=255)
                     cTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAAGAGGc                                                                                                                                    <  1:2082304/146‑1 (MQ=39)
                     cTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCACTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAAGAGGc                                                                                                                                    >  2:2082304/1‑146 (MQ=39)
                        tACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                                                                             <  2:1310776/150‑1 (MQ=255)
                        tACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                                                                             >  1:1895347/1‑150 (MQ=255)
                         acacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAg                                                                                                                            <  1:1094580/150‑1 (MQ=38)
                           acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGcc                                                                                                                                                    <  1:1879901/124‑1 (MQ=255)
                           acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGcc                                                                                                                                                    >  2:1879901/1‑124 (MQ=255)
                                 gAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCtt                                                                                                                                                  >  2:1597637/1‑120 (MQ=255)
                                 gAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCtt                                                                                                                                                  <  1:1597637/120‑1 (MQ=255)
                                    cAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                                                                 <  2:1895347/150‑1 (MQ=255)
                                     aGACGCGCAACCTCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAc                                                                                                                <  1:2165106/150‑1 (MQ=255)
                                      gAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCCAAGAGGCCGGCACAGACACCACAAAAcc                                                                                                               <  1:2130892/150‑1 (MQ=255)
                                       aaGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCt                                                                                                              >  1:972273/1‑150 (MQ=255)
                                         gcgcGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCa                                                                                                            <  2:2112164/150‑1 (MQ=255)
                                             gACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGca                                                                                                                               <  1:699974/127‑1 (MQ=255)
                                             gACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGca                                                                                                                               >  2:699974/1‑127 (MQ=255)
                                             gACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGttt                                                                                                        >  2:853322/1‑150 (MQ=255)
                                                   gTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTcc                                                                                                  <  1:853322/150‑1 (MQ=255)
                                                       cccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgccc                                                                                              >  2:1488457/1‑150 (MQ=255)
                                                        ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcccg                                                                                             >  1:743227/1‑150 (MQ=255)
                                                        ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAAAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccacccc                                                                                             >  1:523707/1‑149 (MQ=255)
                                                         cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcccgc                                                                                            >  1:2407784/1‑150 (MQ=255)
                                                         cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCCCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcccgc                                                                                            >  2:1019470/1‑150 (MQ=255)
                                                          tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAAGAGGCCGGcaca                                                                                                                             >  2:1311890/1‑116 (MQ=37)
                                                          tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAAGAGGCCGGcaca                                                                                                                             <  1:1311890/116‑1 (MQ=37)
                                                          tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                       <  1:2347541/106‑1 (MQ=255)
                                                          tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                       >  2:2347541/1‑106 (MQ=255)
                                                                 ccGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAACACCTCAGTTTTTCTCCGCCCGCTCACGCAg                                                                                    >  1:325644/1‑150 (MQ=255)
                                                                 ccGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAg                                                                                    >  2:2099814/1‑150 (MQ=255)
                                                                  ctcAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCCTCGCGGCCTCGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGc                                                                                   <  2:962039/148‑1 (MQ=255)
                                                                   ccAAGTAAGCTGCAGTAGCCCAACGCAACAAGCCCGCACTCGAACAGCCTTTTTCATTTGTCCTGTGCTGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAACCCTCAGTTTTTCTCCGCCCGCTCTCGCAGcc                                                                                  <  2:2076365/149‑1 (MQ=37)
                                                                    cAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGccc                                                                                 >  1:533454/1‑150 (MQ=255)
                                                                          aGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                                                           <  2:743227/150‑1 (MQ=255)
                                                                          aGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                                                           <  2:533454/150‑1 (MQ=255)
                                                                           gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                                                          <  1:1125750/150‑1 (MQ=255)
                                                                           gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGACCTACTCCCGGGCACAACCGCGAAAATTCCAGCATAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                                                          <  1:2099814/150‑1 (MQ=255)
                                                                               cAGTAGCCGAACGCAAGAAGGCTTCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCTCGGCCTTGCACGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                                                                      <  1:826456/150‑1 (MQ=255)
                                                                                aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtc                                                                                                     <  1:1789558/118‑1 (MQ=255)
                                                                                aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtc                                                                                                     >  2:592127/1‑118 (MQ=255)
                                                                                aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtc                                                                                                     <  1:592127/118‑1 (MQ=255)
                                                                                aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtc                                                                                                     >  2:1789558/1‑118 (MQ=255)
                                                                                               aaagGCTGCCCGAGGGCAGCCTTCTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccac                                                      >  2:1985593/2‑150 (MQ=255)
                                                                                                aagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                         <  2:1353434/66‑1 (MQ=255)
                                                                                                aagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                         >  1:1353434/1‑66 (MQ=255)
                                                                                                   gCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACCGACACCACAAAACCTCAGTTTTTCTCCGCCCGCCCACGCAGCCCTTTCAGGGTATTGACCGGCCCCTccaccccc                                                  >  2:2397419/1‑147 (MQ=255)
                                                                                                   gCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACg                                                  >  2:2397418/1‑150 (MQ=255)
                                                                                                   gCTGCCCGAGGGCAGCCTTCTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACg                                                  >  2:2027907/1‑150 (MQ=255)
                                                                                                     tGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCACCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGaa                                                >  1:2073021/1‑150 (MQ=255)
                                                                                                        ccTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACtt                                             >  1:2052415/1‑150 (MQ=255)
                                                                                                           cGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt                                          <  2:972273/150‑1 (MQ=255)
                                                                                                           cGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTCGACCGGTCTCAGACCTACACACACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt                                          <  1:1488457/150‑1 (MQ=255)
                                                                                                           agggCAGCCTTCTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt                                          <  1:1985593/149‑1 (MQ=255)
                                                                                                            gggCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCCCGAACttgttg                                         >  1:309111/1‑150 (MQ=255)
                                                                                                                    ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt                                          <  2:2304665/141‑1 (MQ=255)
                                                                                                                    ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt                                          >  1:2304665/1‑141 (MQ=255)
                                                                                                                    ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCAt                                 <  2:2052415/150‑1 (MQ=255)
                                                                                                                              tggcTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCTTCCACCACGAACTTGGTGGCCACCATCGCGGGCACg                       >  1:1109290/4‑150 (MQ=37)
                                                                                                                                      ccGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                                                                      >  2:1755319/1‑95 (MQ=255)
                                                                                                                                      ccGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                                                                      <  1:1755319/95‑1 (MQ=255)
                                                                                                                                      ccGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTaccacc               <  2:979526/150‑1 (MQ=255)
                                                                                                                                      ccGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTaccacc               >  1:979526/1‑150 (MQ=255)
                                                                                                                                         gCCTCTTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                                                          <  2:1968324/88‑1 (MQ=25)
                                                                                                                                         gCCTCTTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                                                          >  1:1968324/1‑88 (MQ=25)
                                                                                                                                         gCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                                                              <  2:627479/100‑1 (MQ=255)
                                                                                                                                         gCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                                                              >  1:627479/1‑100 (MQ=255)
                                                                                                                                         gCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTgg            >  1:1798421/1‑150 (MQ=255)
                                                                                                                                         gCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTgg            <  1:1019470/150‑1 (MQ=255)
                                                                                                                                          ccTCTTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTac                   <  2:1575953/142‑1 (MQ=39)
                                                                                                                                          ccTCTTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTac                   >  1:1575953/1‑142 (MQ=39)
                                                                                                                                           ctctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGCGTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCGTCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt          >  2:701740/1‑150 (MQ=39)
                                                                                                                                           ctctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgccc                                                                                              <  2:2213257/66‑1 (MQ=255)
                                                                                                                                           ctctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgccc                                                                                              >  1:2213257/1‑66 (MQ=255)
                                                                                                                                           ctctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                                                          <  1:1703877/86‑1 (MQ=255)
                                                                                                                                           ctctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                                                          >  2:1703877/1‑86 (MQ=255)
                                                                                                                                           ctctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt          >  1:201699/1‑150 (MQ=255)
                                                                                                                                            tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGcca                                     <  1:690186/122‑1 (MQ=38)
                                                                                                                                            tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGcca                                     >  2:690186/1‑122 (MQ=38)
                                                                                                                                              ttCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTac                   >  1:2031960/1‑138 (MQ=255)
                                                                                                                                              ttCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTac                   <  2:2031960/138‑1 (MQ=255)
                                                                                                                                              ttCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg       >  2:58640/1‑150 (MQ=255)
                                                                                                                                               tCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgt      <  1:2027907/150‑1 (MQ=255)
                                                                                                                                                   ggCCTTGCCCGCGAAGAGGCCGGCACAAACACCACAAAACAGAAGCTTTACGCTGCGCGCTCACGCAGCCGTTTCAGGCTAGCGGAAGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGc  >  1:550369/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGC  >  NC_002947/4740735‑4741031

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