Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,494,572 T→G 5.4% K38T (AAG→ACG)  PP_3098 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,494,5720TG 5.4% 296.9 / 6.4 93K38T (AAG→ACG) PP_3098hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/49);  new base G (3/2);  total (42/51)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

GTCGAGCACGCTGCGGCTGATGAAGTCCGGCATGCCGTATTCAAGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTACACAACGGCCTTTCAG  >  NC_002947/3494430‑3494714
                                                                                                                                              |                                                                                                                                              
gTCGAGCACGCTGCGGCTGATGAAGTCCGGCATGCCGTATTCAAGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACt                                                                                                                                          >  4:58688/1‑149 (MQ=255)
    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                                                                                                >  3:168798/1‑149 (MQ=255)
                                          aaGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACt                                                                                                <  8:63388/149‑1 (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                                              >  7:326371/1‑149 (MQ=255)
                                                                                            gggACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCgg                                              >  3:391378/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:358389/1‑149 (MQ=255)
                                                                                                                              gcgGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTACACAAc            >  2:376107/1‑149 (MQ=255)
                                                                                                                               cgGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTACACAACg           <  7:162124/149‑1 (MQ=255)
                                                                                                                                 gACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGt                   <  2:139804/139‑1 (MQ=255)
                                                                                                                                 gACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGt                   >  1:139804/1‑139 (MQ=255)
                                                                                                                                    ctcgCTGGCCGAGTACAAGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCAtgt                              <  7:209777/122‑1 (MQ=255)
                                                                                                                                    ctcgCTGGCCGAGTACAAGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCAtgt                              >  8:209777/4‑125 (MQ=255)
                                                                                                                                     aTGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTACACAACGGCCttt     <  1:20801/149‑1 (MQ=255)
                                                                                                                                     aTGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTACACAACGGCCttt     <  5:384345/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                       >  2:389484/1‑123 (MQ=255)
                                                                                                                                             cTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATcc                       <  1:389484/123‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                              
GTCGAGCACGCTGCGGCTGATGAAGTCCGGCATGCCGTATTCAAGGATGGAACCACCGAGCTGGCTGAAACCGTCGAGCTCTTCGGCGACCAGGGACTGGCGGGTGTTCACCAGGACTTCGAGGTCGCGGACGATGCTGGCCTTGTACTCGGCCAGCGAGCAGGTGCGCTGCGCTGCGGGCTCGGCCGACTGGTAGGGACGATCGTCCAGCAGTCGGTCGAGCAGCGAGGGCAGCAGGCGGGCCTGGTTACTCATGTGTCATCCCTGTACACAACGGCCTTTCAG  >  NC_002947/3494430‑3494714

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: ‑