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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,106,3570GC 8.9% 237.5 / 10.4 78S205R (AGC→AGG) ‡PP_3613L‑sorbosone dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (39/32);  new base C (3/4);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

GTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGCTACATACAGCTTGCCGTTG  >  NC_002947/4106210‑4106504
                                                                                                                                                   |                                                                                                                                                   
gTTACGCAAGCCGGTGGCATAGGGCCGGAGTTCGCCGGACTCGGCATAAACCTGCCAGACCATGGCCCGATCGAGCTACACCTCCCGGCCGCCCTCGGTGATGATGCTCTTAGCGCCGATAACCACATACAAGATGCGAACATACGGGc                                                                                                                                                    >  6:286140/1‑149 (MQ=255)
     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                                                                                                                                     <  2:69766/119‑1 (MQ=255)
                                             aTCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCt                                                                                                                                     >  1:69766/1‑119 (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                                                                  >  3:101129/1‑149 (MQ=255)
                                                                                  tCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGt                                                                  >  2:475/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                          >  8:166367/1‑149 (MQ=255)
                                                                                                                          ccACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTgg                          >  1:243050/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                                                   >  6:369345/1‑119 (MQ=255)
                                                                                                                               tACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTa                                                   <  5:369345/119‑1 (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:43827/149‑1 (MQ=255)
                                                                                                                                                  ggCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGCTACATACAGCTTGCCGTTg  <  4:101129/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCAGGGCCTTGGTCCAATGGTGGTTGATTTCAGCCGGCAGGTCGGTGACCTTGGTTGGCGGGCCGTCGGCCTGGGTCTGGCCTTCCTTGTAGTCAAAGCGCACCAGCGCATCCTGGTTGGCTACATACAGCTTGCCGTTG  >  NC_002947/4106210‑4106504

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: