Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,748,260 A→C 5.3% M240R (ATG→AGG) ‡ etfA ← electron transfer flavoprotein subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,748,2600AC 5.3% 574.9 / 14.9 169M240R (ATG→AGG) ‡etfAelectron transfer flavoprotein subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (87/73);  new base C (3/6);  total (90/79)
Fisher's exact test for biased strand distribution p-value = 3.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 170 total aligned reads displayed.

CACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAACGACGATTTTGGCAGCGGTCAG  >  NC_002947/4748112‑4748408
                                                                                                                                                    |                                                                                                                                                    
cACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCAt                                                                                                                                                     <  1:1958153/150‑1 (MQ=255)
  ccTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGt                                                                                                                                                   <  1:1071215/150‑1 (MQ=255)
   cTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGt                                                                                                                                                   <  1:3583293/149‑1 (MQ=255)
    tGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCg                                                                                                                                                 >  1:3784700/1‑150 (MQ=255)
    tGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCg                                                                                                                                                 <  2:4601495/150‑1 (MQ=255)
    tGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCg                                                                                                                                                 <  2:1032413/150‑1 (MQ=255)
     ggAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGt                                                                                                                                                >  1:3955249/1‑150 (MQ=255)
       aaGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTc                                                                                                                                              <  1:4383730/150‑1 (MQ=255)
        aGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGATGCCGACTGCGATGTAAAGTTCAGGTGCAACCACCTTGCCGGTCTGGCCAACCTGCCAGTCGtttg                                                                                                                                             <  1:1631842/150‑3 (MQ=255)
        aGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCg                                                                                                                                             <  2:3616970/150‑1 (MQ=255)
          aTCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGATGCCGACTGCGATGTAAAGTTCAGGTGCAACCACCTTGCCGGTCTGGCCAACCTGCCAGTCGTTTGGc                                                                                                                                           <  1:4417077/150‑1 (MQ=255)
          aTCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGATGCCGACTGCGATGTAAAGTTCAGGTGCAACCACCTTGCCGGTCTGGCCAACCTGCCAGTCGTTTGGc                                                                                                                                           <  1:2098830/150‑1 (MQ=255)
             ggCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACg                                                                                                                                        >  2:5069046/1‑150 (MQ=255)
              gtgcTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGATGCCGACTGCGATGTAAAGTTCAGGTGCAACCACCTTGCCGGTCTGGCCAACCTGCCAGTCGTTTGGcgcat                                                                                                                                       >  1:1869545/3‑146 (MQ=255)
                gcTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGCCCGTCTGGCCGACCTGCCAGTCGTTCGGcgcatag                                                                                                                                     >  1:3863607/1‑144 (MQ=255)
                   tcttcGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCt                                                                                                                                  <  2:300917/150‑1 (MQ=255)
                   tcttcGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCt                                                                                                                                  <  1:5306879/150‑1 (MQ=255)
                       cGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGt                                                                                                                              <  2:1010043/150‑1 (MQ=255)
                       cGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGt                                                                                                                              <  2:1352046/150‑1 (MQ=255)
                         tCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCg                                                                                                                            >  1:71535/1‑150 (MQ=255)
                            ttgttgATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACt                                                                                                                         >  1:5138427/1‑150 (MQ=255)
                                 gATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTgcggc                                                                                                                    <  2:2155149/150‑1 (MQ=255)
                                  aTCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTgcggcg                                                                                                                   >  2:430398/1‑150 (MQ=255)
                                      ccACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGcgcgcg                                                                                                               >  2:36203/1‑150 (MQ=255)
                                              aCTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCa                                                                                                       >  2:566594/1‑150 (MQ=255)
                                                      gTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAgc                                                                                               >  2:151631/1‑150 (MQ=255)
                                                       tCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAgcg                                                                                              >  2:2714475/1‑150 (MQ=255)
                                                         tttCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTc                                                                                                                                              >  1:1269755/1‑100 (MQ=255)
                                                         tttCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGcc                                                                                            >  1:3172547/1‑150 (MQ=255)
                                                         tttCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGcc                                                                                            >  1:1593941/1‑150 (MQ=255)
                                                             aTGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGc                                                                                        >  2:3952120/1‑150 (MQ=255)
                                                              tGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGTGCGAAGCACCAACCGCAGCGCCGAGCt                                                                                       >  1:3784768/1‑150 (MQ=255)
                                                              tGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCt                                                                                       >  2:2379397/1‑150 (MQ=255)
                                                              tGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCt                                                                                       >  2:2652984/1‑150 (MQ=255)
                                                                ccGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTg                                                                                     <  2:768457/150‑1 (MQ=255)
                                                                ccGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTg                                                                                     >  2:3276055/1‑150 (MQ=255)
                                                                  ggCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCTCGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTc                                                                                   <  2:2210442/150‑1 (MQ=255)
                                                                  ggCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTc                                                                                   <  1:2139218/150‑1 (MQ=255)
                                                                  ggCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTc                                                                                   >  1:4449344/1‑150 (MQ=255)
                                                                  ggCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTAGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTc                                                                                   <  2:2018949/150‑1 (MQ=255)
                                                                    ccAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCgg                                                                                 <  2:4204478/150‑1 (MQ=255)
                                                                     cAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGc                                                                                <  2:5240751/150‑1 (MQ=255)
                                                                       ggTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCa                                                                              <  2:3373091/150‑1 (MQ=255)
                                                                        gTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAg                                                                             >  1:1560700/1‑150 (MQ=255)
                                                                        gTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAg                                                                             >  1:1560710/1‑150 (MQ=255)
                                                                            tGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCTGCCAGGCTg                                                                         >  1:2741729/1‑150 (MQ=255)
                                                                              gATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTa                                                                       >  1:2649467/1‑150 (MQ=255)
                                                                              gATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTa                                                                       >  1:468656/1‑150 (MQ=255)
                                                                              gATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTa                                                                       >  2:1193096/1‑150 (MQ=255)
                                                                              gATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTa                                                                       >  1:989814/1‑150 (MQ=255)
                                                                              gATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTa                                                                       >  1:2793920/1‑150 (MQ=255)
                                                                               aTCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTAc                                                                      <  2:4995495/150‑1 (MQ=255)
                                                                                tCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACa                                                                     <  1:2757719/150‑1 (MQ=255)
                                                                                 cgcgCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAg                                                                    <  1:2208929/150‑1 (MQ=255)
                                                                                     ccGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGc                                                                <  1:4747843/150‑1 (MQ=255)
                                                                                        gagaTACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTg                                                             <  1:1756416/150‑1 (MQ=255)
                                                                                        gagaTACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTg                                                             >  1:4519443/1‑150 (MQ=255)
                                                                                        gagaTACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTg                                                             <  2:4755409/150‑1 (MQ=255)
                                                                                          gaTACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGaa                                                           <  2:1593941/150‑1 (MQ=255)
                                                                                             aCCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGtt                                                        <  2:3500744/150‑1 (MQ=255)
                                                                                             aCCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGtt                                                        <  1:1361449/150‑1 (MQ=255)
                                                                                                gACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTc                                                     >  2:895511/1‑150 (MQ=255)
                                                                                                   ggCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCAcc                                                  >  1:5163734/1‑150 (MQ=255)
                                                                                                    gCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCg                                                 <  1:1841680/150‑1 (MQ=255)
                                                                                                           aCAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCa                                          >  1:2517584/1‑150 (MQ=255)
                                                                                                             aGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATg                                        >  1:3952817/1‑150 (MQ=255)
                                                                                                             aGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATg                                        >  1:3577420/1‑150 (MQ=255)
                                                                                                             aGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATg                                        >  2:5001667/1‑150 (MQ=255)
                                                                                                              gCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGc                                       >  2:1662964/1‑150 (MQ=255)
                                                                                                              gCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGc                                       <  2:1560710/150‑1 (MQ=255)
                                                                                                                 gtgGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCAc                                    <  1:70870/150‑1 (MQ=255)
                                                                                                                   ggCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACgg                                  <  1:3670292/150‑1 (MQ=255)
                                                                                                                     cgcATCGATCTGGCCGGTCTGGCCGTCCTGCATGTCGTTCGGCACGAAGCGTGAGTCGACTGCGGAGCGCGCCGCACCTACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCCCCGTTACCCATGCCACGgcc                                <  1:5202325/150‑1 (MQ=255)
                                                                                                                     cgcAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCtt                                                              >  1:100696/1‑120 (MQ=255)
                                                                                                                      gcAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGgccg                               >  2:195827/1‑150 (MQ=255)
                                                                                                                        aaCGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGgccgcc                             >  2:2261951/1‑150 (MQ=255)
                                                                                                                             tCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGaaa                        <  1:4614605/150‑1 (MQ=255)
                                                                                                                             tCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGaaa                        >  1:4110488/1‑150 (MQ=255)
                                                                                                                              cTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAac                       >  2:1839220/1‑150 (MQ=255)
                                                                                                                               ttGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAacg                      >  1:4705793/1‑150 (MQ=255)
                                                                                                                               ttGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAacg                      >  2:3013983/1‑150 (MQ=255)
                                                                                                                               atGCAGGTCGGCCAGACCGGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAacg                      >  2:1522699/2‑150 (MQ=255)
                                                                                                                                tGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAacga                     <  1:3995233/150‑1 (MQ=255)
                                                                                                                                tGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAacga                     <  2:4387220/150‑1 (MQ=255)
                                                                                                                                tGCAGGTCGGCCAGACCGGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAacga                     >  2:2899690/1‑150 (MQ=255)
                                                                                                                                   cGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAAcgacga                  <  2:2943323/150‑1 (MQ=255)
                                                                                                                                    ggTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAACGACGAt                 >  1:1122394/1‑150 (MQ=255)
                                                                                                                                        tGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAACGACGATTTTg             >  2:4611105/1‑150 (MQ=255)
                                                                                                                                        tGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAACGACGATTTTg             <  1:1535569/150‑1 (MQ=255)
                                                                                                                                        tGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAACGACGATTTTg             <  1:4611105/150‑1 (MQ=255)
                                                                                                                                        tGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAACGACGATTTTg             <  1:4326553/150‑1 (MQ=255)
                                                                                                                                        tGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAACGACGATTTTg             <  1:3710698/150‑1 (MQ=255)
                                                                                                                                        tGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAACGACGATTTTg             >  1:4026452/1‑150 (MQ=255)
                                                                                                                                        tGGCCGACCTGCATGTCGTTCGGCACGAAGCCGGCGTCGACTGCGGCGCGAGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAACGACGATTTTg             >  1:173477/1‑150 (MQ=255)
                                                                                                                                         ggCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAACGACGAGTTTgg            >  2:3526713/1‑150 (MQ=255)
                                                                                                                                              aCCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAACGACGATTTTGGCAGCg       <  1:3013983/150‑1 (MQ=255)
                                                                                                                                              aCCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAACGACGATTTTGGCAGCg       <  1:3013629/150‑1 (MQ=255)
                                                                                                                                                cTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAACGACGATTTTGGCAGCGGt     >  1:494483/1‑150 (MQ=255)
                                                                                                                                                  gCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAACGACGATTTTGGCAGCGGTCa   >  2:4478470/1‑150 (MQ=255)
                                                                                                                                                   cATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAACGACGATTTTGGCAGCGGTCAg  >  2:2456097/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAACGACGATTTTGGCAGCGGTCAG  >  NC_002947/4748112‑4748408

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