Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,494,579 C→A 5.3% E36* (GAG→TAG)  PP_3098 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,494,5790CA 5.3% 307.4 / 4.9 95E36* (GAG→TAG) PP_3098hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (39/51);  new base A (3/2);  total (42/53)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-01

AGCACGCTGCGGCTGATGAAGTCCGGCATGCCGTATTCAAGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTACACAACGGCCTTTCAGCCGCGCGC  >  NC_002947/3494434‑3494722
                                                                                                                                                 |                                                                                                                                               
aGCACGCTGCGGCTGATGAAGTCCGGCATGCCGTATTCAAGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGc                                                                                                                                              <  2:348977/149‑1 (MQ=255)
   aCGCTGCGGCTGATGAAGTCCGGCATGCCGTATTCAAGGTTGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCa                                                                                                                                            <  4:67965/148‑1 (MQ=255)
    cGCTGCGGCTGATGAAGTCCGGCATGCCGTATTCAAGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGATCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGc                                                                                                                                          <  6:190552/149‑1 (MQ=255)
    cGCTGCGGCGGCTGAAGTCCGGCATGCCGTATTCAAGGAGGGAACCACCGAGCTGGCTGAACCCGTCGAGCGCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGc                                                                                                                                          <  3:209321/149‑1 (MQ=255)
     gCTGCGGCTGATGAAGTCCGGCATGCCGTATTCAAGGATGGAACCACCGAGCTGGCTGACACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACAAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCCGc                                                                                                                                          >  5:108416/1‑148 (MQ=255)
       tGCGGCTGATGAAGTCCGGCATGCCGTATTCAAGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAg                                                                                                                                       >  3:337703/1‑149 (MQ=255)
             tgatgaAGTCCGGCATGCCGTATTCAAGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGtg                                                                                                                                 >  7:170921/1‑149 (MQ=255)
               atgaAGTCCGGCATGCCGTATTCAAGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGtgcg                                                                                                                               <  3:58688/149‑1 (MQ=255)
               atgaAGTCCGGCATGCCGTATTCAAGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCCTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGtgcg                                                                                                                               >  2:225824/1‑149 (MQ=255)
                tgaAGTCCGGCATGCCGTATTCAAGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGtgcgc                                                                                                                              >  5:342955/1‑149 (MQ=255)
                   aGTCCGGCATGCCGTATTCAAGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGtgcgctgc                                                                                                                           >  4:233584/1‑149 (MQ=255)
                         gCATGCCGTATTCAAGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGctgcgctgcg                                                                                                                     <  6:342955/149‑1 (MQ=255)
                         gCATGCCGTATTCAAGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGctgcgctgcg                                                                                                                     >  4:252306/1‑149 (MQ=255)
                            tGCCGTATTCAAGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGGTCACCAGGACTTCGAGGGCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGc                                                                                                                  <  5:377555/149‑1 (MQ=255)
                              ccGTATTCAAGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTc                                                                                                                <  1:274999/149‑1 (MQ=255)
                                 tATTCAAGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGCCTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGc                                                                                                             >  2:311713/1‑149 (MQ=255)
                                   ttCAAGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGTGACTGGCGGGTGTTCACCAGGACTTGGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCGAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGcc                                                                                                            <  6:127184/148‑1 (MQ=255)
                                      aaGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGAc                                                                                                         <  4:346669/148‑1 (MQ=255)
                                      aaGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGAc                                                                                                         >  3:346669/1‑148 (MQ=255)
                                      aaGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACt                                                                                                        <  8:63388/149‑1 (MQ=255)
                                      aaGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACt                                                                                                        >  3:168798/1‑149 (MQ=255)
                                         gATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGt                                                                                                     <  4:134086/149‑1 (MQ=255)
                                         gATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCGAGTACAAGGCCAGcctc                                                                                                                                       >  7:311294/1‑112 (MQ=255)
                                         gATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCGAGTACAAGGCCAGcatc                                                                                                                                       <  8:311294/115‑4 (MQ=255)
                                              aaCCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCGCCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGa                                                                                                <  4:61845/149‑1 (MQ=255)
                                              aaCCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGa                                                                                                >  5:262157/1‑149 (MQ=255)
                                              aaCCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGa                                                                                                >  1:23189/1‑149 (MQ=255)
                                              aaCCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCACGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGa                                                                                                >  8:280260/1‑149 (MQ=255)
                                                     gAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGt                                                                                         >  1:44309/1‑149 (MQ=255)
                                                      aGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTc                                                                                        <  3:233584/149‑1 (MQ=255)
                                                        ctggctgAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCca                                                                                      <  4:194584/149‑1 (MQ=255)
                                                           gctgAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCcagca                                                                                   <  4:277735/149‑1 (MQ=255)
                                                           gctgAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCcag                                                                                     >  4:275853/1‑147 (MQ=255)
                                                           gctgAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCcag                                                                                     <  1:375909/147‑1 (MQ=255)
                                                            ctgAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCcagcag                                                                                  >  6:25523/1‑149 (MQ=255)
                                                                     tCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGag                                                                         <  2:368494/149‑1 (MQ=255)
                                                                       gAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGagc                                                                        <  1:351716/148‑1 (MQ=255)
                                                                        aGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGagcag                                                                      >  8:162124/1‑149 (MQ=255)
                                                                          ctctTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGagcagc                                                                     >  4:213500/1‑148 (MQ=255)
                                                                          ctctTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCg                                                                    >  8:644/1‑149 (MQ=255)
                                                                             ttCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGa                                                                   <  7:239514/147‑1 (MQ=255)
                                                                             ttCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAgg                                                                 >  1:254732/1‑149 (MQ=255)
                                                                                ggCGTCCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGgc                                                               <  4:168798/148‑1 (MQ=255)
                                                                                ggCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGgca                                                              <  2:254732/149‑1 (MQ=255)
                                                                                 gCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGgcag                                                             <  1:356313/149‑1 (MQ=255)
                                                                                      cAGGGACTGGCGGGTGTTCACCAGGACTACGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGc                                                        >  1:224220/1‑149 (MQ=255)
                                                                                        gggACTGGCGGGTGTTCCCCCGGACTTCGAGGGCGCGGACGATGCTGGCCTTGTCCTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCGCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCgg                                                      <  6:108416/149‑1 (MQ=255)
                                                                                        gggACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCgg                                                      >  3:391378/1‑149 (MQ=255)
                                                                                        gggACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCgg                                                      >  7:326371/1‑149 (MQ=255)
                                                                                         ggACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCggg                                                     >  6:205749/1‑149 (MQ=255)
                                                                                            cTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCt                                                  <  4:337703/149‑1 (MQ=255)
                                                                                              ggCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCGGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTgg                                                <  7:146706/149‑1 (MQ=255)
                                                                                               gCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTgg                                                <  8:326371/148‑1 (MQ=255)
                                                                                               gCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGt                                               >  6:32603/1‑149 (MQ=255)
                                                                                                agggTGTTCACCAGTAATTCGAGGTCTCGGACGATGCTGGCCTTGTACTCTGCCAGCGAGCAGGTTCGCTGCGCTGCTGGATCGGCAGACTGGTTTGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGtt                                              <  4:67736/148‑1 (MQ=255)
                                                                                                    tgtTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTc                                          <  5:32603/149‑1 (MQ=255)
                                                                                                      ttCACCAGGACTTCGATGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCAt                                        <  1:225824/149‑1 (MQ=255)
                                                                                                      ttCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCAt                                        <  3:213500/149‑1 (MQ=255)
                                                                                                        cACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGTAGGCGGGCCTGGTTACTCAtgt                                      >  6:189548/1‑149 (MQ=255)
                                                                                                         aCCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCAtgtg                                     >  4:16035/1‑149 (MQ=255)
                                                                                                         aCCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCAtgtg                                     <  2:23189/149‑1 (MQ=255)
                                                                                                         aCCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCAtgt                                      <  8:170921/148‑1 (MQ=255)
                                                                                                          ccAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCAtgtg                                     >  3:364977/1‑148 (MQ=255)
                                                                                                           cAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATgtgt                                    >  6:225834/1‑148 (MQ=255)
                                                                                                             ggACTTCGAGGTCGCGTACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCAt                                 >  7:144009/1‑149 (MQ=255)
                                                                                                              gACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGt                                                                                                     >  2:302370/1‑80 (MQ=255)
                                                                                                              gACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGt                                                                                                     <  1:302370/80‑1 (MQ=255)
                                                                                                               aCTTCGAGGTCGCGGACGATGCTTGCCTTGTACTCTGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATcc                               <  2:224220/149‑1 (MQ=255)
                                                                                                               aCTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATcc                               <  2:44309/149‑1 (MQ=255)
                                                                                                                 ttCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTgg                                                >  8:229231/1‑130 (MQ=255)
                                                                                                                 ttCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTgg                                                <  7:229231/130‑1 (MQ=255)
                                                                                                                 ttCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCt                             >  3:195012/1‑149 (MQ=255)
                                                                                                                    gAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGt                           <  1:311713/148‑1 (MQ=255)
                                                                                                                    gAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTa                          <  5:205749/149‑1 (MQ=255)
                                                                                                                    gAGGTCACGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTa                          <  7:280260/149‑1 (MQ=255)
                                                                                                                          gcgGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTACACAAc                    >  2:376107/1‑149 (MQ=255)
                                                                                                                          gcgGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTACACAAc                    >  2:358389/1‑149 (MQ=255)
                                                                                                                           cgGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTACACAACg                   <  7:162124/149‑1 (MQ=255)
                                                                                                                             gACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGt                           >  1:139804/1‑139 (MQ=255)
                                                                                                                             gACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGt                           <  2:139804/139‑1 (MQ=255)
                                                                                                                                ctcgCTGGCCGAGTACAAGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCAtgt                                      <  7:209777/122‑1 (MQ=255)
                                                                                                                                ctcgCTGGCCGAGTACAAGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCAtgt                                      >  8:209777/4‑125 (MQ=255)
                                                                                                                                 aTGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTACACAACGGCCttt             <  5:384345/149‑1 (MQ=255)
                                                                                                                                 aTGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTACACAACGGCCttt             <  1:20801/149‑1 (MQ=255)
                                                                                                                                  tGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTACACAACGGCCTTTc            <  3:252306/149‑1 (MQ=255)
                                                                                                                                    cTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTACACAACGGCCTTTCa           <  7:644/148‑1 (MQ=255)
                                                                                                                                    cTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTACACAACGGCCTTTCAg          <  5:225834/149‑1 (MQ=255)
                                                                                                                                    cTGGCCGAGTACAAGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTACACAACGGCCTTTCAg          >  5:31148/1‑149 (MQ=255)
                                                                                                                                     tGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTac                         <  6:254040/133‑1 (MQ=255)
                                                                                                                                     tGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTac                         >  5:254040/1‑133 (MQ=255)
                                                                                                                                         cTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATcc                               <  1:389484/123‑1 (MQ=255)
                                                                                                                                         cTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATcc                               >  2:389484/1‑123 (MQ=255)
                                                                                                                                            ggaCTCGGCCAGCGAGCCGGCGCGCTCCGCTCCGGGCGCGGCCGACTGGCAGGGACGACCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCCTGTGGCATCCCTCTACACAACGGCCGTTCAGCcgcgcg   <  4:391378/146‑1 (MQ=255)
                                                                                                                                            gTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTACACAACGGCCTTTCAGCcgcgc    <  3:275853/147‑1 (MQ=255)
                                                                                                                                            gTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGGCATCCCTGTACACACCGGCCGTTCCGCcgcgcgc  >  6:230812/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
AGCACGCTGCGGCTGATGAAGTCCGGCATGCCGTATTCAAGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTACACAACGGCCTTTCAGCCGCGCGC  >  NC_002947/3494434‑3494722

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