Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,198,992 G→A 5.3% G434G (GGG→GGA gspE → protein secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,198,9920GA 5.3% 231.9 / 2.3 75G434G (GGG→GGAgspEprotein secretion protein
Reads supporting (aligned to +/- strand):  ref base G (37/34);  new base A (2/2);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

CCCGCTACAGGCTGCCGAGCGCGAGCTATGGCCTGAGCTTGCGGCGTTGGGCTGCTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCGGTGAGACCAGCCTGGCCGAAGTGCTGCGCGTGGTGCAG  >  NC_002947/1198848‑1199136
                                                                                                                                                |                                                                                                                                                
cccGCCACAGGCTGCCGAGCGCGAGCTATGGCCTGAGCTTGCGGCGTTGGGCTGCTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCa                                                                                                                                              >  7:291362/1‑149 (MQ=255)
 ccGCTACAGGCTGCCGAGCGCGAGCTATGGCCTGAGCTTGCAGCGTTGGGCTGCTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTcgt                                                                                                                                             >  8:364083/1‑147 (MQ=255)
  cGCTACAGGCTGCCGAGCGCGAGCTATGGCCTGAGCTTGCGGCGTTGGGCTGCTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATc                                                                                                                                            <  4:263313/149‑1 (MQ=255)
   gCTACAGGCTGCCGAGCGCGAGCTATGGTCTGAGGTTGCGGCGTTGGGCTTCTCTTATCACGCCGTTTGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTTCGAATTCATCGAGCTGGACGCTGGGCTCATCg                                                                                                                                           <  1:172120/149‑1 (MQ=255)
              ccGAGCGCGAGCTATGGCCTGAGCTTGCGGCGTTGGGCTGCTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTac                                                                                                                                >  1:70514/1‑149 (MQ=255)
                     cgAGCTATGGCCTGAGCTTGCGGCGTTGGGCTGCTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGgcg                                                                                                                         >  5:345360/1‑149 (MQ=255)
                      gAGCTATGGCCTGAGCTTGCGGCGTTGGGCTGCTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGgcgc                                                                                                                        >  3:404494/1‑149 (MQ=255)
                        gCTATGGCCTGTGCGTGCGGCGTTGGGCTGGTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCg                                                                                                                                           <  6:62510/128‑1 (MQ=255)
                        gCTATGGCCTGAGCTTGCGGCGTTGGGCTGCTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCg                                                                                                                                           >  5:62510/1‑128 (MQ=255)
                           aTGGCCTGAGCTTGCGGCGTTGGGCTGCTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTTATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGt                                                                                                                    >  3:373684/1‑148 (MQ=255)
                            tGGCCTGAGCTTGCGGCGTTGGGCTGCTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGt                                                                                                                    <  2:96156/147‑1 (MQ=255)
                            tGGCCTGAGCTTGCAGCGTTGGGCTGCTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCGTCGGCCTGTTGTACGACGGCGCCAGt                                                                                                                    <  7:364083/147‑1 (MQ=255)
                                cTGAGCTTGCGGCGTTGGGCTGCTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTg                                                                                                              <  6:345360/149‑1 (MQ=255)
                                 tGAGCTTGCGGCGTTGGGCTGCTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGACTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTgg                                                                                                             >  6:25640/1‑149 (MQ=255)
                                 tGAGCTTGCGGCGTTGGGCTGCTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGACTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTgg                                                                                                             >  6:124539/1‑149 (MQ=255)
                                   aGCTTGCGGCGTTGGGCTGCTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTgg                                                                                                             <  7:16600/147‑1 (MQ=255)
                                         cggcgTTGGGCTGCTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCa                                                                                                      <  5:335896/148‑1 (MQ=255)
                                          ggcgTTGGGCTGCTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGAGTGGGGTTGTACGAATTCATCGAGCTGGACGGTGGGCTCATCGGGCTGTTGTACGCCGGCGCCCGTGAGCTGGCCATGGAgg                                                                                                    >  8:390898/1‑149 (MQ=255)
                                          ggcgTTGGGCTGCTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAg                                                                                                     >  5:215442/1‑148 (MQ=255)
                                                                 ccGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCa                                                                              >  8:373108/1‑148 (MQ=255)
                                                                  cGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCaaa                                                                            <  7:352582/149‑1 (MQ=255)
                                                                    tGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGc                                                                          >  6:375726/1‑149 (MQ=255)
                                                                        ttGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACTAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTTTCTGGCTGACCGGCGCCAAAGCCtgg                                                                      >  3:109270/1‑149 (MQ=255)
                                                                           cGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTATGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCtggtgg                                                                   >  8:82514/1‑149 (MQ=255)
                                                                           cGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCtggtgg                                                                   >  8:129148/1‑149 (MQ=255)
                                                                            gAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCCCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGc                                                                  >  2:326021/1‑149 (MQ=255)
                                                                              gCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGATCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCAATGAAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGa                                                                >  7:306004/1‑149 (MQ=255)
                                                                               cAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGAt                                                               <  6:215442/149‑1 (MQ=255)
                                                                                    ccAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGc                                                            >  4:170702/1‑147 (MQ=255)
                                                                                    ccAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCa                                                          <  1:411685/149‑1 (MQ=255)
                                                                                     cAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGgcgc                                                                                                                        >  7:270338/1‑86 (MQ=255)
                                                                                     cAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGgcgc                                                                                                                        <  8:270338/86‑1 (MQ=255)
                                                                                     cAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAg                                                         <  7:373108/149‑1 (MQ=255)
                                                                                      aGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGc                                                        >  4:106087/1‑149 (MQ=255)
                                                                                       gggCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGt                                                                                                                    >  1:103721/1‑88 (MQ=255)
                                                                                       gggCAGCGGTTATATCGGGCGACGGGGGTTGTACGAATTCTTCGCGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGt                                                                                                                    <  2:103721/88‑1 (MQ=255)
                                                                                        ggCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTTATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGa                                                      <  4:373684/149‑1 (MQ=255)
                                                                                         gCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGATGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAATCCtgg                                                                      <  5:289630/132‑1 (MQ=255)
                                                                                         gCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGATGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAATCAtgg                                                                      >  6:289630/1‑132 (MQ=255)
                                                                                         gCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGc                                                            <  2:245735/142‑1 (MQ=255)
                                                                                         gCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGc                                                            >  1:245735/1‑142 (MQ=255)
                                                                                            gCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGc                                                  <  8:291362/149‑1 (MQ=255)
                                                                                              ggTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCt                                                <  4:404494/149‑1 (MQ=255)
                                                                                              ggTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCt                                                >  6:71502/1‑149 (MQ=255)
                                                                                                  atatCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCtggtgg                                                                   >  5:55606/1‑126 (MQ=255)
                                                                                                  atatCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCtggtgg                                                                   <  6:55606/126‑1 (MQ=255)
                                                                                                     tCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGcc                                                                               <  6:17328/111‑1 (MQ=255)
                                                                                                     tCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGcc                                                                               >  5:17328/1‑111 (MQ=255)
                                                                                                       gggCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAg                                                         >  7:362966/1‑131 (MQ=255)
                                                                                                       gggCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAg                                                         <  8:362966/131‑1 (MQ=255)
                                                                                                       gggCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCgg                                       >  6:201770/1‑149 (MQ=255)
                                                                                                       gggCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCgg                                       >  7:164/1‑149 (MQ=255)
                                                                                                              tGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCGGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCGGTGAGAc                                 >  1:187008/1‑148 (MQ=255)
                                                                                                               ggggTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCGGTGAGAc                                 >  1:45359/1‑147 (MQ=255)
                                                                                                                  gTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCGGTGAGACCa                               >  8:312765/1‑146 (MQ=255)
                                                                                                                  gTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCGGTGAGACCa                               <  7:312765/146‑1 (MQ=255)
                                                                                                                     gTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCGGTGAGACCAGCCTgg                         >  5:72586/1‑149 (MQ=255)
                                                                                                                      tACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCa                                                                                                      <  8:38275/71‑1 (MQ=255)
                                                                                                                      tACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCa                                                                                                      >  7:38275/1‑71 (MQ=255)
                                                                                                                      tACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCGGTGAGACCAGCCTgg                         <  8:164/148‑1 (MQ=255)
                                                                                                                        cGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAg                                                                                                                 <  8:71246/58‑1 (MQ=255)
                                                                                                                        cGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAg                                                                                                                 >  7:71246/1‑58 (MQ=255)
                                                                                                                           aTTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCGGTGAGACCAGCCTGGCCGaa                    <  2:45359/148‑1 (MQ=255)
                                                                                                                            ttCATCGAGCTGGACGCTGGACTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCGGTGAGACCAGCCTGGCCGAAGt                  <  5:124539/149‑1 (MQ=255)
                                                                                                                            ttCATCGAGCTGGACGCTGGACTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCGGTGAGACCAGCCTGGCCGAAGt                  <  5:25640/149‑1 (MQ=255)
                                                                                                                                 cGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCGGTGAGAc                                 >  6:52399/1‑129 (MQ=255)
                                                                                                                                 cGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCGGTGAGAc                                 <  5:52399/129‑1 (MQ=255)
                                                                                                                                  gAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCGGTGAGACCAGCCTGGCCGAAGtgctgc             <  7:129148/148‑1 (MQ=255)
                                                                                                                                  gAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCGGTGAGACCAGCCTGGCCGAAGtgctgc             <  6:72586/148‑1 (MQ=255)
                                                                                                                                     cTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCgg                                       <  6:108801/119‑1 (MQ=255)
                                                                                                                                     cTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCgg                                       >  5:108801/1‑119 (MQ=255)
                                                                                                                                       ggACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTgg                                              >  2:163013/1‑110 (MQ=255)
                                                                                                                                       ggACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTgg                                              <  1:163013/110‑1 (MQ=255)
                                                                                                                                        gACGCTGGGCCCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCGGTGAGACCAGCCTGGCCGAAGTGCTGCGCgtg        <  1:326021/147‑1 (MQ=255)
                                                                                                                                            cTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCGGTGAGACCAGCCTGGCCGAAGTGCTGCGCGTGGTGCAg  <  3:106087/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CCCGCTACAGGCTGCCGAGCGCGAGCTATGGCCTGAGCTTGCGGCGTTGGGCTGCTCTTATCACGCCGTGGGTTGCGAGCAGTGCCAGGGCAGCGGTTATATCGGGCGACTGGGGTTGTACGAATTCATCGAGCTGGACGCTGGGCTCATCGGCCTGTTGTACGACGGCGCCAGTGAGCTGGCCATGCAGGACTATCTGGCTGACCGGCGCCAAAGCCTGGTGGCGATGGCCAGCGATTGCCTGGCTCGCGGTGAGACCAGCCTGGCCGAAGTGCTGCGCGTGGTGCAG  >  NC_002947/1198848‑1199136

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: