Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 873,418 G→A 6.2% G18S (GGC→AGC)  PP_5457 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947873,4180GA 6.2% 243.5 / 4.1 80G18S (GGC→AGC) PP_5457hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (39/36);  new base A (3/2);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

CCGAGGATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTAC  >  NC_002947/873270‑873564
                                                                                                                                                    |                                                                                                                                                   
ccGAGGATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTggcag                                                                                                                                                     >  3:219026/1‑149 (MQ=255)
  gAGGATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTggcagg                                                                                                                                                    >  2:358463/1‑148 (MQ=255)
    ggATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTTTCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCAt                                                                                                                                                 >  6:124195/1‑149 (MQ=255)
     gATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATg                                                                                                                                                >  4:323735/1‑149 (MQ=255)
         gCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCtt                                                                                                                                             <  3:160046/148‑1 (MQ=255)
         gCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCtt                                                                                                                                             <  3:241319/148‑1 (MQ=255)
          ccAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGc                                                                                                                                           >  7:148801/1‑149 (MQ=255)
            aGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGCTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGAATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAAAGCTGGCAGGCATGCTTCCCa                                                                                                                                         >  3:326467/1‑149 (MQ=255)
              gtggGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAAGCATGCCTGCCAGc                                                                                                                                       <  6:201179/149‑1 (MQ=255)
                 tgCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAg                                                                                                                                    >  4:63499/1‑149 (MQ=255)
                   cTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAggg                                                                                                                                  >  4:140604/1‑149 (MQ=255)
                   cTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCAGCGAGGACGACACCTGATGCAGCGTATGTACGAACCGGATAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAgg                                                                                                                                   <  8:15436/148‑1 (MQ=255)
                     gTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGa                                                                                                                                >  4:368413/1‑149 (MQ=255)
                        aGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATc                                                                                                                              >  2:103981/1‑148 (MQ=255)
                         gCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGa                                                                                                                            >  5:227074/1‑149 (MQ=255)
                          gcTGTGGTGTGCATTCTGCAGGTGCATTGTTAGGATTGGTTGCTTGAAGTCCTGCGGGAGGTCGACTCCTGATGGTGGGTATCTACGAACCGGAAAGCGTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGAt                                                                                                                           <  4:326467/148‑1 (MQ=255)
                          ccTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGAt                                                                                                                           <  8:148801/149‑1 (MQ=255)
                             gtggtgTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGt                                                                                                                         >  3:147030/1‑148 (MQ=255)
                                gtgtGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTacac                                                                                                                     >  8:37355/1‑149 (MQ=255)
                                gtgtGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTacac                                                                                                                     >  6:194689/1‑149 (MQ=255)
                                  gtgCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGt                                                                                                                   >  2:274561/1‑149 (MQ=255)
                                   tgCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGt                                                                                                                   <  3:140604/148‑1 (MQ=255)
                                   tgCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTg                                                                                                                  >  1:374041/1‑149 (MQ=255)
                                      aTTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTc                                                                                                               >  1:233327/1‑149 (MQ=255)
                                          tGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGcc                                                                                                           >  7:117588/1‑149 (MQ=255)
                                            cAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGcc                                                                                                           <  2:360978/147‑1 (MQ=255)
                                            cAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGcc                                                                                                           <  5:15498/147‑1 (MQ=255)
                                            cAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGt                                                                                                         <  1:103981/149‑1 (MQ=255)
                                            cAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCCGGCAGGCATGCTTGCCAGCGa                                                                                                                                     >  4:98529/1‑121 (MQ=255)
                                                tGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGAcc                                                                                                     <  2:374041/149‑1 (MQ=255)
                                                 gCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCt                                                                                                    >  1:20033/1‑149 (MQ=255)
                                                    tGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGt                                                                                                >  3:295240/1‑149 (MQ=255)
                                                     ggTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGCACACGTGGTCGGCCGTGACCTTGt                                                                                                <  8:260962/148‑1 (MQ=255)
                                                     ggTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGCACACGTGGTCGGCCGTGACCTTGt                                                                                                >  7:260962/1‑148 (MQ=255)
                                                      gTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGGCGGCCGTGACCTTGTc                                                                                               <  1:358463/148‑1 (MQ=255)
                                                        gACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGt                                                                                            >  4:169632/1‑149 (MQ=255)
                                                         aCGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTg                                                                                           <  2:233327/149‑1 (MQ=255)
                                                           gAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCg                                                                                              >  4:57290/1‑144 (MQ=255)
                                                           gAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCCTGACCTTGTCg                                                                                              <  3:57290/144‑1 (MQ=255)
                                                            aaTGGATGCTTGAAGTCCGCCGCGAGGACGTCACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTgccg                                                                                        >  1:33750/1‑149 (MQ=255)
                                                                gATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCgg                                                                                    <  6:266493/149‑1 (MQ=255)
                                                                 aTGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGc                                                                                   >  6:74298/1‑149 (MQ=255)
                                                                 aTGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGc                                                                                   >  8:339028/1‑149 (MQ=255)
                                                                    cTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGa                                                                                 <  3:323735/148‑1 (MQ=255)
                                                                    cTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAg                                                                                <  3:368413/149‑1 (MQ=255)
                                                                        aaGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGc                                                                            >  5:329376/1‑149 (MQ=255)
                                                                          gTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGccc                                                                          <  3:63499/149‑1 (MQ=255)
                                                                          gTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGccc                                                                          >  3:31933/1‑149 (MQ=255)
                                                                          gTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGAGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGccc                                                                          <  2:166979/149‑1 (MQ=255)
                                                                            cagctgagAGGACTAGACATGATGCATCTTATCTAAGAAACGGACAGACTGATTGATGAGCAAATGCTGGCAGGCATGCTTGCATGAGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGACCt                                                                         <  2:256034/141‑1 (MQ=255)
                                                                              gccgcGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGc                                                                      <  6:227074/149‑1 (MQ=255)
                                                                                   gAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGc                                                                 <  2:20033/149‑1 (MQ=255)
                                                                                          acacCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTgg                                                          <  4:295240/149‑1 (MQ=255)
                                                                                          acacCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTgg                                                          <  7:339028/149‑1 (MQ=255)
                                                                                              cTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCt                                                      >  6:352991/1‑149 (MQ=255)
                                                                                              cTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCt                                                      <  4:147030/149‑1 (MQ=255)
                                                                                              cTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCAAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCt                                                      >  3:354491/1‑149 (MQ=255)
                                                                                               tGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCtc                                                     >  8:3201/1‑149 (MQ=255)
                                                                                                 aTGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCtcgt                                                   >  1:272356/1‑149 (MQ=255)
                                                                                                    cAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGCTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCtcgtcg                                                 <  1:274561/149‑1 (MQ=255)
                                                                                                    cAGCGTAGCTACGAACCGGAAAGCCTGCGTCTGGCGCAAATGCTGGCAGGCATGCGTGCCAGCGCGGGGCTCGAGGTACACGTGGTCGGCCGTGACCTTTTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTcgtcgt                                                <  7:270456/149‑1 (MQ=255)
                                                                                                      gCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTgg                                              <  7:37355/149‑1 (MQ=255)
                                                                                                         tATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCtg                                           <  4:31933/149‑1 (MQ=255)
                                                                                                             tACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCtgat                                         >  1:173842/1‑147 (MQ=255)
                                                                                                               cGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCtgat                                         <  2:33750/145‑1 (MQ=255)
                                                                                                               cGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAg                                      >  2:364118/1‑148 (MQ=255)
                                                                                                                gAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAg                                      >  7:87076/1‑147 (MQ=255)
                                                                                                                    cGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGcc                                <  1:364118/149‑1 (MQ=255)
                                                                                                                           ccTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATgcg                         >  5:53434/1‑149 (MQ=255)
                                                                                                                           ccTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATgcg                         <  5:194689/149‑1 (MQ=255)
                                                                                                                                  gAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCAGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACgg                    >  8:283628/1‑147 (MQ=255)
                                                                                                                                   aGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCa                  <  7:315384/148‑1 (MQ=255)
                                                                                                                                     gcgcAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCt               >  7:196423/1‑149 (MQ=255)
                                                                                                                                       gcAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGa             >  7:22801/1‑149 (MQ=255)
                                                                                                                                         cctcgatGGCAAGCATGCCTGCAAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGACGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGAACGGAACCTGATc           >  6:44333/7‑149 (MQ=255)
                                                                                                                                         aaaTGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATc           >  6:118707/1‑149 (MQ=255)
                                                                                                                                           tcgCTGGCAAGCATGCCTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATc           <  4:330656/145‑1 (MQ=255)
                                                                                                                                           tcgCTGGCAAGCATGCCTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATc           >  3:330656/3‑147 (MQ=255)
                                                                                                                                            tGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGAt        <  8:117588/149‑1 (MQ=255)
                                                                                                                                            cgCTGGCAAGCATGCCTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGAt        >  8:55062/2‑149 (MQ=255)
                                                                                                                                              cTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATg       <  5:352991/148‑1 (MQ=255)
                                                                                                                                              cTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGc      <  4:219026/149‑1 (MQ=255)
                                                                                                                                                 gcaggcaTGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGt    <  5:74298/148‑1 (MQ=255)
                                                                                                                                                  caggcaTGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTAc  <  8:87076/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCGAGGATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTAC  >  NC_002947/873270‑873564

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