Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,742,809 T→A 23.5% 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,742,8090TA23.5% 186.2 / 60.6 117intergenic (+78/+674)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (45/37);  new base A (13/14);  total (59/58)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 117 total aligned reads displayed.

CGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCAC  >  NC_002947/4742662‑4742955
                                                                                                                                                    |                                                                                                                                                  
cGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGc                                                                                                                                                  <  1:2598180/150‑1 (MQ=255)
   ccGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcg                                                                                                                                               <  2:5155174/150‑1 (MQ=255)
       ccGCGCTAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTATTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                                                                           <  2:1258577/150‑1 (MQ=255)
       ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                                                                           <  2:2030629/150‑1 (MQ=255)
       ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                                                                           <  2:1420550/150‑1 (MQ=255)
       ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                                                                           <  1:5054666/150‑1 (MQ=255)
       ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                           <  2:4472827/150‑1 (MQ=255)
       ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                           <  1:1598697/150‑1 (MQ=255)
       ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                           <  1:3619327/150‑1 (MQ=255)
         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                         <  2:4202198/150‑1 (MQ=255)
         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                         <  2:1053479/150‑1 (MQ=255)
         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                         <  2:351978/150‑1 (MQ=255)
         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                         <  1:854741/150‑1 (MQ=255)
         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                         <  2:4038966/150‑1 (MQ=255)
         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                         <  2:1053218/150‑1 (MQ=255)
         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                         <  2:4250791/150‑1 (MQ=255)
         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                         <  1:2544480/150‑1 (MQ=255)
         gcgcCAGCGGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                         <  2:5317469/150‑1 (MQ=255)
          cgcCAGCTGTACACCGGCGAGCAGCAGCGCGACATCGTTGCCCTGACAGACCGCAAGTAAGCTGCAGTAGCCGAACGCACGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGGTCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                                                                        <  2:834122/150‑1 (MQ=38)
          cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                                                                        <  1:5012289/150‑1 (MQ=255)
            ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAAGAGGc                                                                                                                                      <  1:5060325/150‑1 (MQ=255)
                cTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCGTCGCCGCCCTGCCCGCCCAGACGCCGGc                                                                                                                                  >  1:654862/1‑150 (MQ=255)
                  gTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCGTGCCCGCGCAGCGGCCGGcac                                                                                                                                >  2:2303458/1‑150 (MQ=255)
                   tACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                                                                               >  2:3144833/1‑150 (MQ=255)
                    acacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCTTCGCGGCCGTGCCCGCCCAGAGGCCCGcacac                                                                                                                              >  1:673293/1‑149 (MQ=255)
                    acacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCTTCGCGGCCGGGCCCCCGCACAGGCCCGcacac                                                                                                                              >  1:2016714/1‑149 (MQ=255)
                     cacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGGCTCTTAGCGGCCGTGCCCGCGAAGATGCCGGCACAGa                                                                                                                             >  1:3428697/1‑150 (MQ=255)
                                aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCATCGCGGCCGTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAc                                                                                                                  >  1:3593476/1‑150 (MQ=255)
                                    gcgcGACATCGTTGCCCTGAAAGACCGCAATTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCa                                                                                                              <  2:1783646/150‑1 (MQ=255)
                                       cgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGtt                                                                                                           >  1:4951349/1‑150 (MQ=255)
                                        taCATCGTGGCCCTGAAAGACCGCACGTATGCTGCTGTAGCCTAAGGCAAGAAGGGTGCCCTCGGGCAGCCTTTTTTCTTTGCCGGTGGCGGGCCTCTTGGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGttt                                                                                                          <  2:3747272/149‑1 (MQ=255)
                                          cATGGTTGCCCTTAAAGACCGAAATTAAGCTGCAGTAGCCGACCGCAACATGGCGGCCCTCGGGCAGCATTTTTTTTTTGGCTTGTGCCGGCATCTTCGGGGCCTTGCCCGCTATGAGGCCGGCACAGACACCACAAAACCTCAGttttt                                                                                                        <  1:3993607/150‑1 (MQ=255)
                                                   ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGa                                                                                                                                            <  2:1507552/105‑1 (MQ=255)
                                                   ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGa                                                                                                                                            >  1:1507552/1‑105 (MQ=255)
                                                     tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACTTCGCGGCCATGCCCGCGAAGAGGCCCGCACAGACACCACAAAAACGCAGTTTTTATccgcccgcg                                                                                             >  2:489564/1‑149 (MQ=255)
                                                            ccGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGaa                                                                                                                            >  1:4515749/1‑111 (MQ=255)
                                                            ccGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGaa                                                                                                                            <  2:4515749/112‑2 (MQ=255)
                                                             cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGc                                                                                     <  2:4485774/150‑1 (MQ=255)
                                                             cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACACACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGc                                                                                     >  2:2006948/1‑150 (MQ=255)
                                                              gCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGcc                                                                                    >  2:5447417/1‑150 (MQ=255)
                                                               cAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGccc                                                                                   >  2:4841097/1‑150 (MQ=255)
                                                                 aGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACCGACACCACAAAACCGCAGTTTTTCTCCGCCCGCGCACGCAGcccgg                                                                                 >  2:2468809/1‑148 (MQ=255)
                                                                     cgCTGCAGTAGCCGACCGCAAGACGGCGGCCCGCGGGCAGCCTTTTTTATTTGTCCTGGGCCGGCCTCTTCGCGGCCTTGCCCGCGACGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                                                             <  2:837721/149‑1 (MQ=255)
                                                                           aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                         <  1:4535488/84‑1 (MQ=25)
                                                                           aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                         >  2:4535488/1‑84 (MQ=25)
                                                                               gCCGACCGCAAGCAGGCTGCCCTCGGGCTGCCTTTTTTATTTGTCCTGGGCCGGCCTCGTGGGGGCCTGGCCCGCGACGAGGCCGGCACAGACACCACAACACCTCCGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAAc                                                                   <  2:2723176/150‑1 (MQ=255)
                                                                                ccGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                         >  1:1143734/1‑79 (MQ=255)
                                                                                ccGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                         <  2:1143734/79‑1 (MQ=255)
                                                                                ccGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACg                                                                  <  2:4951349/150‑1 (MQ=255)
                                                                                  gaccGCAAGATGGCTGCCCTCGGGCAGCCTTTTTTATTTGGCCTGGGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                                                                <  2:673293/147‑1 (MQ=255)
                                                                                   aaCGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGGGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCACCCCGTTCAGGGTATTGAACGgcg                                                               >  2:3468403/1‑150 (MQ=255)
                                                                                           aagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCGTGCCCGCGAAGAGGCCGGCACGGACACCACAAAACCTCAGTTTTTATCCGCCCGCTCACGCAGCCCGTTCAGGGTATTGAAAGGCGCATccacc                                                       >  2:593221/1‑150 (MQ=255)
                                                                                             ggCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCcaccac                                                     >  2:4519881/1‑150 (MQ=255)
                                                                                                  cccTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACt                                                >  2:2417792/1‑150 (MQ=255)
                                                                                                   ccTCGGCCAGCCTTTTTTTTTTGTCCTGTGCCGGCCTCTTGGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACtt                                               <  2:2016714/150‑1 (MQ=255)
                                                                                                     gggggCTGCTTTTTTTTTTTGCCGTGGGCCGGCCCCTTCGGGGCCGTCCCCGCTACGCGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGACCttgt                                             <  1:323146/148‑1 (MQ=255)
                                                                                                      cGGGCAGCCTTTTTTATGTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt                                            <  2:3096777/150‑1 (MQ=255)
                                                                                                             ccTTTGTTATTTGTCCTGTGCCGGCCGCGTCGCCGCCATGCCCCCGCCGAGGCCGGCACAGACACCACAAAACCACAGTTTTTAGCCCCCCCCGCACGCAGCCCGGTCACGGGAGTGAACCGCCCATCCACCACCAACGTGTTGGccccc                                     >  2:4108139/1‑150 (MQ=255)
                                                                                                               ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGGCGGCACAGGCACCACAAAAACTCAGTTTTTCTGCGCCCGCTCACGCAGCCCGGTGCGGGTATTGAACGGCGCATCCACCACGCACTTGTTGGccaccag                                   >  2:38715/1‑149 (MQ=255)
                                                                                                               ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCGTTCAGGGTATTGAACGGCGCATCCACCACGAACATGTTGGCCACCAt                                   >  2:1293414/1‑150 (MQ=255)
                                                                                                                     ttttGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCgg                             <  2:4852608/149‑1 (MQ=255)
                                                                                                                     ttttGTCCTGGGCCGGCCTCTTCGGGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCgg                             <  2:1574969/149‑1 (MQ=255)
                                                                                                                       ttGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGATGCCGGCACATACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCCCGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCa                           >  2:454765/1‑150 (MQ=255)
                                                                                                                            cTGTGCCGGCCTCTTCGCGGCCTGGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTa                      <  2:2099743/150‑1 (MQ=255)
                                                                                                                               gggCCGGCCTCTTCGCGGCCTGGCCCGCGATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCCTCGCCGGCACGCTacc                    <  1:4926020/148‑1 (MQ=255)
                                                                                                                                ggCCGGCCTCTGCGGGGCCTGGCCCGGGACGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTacca                   <  2:3898998/149‑1 (MQ=255)
                                                                                                                                  ccGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTaccacc                 >  2:1181282/1‑150 (MQ=255)
                                                                                                                                  ccGGCCTCTTCGCGGCCCTGCCCCCGAAGAGGCCGGCACATACACCACAAAACCGCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCCCCAACGCCGGCACGCTaccacc                 >  2:1570345/1‑150 (MQ=255)
                                                                                                                                   cGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCt                >  2:810861/1‑150 (MQ=255)
                                                                                                                                    ggCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTg               >  2:2312211/1‑150 (MQ=255)
                                                                                                                                       ctctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt            >  1:2654446/1‑150 (MQ=255)
                                                                                                                                       ctctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCAGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt            >  2:5371123/1‑150 (MQ=255)
                                                                                                                                        tctTCGCGGGCACGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAggg                                                                           <  1:4646473/86‑1 (MQ=37)
                                                                                                                                        tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGGCACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc           >  1:4617868/1‑150 (MQ=39)
                                                                                                                                        tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGGCACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc           >  1:4617843/1‑150 (MQ=39)
                                                                                                                                        tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                                                            >  2:2885708/1‑85 (MQ=37)
                                                                                                                                        tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                                                            <  1:2885708/85‑1 (MQ=37)
                                                                                                                                        tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                                                                >  2:940440/1‑97 (MQ=37)
                                                                                                                                        tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc           >  1:4315988/1‑150 (MQ=255)
                                                                                                                                        tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc           >  2:4766586/1‑150 (MQ=255)
                                                                                                                                        tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCgg                             <  2:934705/132‑1 (MQ=255)
                                                                                                                                        tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCgg                             >  1:934705/1‑132 (MQ=255)
                                                                                                                                        tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc           >  1:166328/1‑150 (MQ=255)
                                                                                                                                        tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc           >  2:2374533/1‑150 (MQ=255)
                                                                                                                                        tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc           >  1:4574965/1‑150 (MQ=255)
                                                                                                                                         ctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCg          <  2:210907/150‑1 (MQ=255)
                                                                                                                                         ctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCg          >  1:4859225/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg         >  2:3266741/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg         >  1:4319323/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg         >  2:1721814/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg         >  2:4963477/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg         >  1:717323/1‑150 (MQ=255)
                                                                                                                                           tCGCGGCCTTGCCCGCGAATAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgt        <  1:1266174/150‑1 (MQ=255)
                                                                                                                                           tCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACACAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgt        <  1:3144833/150‑1 (MQ=255)
                                                                                                                                           tCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgt        <  2:3016651/150‑1 (MQ=255)
                                                                                                                                           tCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgt        <  2:4371047/150‑1 (MQ=255)
                                                                                                                                           gcgcgGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgt        <  2:2366706/149‑1 (MQ=255)
                                                                                                                                             gcgGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgtgt      >  2:3857858/1‑150 (MQ=255)
                                                                                                                                               ggCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGc    >  1:689971/1‑150 (MQ=255)
                                                                                                                                                 ccTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGATTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCAc  <  2:3525443/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCAC  >  NC_002947/4742662‑4742955

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: