Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,106,359 T→G 8.9% S205R (AGC→CGC) ‡ PP_3613 ← L‑sorbosone dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,106,3590TG 8.9% 233.8 / 13.1 79S205R (AGC→CGC) ‡PP_3613L‑sorbosone dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (38/34);  new base G (3/4);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.74e-01

GCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGCTACATACAGCTTGCCGTTGGCA  >  NC_002947/4106215‑4106507
                                                                                                                                                |                                                                                                                                                    
gCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGCTGATGTTGCTGTTCGAGCCGATACCCACATACAAGACGCGACCATCCGGGCTGATg                                                                                                                                                  >  7:153771/1‑149 (MQ=255)
gCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGAGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATg                                                                                                                                                  >  2:83879/1‑149 (MQ=255)
  aaGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGt                                                                                                                                                >  8:31197/1‑149 (MQ=255)
   aGCCTGTGGCATAGGGCCTGTGTGCGCCGTTCTCGGCATCAACCTGGCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGTTGTTGCTGTTCGAGCGGATGCCCACATACAAGAAGCGACCATCAGCCCGGATGGTc                                                                                                                                               <  7:159549/149‑1 (MQ=255)
   aGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCTATACCCACATACAAGAAGCGACCATCAGCCCGGATGGTc                                                                                                                                               <  8:206866/149‑1 (MQ=255)
    gCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGt                                                                                                                                                >  3:43827/1‑147 (MQ=255)
    gCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCa                                                                                                                                              <  3:340927/149‑1 (MQ=255)
     ccGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAg                                                                                                                                             <  5:392204/149‑1 (MQ=255)
     ccGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCAGCCCGGATGGTcgc                                                                                                                                             <  4:280712/149‑3 (MQ=255)
      cGGTGTCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAgg                                                                                                                                            <  6:252954/149‑1 (MQ=255)
      cGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAgg                                                                                                                                            <  5:351087/149‑1 (MQ=255)
      cGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCAGCCCGGATGGTcgct                                                                                                                                            <  1:99687/149‑4 (MQ=255)
         tGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGcc                                                                                                                                         <  5:122019/149‑1 (MQ=255)
            cATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTg                                                                                                                                      <  7:200317/149‑1 (MQ=255)
             aTAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTgg                                                                                                                                     >  6:226690/1‑149 (MQ=255)
                 ggCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTcc                                                                                                                                  >  2:286309/1‑148 (MQ=255)
                  gCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCaa                                                                                                                                <  5:226690/149‑1 (MQ=255)
                     ggTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAAtgg                                                                                                                             >  4:357005/1‑149 (MQ=255)
                                 tctcGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCa                                                                                                                 <  2:234291/149‑1 (MQ=255)
                                    cGGCATCACCCTGCCAGACCATGGCCCGATCGAGCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGcc                                                                                                              <  8:318417/149‑1 (MQ=255)
                                      gCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCgg                                                                                                            <  6:238183/149‑1 (MQ=255)
                                        aTCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCt                                                                                                                                        >  1:69766/1‑119 (MQ=255)
                                        aTCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCt                                                                                                                                        <  2:69766/119‑1 (MQ=255)
                                           aaCCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGt                                                                                                       >  4:310725/1‑149 (MQ=255)
                                             ccTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCg                                                                                                     >  8:72774/1‑149 (MQ=255)
                                              cTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCgg                                                                                                    >  2:72710/1‑149 (MQ=255)
                                                 ccAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGa                                                                                                 <  3:357005/149‑1 (MQ=255)
                                                       cATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCttgg                                                                                           >  1:279439/1‑149 (MQ=255)
                                                         gggCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGATCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGGGGTTGATTTCAGCCGGCAGGTCGGTGACCttggtt                                                                                         <  2:287571/148‑1 (MQ=255)
                                                           gCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCttggttgg                                                                                       >  2:63838/1‑149 (MQ=255)
                                                              cgatcgatCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCgg                                                                                    >  5:86575/1‑149 (MQ=255)
                                                               gatcgatcTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCggg                                                                                   >  7:17244/1‑149 (MQ=255)
                                                                   gatcTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGt                                                                               >  6:279435/1‑149 (MQ=255)
                                                                     tctcCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCg                                                                             >  5:11823/1‑149 (MQ=255)
                                                                         cACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTAGGTGACCTTGGTTGGCGGGCCGTCGGCCt                                                                         >  3:16347/1‑149 (MQ=255)
                                                                          aCCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTg                                                                        >  4:237985/1‑149 (MQ=255)
                                                                           ccTCCAGGCCGAGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTgg                                                                       <  1:83879/149‑1 (MQ=255)
                                                                             tCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGt                                                                     >  2:475/1‑149 (MQ=255)
                                                                             tCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGt                                                                     >  3:101129/1‑149 (MQ=255)
                                                                                     gcgcTCGGTGATGTTGATGTTCGATCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCATTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGAGGGCCGTCgg                                                                            >  3:139484/1‑134 (MQ=255)
                                                                                         tCGGCGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGccttcctt                                                         <  5:201220/149‑1 (MQ=255)
                                                                                           ggTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGccttcctt                                                         <  7:72774/147‑1 (MQ=255)
                                                                                           ggTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTg                                                        <  5:279435/148‑1 (MQ=255)
                                                                                            gTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTa                                                      <  1:475/149‑1 (MQ=255)
                                                                                               aTGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGt                                                    <  1:72710/148‑1 (MQ=255)
                                                                                               aTGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGt                                                    >  2:148585/1‑148 (MQ=255)
                                                                                                  ttGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCaaa                                                <  8:17244/149‑1 (MQ=255)
                                                                                                        ttCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGTGTCTGGCCTTCCTTGTAGTCAAAGCGCAc                                          <  1:148585/149‑1 (MQ=255)
                                                                                                           gAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAg                                       >  6:345841/1‑149 (MQ=255)
                                                                                                           gAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCCGGCCTTCCTTGTAGTCAAAGCGCACCa                                        <  2:279439/148‑1 (MQ=255)
                                                                                                                 aTACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATc                                 >  7:238807/1‑149 (MQ=255)
                                                                                                                  tACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATcc                                <  5:345841/149‑1 (MQ=255)
                                                                                                                    cccACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTg                              <  7:31197/149‑1 (MQ=255)
                                                                                                                     ccACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCGTTGTAGTCAAAGCGCACCAGCGCATCCTgg                             >  4:248994/1‑149 (MQ=255)
                                                                                                                     ccACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTgg                             >  1:243050/1‑149 (MQ=255)
                                                                                                                     ccACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTgg                             >  8:166367/1‑149 (MQ=255)
                                                                                                                     ccACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTgg                             >  3:104843/1‑149 (MQ=255)
                                                                                                                      cacaTACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGt                            <  3:310725/149‑1 (MQ=255)
                                                                                                                       acaTACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCTTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGtt                           <  1:63838/149‑1 (MQ=255)
                                                                                                                       acaTACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATCGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTCGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCTAAGCGAACCAGCGCATCCTGGtt                           >  8:238717/1‑149 (MQ=255)
                                                                                                                          tACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTa                                                      <  5:369345/119‑1 (MQ=255)
                                                                                                                          tACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTa                                                      >  6:369345/1‑119 (MQ=255)
                                                                                                                            cAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGCta                      <  3:237985/149‑1 (MQ=255)
                                                                                                                             aagaagCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGtt                           <  4:175278/143‑1 (MQ=255)
                                                                                                                             aagaagCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGtt                           >  3:175278/1‑143 (MQ=255)
                                                                                                                              agaagCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGc                        <  4:49016/145‑1 (MQ=255)
                                                                                                                              agaagCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGc                        >  3:49016/1‑145 (MQ=255)
                                                                                                                                aagCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGCtacata                  >  7:47459/1‑149 (MQ=255)
                                                                                                                                  gCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGCtacata                  >  3:73325/1‑147 (MQ=255)
                                                                                                                                  ctgACCATCAGCCCGGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGCtacata                  >  7:212620/3‑147 (MQ=255)
                                                                                                                                  ctgACCATCAGCCCGGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGCtacata                  >  3:205417/3‑147 (MQ=255)
                                                                                                                                  ctgACCATCAGCCCGGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTGGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGCtacata                  >  4:72465/3‑147 (MQ=255)
                                                                                                                                    gACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGCta                      <  8:223717/141‑1 (MQ=255)
                                                                                                                                    gACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGCta                      >  7:223717/1‑141 (MQ=255)
                                                                                                                                     aCCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTgg                         <  7:251934/137‑1 (MQ=255)
                                                                                                                                     aCCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTgg                         >  8:251934/1‑137 (MQ=255)
                                                                                                                                       cATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGCTACATACAGCTTg           >  7:265158/1‑149 (MQ=255)
                                                                                                                                             ggCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGCTACATACAGCTTGCCGTTg     <  4:101129/149‑1 (MQ=255)
                                                                                                                                             ggCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGCTACATACAGCTTGCCGTTg     <  4:43827/149‑1 (MQ=255)
                                                                                                                                                tGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGCTACATACAGCTTGCCGTTGGCa  >  7:362631/1‑149 (MQ=255)
                                                                                                                                                tGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGCTACATACAGCTTGCCGTTGGCa  >  3:52613/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGCTACATACAGCTTGCCGTTGGCA  >  NC_002947/4106215‑4106507

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: