Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,009,577 G→A 5.5% Q634* (CAG→TAG)  PP_1792 ← group 2 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,009,5770GA 5.5% 212.4 / 2.2 72Q634* (CAG→TAG) PP_1792group 2 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (38/30);  new base A (2/2);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

ACTTGCGTAGCAGCAGGCACGCCGCCGTTACCGCGCTCTGCCTCCTGGTAACCTGGTTCATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCGTTGTTGACTAGGCCGATCAGCTCTCCCGTGGCAGCGCTC  >  NC_002947/2009430‑2009725
                                                                                                                                                   |                                                                                                                                                    
aCTTGCGTAGCAGCAGGCACGCCGCCGTTACCGCGCTCTGCCTCCTGGTAACCTGGTTCATCCCAAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGa                                                                                                                                                    <  2:1455613/150‑1 (MQ=255)
  ttGCGTAGCAGCAGGCACGCCGCCGTTACCGCGCTCTGCCTCCTGGTAACCTGGTTCATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGCCGCAAACAACTACGCCACCATGCTGGAcc                                                                                                                                                  <  1:258835/150‑1 (MQ=255)
  ttGCGTAGCAGCAGGCACGCCGCCGTTACCGCGCTCCGCACCATGGTAACCTGGTTCATCCCCAGGTACCCGACATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGAcc                                                                                                                                                  <  1:1068380/150‑1 (MQ=255)
     cGTAGCAGCAGGCACGCCGCCGTTACCGCGCTCTGCCTCCTGGTAACCTGGTTCATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCGGGACCATc                                                                                                                                               <  2:78876/150‑1 (MQ=255)
         gcagcaGGCACGCCGCCGTCACCGCGCTCTGCCTACTGGTAACCTGGTTCATGCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGt                                                                                                                                           >  2:222289/1‑150 (MQ=255)
           agcagGCACGCCGCCGTTACCGCGCTCTGCCTCCTGGTAACCTGGTTCATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTg                                                                                                                                         <  1:279940/150‑1 (MQ=255)
            gcagGCACGCCGCCGTTACCGCGCTCTGCCTCCTGGTAACCTGGTTCATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTgg                                                                                                                                        >  1:977501/1‑150 (MQ=255)
            gcagGCACGCCGCCGTTACCGCGCTCTGCCTCCTGGTAACCTGGTTCATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTgg                                                                                                                                        <  2:1276895/150‑1 (MQ=255)
               ggCACGCCGCCGTTACCGCGCTCTGCCTCCTGGTAACCTGGTTCATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGcc                                                                                                                                     >  2:712367/1‑150 (MQ=255)
                  aCGCCGCCGTTACCGCGCTCTGCCTCCTGGTAACCTGGTTCATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCcaca                                                                                                                                  <  1:821787/150‑1 (MQ=255)
                        ccgTTACCGCGCTCTGCCTCCTGGTAACCTGGTTCATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACt                                                                                                                            >  2:1128187/1‑150 (MQ=255)
                               cgcnCTCTGCCTCCTGGTAACCTGGTTCATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCAcc                                                                                                                     >  2:1204148/1‑150 (MQ=255)
                               cgcgCTCTGCCTCCTGGTAACCTGGTTCATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCTGACGCCCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCAcc                                                                                                                     >  1:1059089/1‑150 (MQ=255)
                               cgcgCTCTGCCTCCTGGTAACCTGGTTCATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCAcc                                                                                                                     >  2:1678205/1‑150 (MQ=255)
                                   ctctGCCTCCTGGTAACCTGGTTCATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACt                                                                                                                 <  2:56813/150‑1 (MQ=255)
                                         ctcctGGTAACCTGGTTCATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCg                                                                                                           >  1:1916719/1‑150 (MQ=255)
                                          tcctGGTAACCTGGTTCATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGAACATCCTGTTGGGCCACAGCAACTTAGCACCCACTGCGCCGa                                                                                                          >  2:587128/1‑150 (MQ=255)
                                                 aaCCTGGTTCATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCggg                                                                                                   <  1:1653915/150‑1 (MQ=255)
                                                   ccTGGTTCATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGgc                                                                                                 <  2:1108844/150‑1 (MQ=255)
                                                     tGGTTCATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGgcgc                                                                                               >  1:226695/1‑150 (MQ=255)
                                                          cATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGc                                                                                          <  2:1266471/150‑1 (MQ=255)
                                                            tCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCt                                                                                        >  1:534646/1‑150 (MQ=255)
                                                            tCCCCAGGTACCCGGCATCGCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCt                                                                                        <  2:195887/150‑1 (MQ=255)
                                                               ccAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCg                                                                                     >  1:1142648/1‑150 (MQ=255)
                                                                   gTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAAc                                                                                 <  2:226695/150‑1 (MQ=255)
                                                                       ccGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCAtt                                                                             >  2:1086452/1‑150 (MQ=255)
                                                                       ccGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCAtt                                                                             >  2:894847/1‑150 (MQ=255)
                                                                         ggCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTc                                                                           >  2:1885247/1‑150 (MQ=255)
                                                                          gCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCt                                                                          <  1:1128187/150‑1 (MQ=255)
                                                                              ctcGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGCTGGGCCACAGCAACTTGGCACGCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTaa                                                                      <  2:977501/150‑1 (MQ=255)
                                                                              ctcGTTGCTCCAACGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTaa                                                                      <  1:381497/150‑1 (MQ=255)
                                                                               tcGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGcc                                                                                                                                     >  2:1406053/1‑86 (MQ=255)
                                                                               tcGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGcc                                                                                                                                     <  1:1406053/86‑1 (MQ=255)
                                                                                 gTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGAATCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAagcc                                                                   >  1:16534/1‑150 (MQ=255)
                                                                                      tCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCg                                                              >  1:909376/1‑150 (MQ=255)
                                                                                      tCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCg                                                              >  1:531444/1‑150 (MQ=255)
                                                                                        cAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCCTCCTATTAGCCCAAAGCAACGGGGCCCGCACTTCGCGGATGCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCt                                                            >  2:1159970/1‑150 (MQ=255)
                                                                                          aGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACTACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCttt                                                          <  1:1260576/150‑1 (MQ=255)
                                                                                                tGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCgatga                                                    <  2:1639370/150‑1 (MQ=255)
                                                                                                   cAGCATGGGCCGCGAGGTCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGAttt                                                 >  1:15519/1‑150 (MQ=255)
                                                                                                   cAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGACTGGATAACCAGTTCTTTAAGCCAGCCGCTTTCGATGAttt                                                 >  2:1400345/1‑150 (MQ=255)
                                                                                                    aGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAagccag                                                                 <  1:969970/133‑1 (MQ=255)
                                                                                                    aGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAagccag                                                                 >  2:969970/1‑133 (MQ=255)
                                                                                                      cATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCCCTTTCGATGATTTCGa                                              <  2:1059089/150‑1 (MQ=255)
                                                                                                             cgcgAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCgtt                                       <  2:1142648/150‑1 (MQ=255)
                                                                                                               cgAGGCCATTGAGGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCgttgt                                     >  1:919634/1‑150 (MQ=255)
                                                                                                                gAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCgttgtt                                    <  1:712367/150‑1 (MQ=255)
                                                                                                                 aGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCGGTTGGGCCACAGCAACTTGGCACCCGCTTGGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCGttgttg                                   >  1:1076321/1‑150 (MQ=255)
                                                                                                                          aCGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCGTTGTTGACTAGGCCg                          >  1:767734/1‑150 (MQ=255)
                                                                                                                          aCGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCGTTGTTGACTAGGCCg                          >  2:507833/1‑150 (MQ=255)
                                                                                                                            gCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCGTTGTTGACTAGGCCGAt                        <  1:1086452/150‑1 (MQ=255)
                                                                                                                             ccAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCGTTGTTGACTAGGCCGATc                       >  1:1546517/1‑150 (MQ=255)
                                                                                                                             ccAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCGTTGTTGACTAGGCCGATc                       >  1:282030/1‑150 (MQ=255)
                                                                                                                             ccAACAACTACGCCACCATGCTAGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCgttg                                      >  1:282476/1‑135 (MQ=255)
                                                                                                                             ccAACAACTACGCCACCATGCTAGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCgttg                                      <  2:283004/135‑1 (MQ=255)
                                                                                                                             ccAACAACTACGCCACCATGCTAGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCgttg                                      <  2:282476/135‑1 (MQ=255)
                                                                                                                             ccAACAACTACGCCACCATGCTAGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCgttg                                      >  1:283004/1‑135 (MQ=255)
                                                                                                                                 caacTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCCCTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGCTTTCGATGTCGTTGTTGACTAGGCCGATCAGct                   >  1:864697/1‑150 (MQ=255)
                                                                                                                                  aacTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGc                                                                                          <  2:732334/78‑1 (MQ=255)
                                                                                                                                  aacTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGc                                                                                          >  1:732334/1‑78 (MQ=255)
                                                                                                                                  aacTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCGTTGTTGACTAGGCCGATCAGctc                  >  1:859941/1‑150 (MQ=255)
                                                                                                                                        gCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCGTTGTTGACTAGGCCGATCAGCTCTCCCGt            <  2:767734/150‑1 (MQ=255)
                                                                                                                                          caccaTGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCAACCACTGCGCAGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCGTTGTTGACTAGGCCGATCAGCTCTCCCGTgg          >  1:234811/1‑150 (MQ=255)
                                                                                                                                            ccaTGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCGTTGTTGACTAGGCCGATCAGCTCTCCCGTGGCa        >  1:56759/1‑150 (MQ=255)
                                                                                                                                            ccaTGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCGTTGTTGACTAGGCCGATCAGCTCTCCCGTGGCa        <  2:531444/150‑1 (MQ=255)
                                                                                                                                             caTGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCGTTGTTGACTAGGCCGATCAg                     >  2:920801/1‑136 (MQ=255)
                                                                                                                                             caTGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCGTTGTTGACTAGGCCGATCAg                     <  1:920801/136‑1 (MQ=255)
                                                                                                                                              aTGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCGTTGTTGACTAGGCCGATCAGCTCTCCCGTagcagc      <  2:534646/150‑1 (MQ=255)
                                                                                                                                              aTGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCGTTGTTGACTAGGCCGATCAGCTCTCCCGTGGCAgc      >  1:984638/1‑150 (MQ=255)
                                                                                                                                                  tGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTaa                                                                      >  2:919282/1‑82 (MQ=255)
                                                                                                                                                  tGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTaa                                                                      <  1:919282/82‑1 (MQ=255)
                                                                                                                                                  tGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCGTTGTTGACTAGGCCGATCAGCTCTCCCGTGGCAGCGCTc  >  2:334994/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACTTGCGTAGCAGCAGGCACGCCGCCGTTACCGCGCTCTGCCTCCTGGTAACCTGGTTCATCCCCAGGTACCCGGCATCTCGTTGCTCCAAGGTGTTGCCAGCATGGGCCGCGAGGCCATTGACGCCAACAACTACGCCACCATGCTGGACCATCCTGTTGGGCCACAGCAACTTGGCACCCACTGCGCCGATTCCGGGGCGCAGTGCCTGCGAAACCATTTCTTTAAGCCAGCCGCTTTCGATGATTTCGATATCGTTGTTGACTAGGCCGATCAGCTCTCCCGTGGCAGCGCTC  >  NC_002947/2009430‑2009725

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