Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,652,200 A→G 5.2% intergenic (+119/+42) PP_1448 → / ← PP_1449 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,652,2000AG 5.2% 230.4 / 2.1 76intergenic (+119/+42)PP_1448/PP_1449membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/39);  new base G (2/2);  total (35/41)
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.88e-01

GCGGGGGCTGCACAGGCTGGCTTGACCTGTAGCGGCCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTTCTCAACACGCCCAGGTCCAGCGCATTTTCCCCAT  >  NC_002947/1652057‑1652336
                                                                                                                                               |                                                                                                                                        
gCGGGGGCTGCACAGGCTGGCTTGACCTGTAGCGGCCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTg                                                                                                                                     >  5:266893/1‑149 (MQ=255)
  gggggCTGCACAGGCTGGCTTGACCTGTAGCGGGCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCAAAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCt                                                                                                                                   <  8:252352/149‑1 (MQ=255)
        tGCACAGGCTGGCTTGACCTGTAGCGGCCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGCGTGTCAGCAGCAGGTCACAAGATTGCTGATATc                                                                                                                             >  6:24558/1‑149 (MQ=255)
          cacaGGCTGGCTTGACCTGTAGCGGCCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCTCAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATcc                                                                                                                            <  7:367801/148‑1 (MQ=255)
           acaGGCTGGCTTGACCTGTAGCGGCCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATccc                                                                                                                           >  3:213016/1‑148 (MQ=255)
                 tggctTGACCTGTAGCGGCCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATgg                                                                                                                    >  1:380931/1‑149 (MQ=255)
                       gACCTGTAGCGGCCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATcccc                                                                                                                          >  5:78964/1‑137 (MQ=255)
                       gACCTGTAGCGGCCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATcccc                                                                                                                          <  6:78964/137‑1 (MQ=255)
                         ccTGTAGCGGCCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGcc                                                                                                            >  8:55680/1‑149 (MQ=255)
                          cTGTAGCGGCCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCcg                                                                                                           >  3:13157/1‑149 (MQ=255)
                           tGTAGCGGCCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCcgc                                                                                                          >  8:159491/1‑149 (MQ=255)
                           tGTAGCGGCCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCACAAGGGCACTGGAATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGAACAAGATTGCTGATATACCCAAAATGCCATGGCcgc                                                                                                          >  6:13004/1‑149 (MQ=255)
                            gTAGCGGCCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCcgc                                                                                                          >  6:303381/1‑148 (MQ=255)
                             tAGCGGCCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCcgcg                                                                                                         <  5:143466/148‑1 (MQ=255)
                             tAGCGGCCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCcgc                                                                                                          <  4:213016/147‑1 (MQ=255)
                             tAGCGGCAGCTTCGCGGGCAAGCCCGCCACCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCcgc                                                                                                          <  6:22604/147‑1 (MQ=255)
                               gCGGCCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATg                                                                                                      <  5:410563/149‑1 (MQ=255)
                                cGGCCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATcc                                                                                                                            >  2:346952/1‑126 (MQ=255)
                                cGGCCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATcc                                                                                                                            <  1:346952/126‑1 (MQ=255)
                                   ccTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGt                                                                                                  <  2:171154/149‑1 (MQ=255)
                                    ctctTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTc                                                                                                 <  6:294961/149‑1 (MQ=255)
                                     tctTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATgg                                                                                                     >  8:129001/1‑144 (MQ=255)
                                     tctTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATgg                                                                                                     <  7:129001/144‑1 (MQ=255)
                                     tctTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCt                                                                                                >  8:257659/1‑149 (MQ=255)
                                          gcgGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGt                                                                                            <  5:303381/148‑1 (MQ=255)
                                          gcgGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGt                                                                                            <  4:217295/148‑1 (MQ=255)
                                             ggCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGctt                                                                                        <  5:387843/149‑1 (MQ=255)
                                                aaGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGcttctt                                                                                     <  7:257659/149‑1 (MQ=255)
                                                       cccccGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGTCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGc                                                                              <  2:112370/149‑1 (MQ=255)
                                                       cccccGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGc                                                                              <  7:159491/149‑1 (MQ=255)
                                                        ccccGCAGTGAAATCACAAGCGCACGCGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCAGTAGATTGCTGATATCCCGAAATGGCCAGGGACGCGATGGCGTCTTAGTGCTTCTTCCCTTGc                                                                              <  5:13004/148‑1 (MQ=255)
                                                        ccccGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCa                                                                             <  1:201004/149‑1 (MQ=255)
                                                               gagaAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTTTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGa                                                                      <  3:246480/149‑1 (MQ=255)
                                                               gagaAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCg                                                                       >  2:3301/1‑148 (MQ=255)
                                                                     tCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCa                                                                >  3:118832/1‑149 (MQ=255)
                                                                       acaAGCCCAAGGGCTGGGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCGAGGGCCGCGCTGGCGTCTTAGTGCTTCTGCCCTTGGCTGCTCGAAACCCCtt                                                              >  8:293305/1‑149 (MQ=255)
                                                                         aaGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGc                                                            >  6:228322/1‑149 (MQ=255)
                                                                              cAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCt                                                       >  5:193600/1‑149 (MQ=255)
                                                                                aGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGc                                                     <  7:164047/149‑1 (MQ=255)
                                                                                     tgtgGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTa                                                >  1:275480/1‑149 (MQ=255)
                                                                                       tgGTAGTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTAt                                               <  2:115528/148‑1 (MQ=255)
                                                                                        ggTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGc                                             >  5:328632/1‑149 (MQ=255)
                                                                                        ggTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGc                                             >  7:294004/1‑149 (MQ=255)
                                                                                         gTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGcc                                            <  6:31936/149‑1 (MQ=255)
                                                                                            ttCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCAc                                         <  6:328632/149‑1 (MQ=255)
                                                                                              cGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACt                                        <  2:380931/148‑1 (MQ=255)
                                                                                              cGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCATGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTg                                       >  3:342095/1‑149 (MQ=255)
                                                                                               gCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAGGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAAc                                                                   <  5:313712/120‑1 (MQ=255)
                                                                                               gCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAGGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAAc                                                                   <  5:308051/120‑1 (MQ=255)
                                                                                               gCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAGGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAAc                                                                   >  6:308051/1‑120 (MQ=255)
                                                                                               gCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAGGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAAc                                                                   >  6:313712/1‑120 (MQ=255)
                                                                                               gCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTg                                       <  7:365705/148‑1 (MQ=255)
                                                                                                     gAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTtctc                                 >  8:24230/1‑148 (MQ=255)
                                                                                                     gAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTtctc                                 >  2:68910/1‑148 (MQ=255)
                                                                                                     gAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTTCTCa                                <  8:294004/149‑1 (MQ=255)
                                                                                                       aGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTTCTCAac                              >  4:348598/1‑149 (MQ=255)
                                                                                                       aGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTTCTCAac                              <  4:13157/149‑1 (MQ=255)
                                                                                                         cGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTTCTCAacac                            <  7:55680/149‑1 (MQ=255)
                                                                                                            gAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGa                                                   >  4:155430/1‑123 (MQ=255)
                                                                                                            gAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGa                                                   <  3:155430/123‑1 (MQ=255)
                                                                                                             aGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTTCTCAACACGccc                        <  6:266893/149‑1 (MQ=255)
                                                                                                               gCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGGCTATGCCCACTGGTTCTCAACACGCCCAg                      >  5:363184/1‑149 (MQ=255)
                                                                                                               gCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTTCTCAACACGCCCAg                      >  6:249759/1‑149 (MQ=255)
                                                                                                                cACGCGGCCCGGGTGTCAGCAGCAGGGCACAACATTGCTGATATCCCCAAATGCCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTGCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTAAGCCCACTGGTTCTCACCACGCCCAgg                     >  5:69456/1‑149 (MQ=255)
                                                                                                                   gcgGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTtctc                                 >  7:383178/1‑134 (MQ=255)
                                                                                                                   gcgGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTtctc                                 <  8:383178/134‑1 (MQ=255)
                                                                                                                    cgGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTTCTCAACACGCCCAGGTCCa                 <  2:275480/149‑1 (MQ=255)
                                                                                                                    cgGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTTCTCAACACGCCCAGGTCCa                 <  1:68910/149‑1 (MQ=255)
                                                                                                                    cgGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCATGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTTCTCAACACGCCCAGGTCCa                 <  4:342095/149‑1 (MQ=255)
                                                                                                                      gCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTTCTCAACACGCCCAGGTCCAgc               <  1:3301/149‑1 (MQ=255)
                                                                                                                           ggTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTTCTCAACACGCCCAGGTCCAGCGCAtt          >  3:151390/1‑149 (MQ=255)
                                                                                                                               tCAGCATCAGGGCACAAGATTGCTGATATCCGCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAAGCCTTTGCTCGCTGCGACTATGCGCACTGGTTCTCATCACGCCCTGGTCCAGCGCATTTTcc      <  5:24558/149‑1 (MQ=255)
                                                                                                                                   cagcagGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTTCTCAACACGCCCAGGTCCAGCGCATTTTCCCCAt  <  6:408374/149‑1 (MQ=255)
                                                                                                                                   cagcagGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTTCTCAACACGCCCAGGTCCAGCGCATTTTCCCCAt  >  5:408374/1‑149 (MQ=255)
                                                                                                                                            acaAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTTCTCAACACGCCCAGGTCCAGCGCAt           >  3:356025/1‑131 (MQ=255)
                                                                                                                                            acaAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTTCTCAACACGCCCAGGTCCAGCGCAt           <  4:356025/131‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                        
GCGGGGGCTGCACAGGCTGGCTTGACCTGTAGCGGCCTCTTCGCGGGCAAGCCCGCCCCCGCAGAGAAATCACAAGCCCAAGGGCTGTGGTATTCGCTGTGGAAGCGGGAGGCACGCGGCCCGGGTGTCAGCAGCAGGGCACAAGATTGCTGATATCCCCAAATGGCCAGGGCCGCGATGGCGTCTTAGTGCTTCTTCCCTTGCATGCTCGAAACCCATTGCTCGCTGCGACTATGCCCACTGGTTCTCAACACGCCCAGGTCCAGCGCATTTTCCCCAT  >  NC_002947/1652057‑1652336

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: