Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,106,216 C→T 7.1% L252L (TTG→TTA)  PP_3613 ← L‑sorbosone dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,106,2160CT 7.1% 208.9 / 4.6 70L252L (TTG→TTA) PP_3613L‑sorbosone dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (32/33);  new base T (2/3);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

TGTAAGGCCAACCATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATG  >  NC_002947/4106074‑4106363
                                                                                                                                              |                                                                                                                                                   
tGTAAGGCCAACCATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCg                                                                                                                                               >  4:247855/1‑149 (MQ=255)
  tAAGGCCAACCATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCgg                                                                                                                                              <  5:350018/148‑1 (MQ=255)
  tAAGGCCAACCATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCg                                                                                                                                               <  3:164231/147‑1 (MQ=255)
    aGGCCAACCATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCACCGGCTTGCGTAACCCgacgg                                                                                                                                           <  2:197856/149‑5 (MQ=255)
       ccAACCATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCTATGCCACCGGCTTGCGTAACCCGACGGCAt                                                                                                                                        >  7:87410/1‑149 (MQ=255)
        cAACCATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCa                                                                                                                                         >  1:13497/1‑147 (MQ=255)
          aCCATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAgg                                                                                                                                     <  7:43791/149‑1 (MQ=255)
           ccATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCTATGCCACCGGCTTGCGTAACCCGACGGCAttgac                                                                                                                                    <  8:87410/149‑5 (MQ=38)
           ccATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAggg                                                                                                                                    >  4:309010/1‑149 (MQ=255)
           ccATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAggg                                                                                                                                    >  1:281761/1‑149 (MQ=255)
           acATAGAAGACGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGAGGCTTGATGGTCAATTCAGTCTGGTTACGCAAGCCGGTGGCATattg                                                                                                                                    >  1:286528/2‑146 (MQ=255)
            cATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCACCGGCTTGCGTAACCCGACGGCAttg                                                                                                                                      >  5:306912/1‑144 (MQ=255)
            cATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCACCGGCTTGCGTAACCCGACGGCAttg                                                                                                                                      <  6:306912/146‑3 (MQ=255)
              tAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGc                                                                                                                                                     >  4:302112/1‑129 (MQ=255)
              tAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGc                                                                                                                                                     >  4:296949/1‑129 (MQ=255)
              tAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCg                                                                                                                                 <  2:281761/149‑1 (MQ=255)
              tAGAAGCCGCCTTCACGTACCGAAGTGAGGGAGTCTGGAACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGAACACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGTTTACGCAAGCCGGTGGCATAGGGCCg                                                                                                                                 <  2:286528/149‑1 (MQ=255)
                gAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCgg                                                                                                                                <  6:367667/148‑1 (MQ=255)
                gAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGgt                                                                                                                               <  3:247855/149‑1 (MQ=255)
                 aaGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGgt                                                                                                                               >  5:57905/1‑148 (MQ=255)
                 aaGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGgt                                                                                                                               <  6:57905/148‑1 (MQ=255)
                   gccgccTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGtgtg                                                                                                                            >  2:371782/1‑149 (MQ=255)
                      gccTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTgcgc                                                                                                                         >  5:64680/1‑149 (MQ=255)
                           cACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGt                                                                                                                     <  1:41025/148‑1 (MQ=255)
                            aCGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGtct                                                                                                                   >  7:217876/1‑149 (MQ=255)
                                   gAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACACCTGCCCGGTCTGCGGCTGGATGGGCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATc                                                                                                            <  7:301288/149‑1 (MQ=255)
                                         aGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTg                                                                                                                                            <  8:202309/111‑1 (MQ=255)
                                          ggTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGc                                                                                                     >  7:361239/1‑149 (MQ=255)
                                           gTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTgc                                                                                                                           <  2:112184/126‑1 (MQ=255)
                                           gTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTgc                                                                                                                           >  1:112184/1‑126 (MQ=255)
                                              tgcTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTGGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATTCCGGCGGGTTACGCAAGCCGGTGGCATAGGGCCTGGGTGCGCCGGTCTCGGCATCAACCTGCCAg                                                                                                  <  8:361239/146‑1 (MQ=255)
                                              gTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGa                                                                                                 >  3:24146/1‑149 (MQ=255)
                                                   gCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATg                                                                                            >  6:368113/1‑149 (MQ=255)
                                                   gCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATg                                                                                            >  3:344886/1‑149 (MQ=255)
                                                     aCCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGc                                                                                          <  2:13497/149‑1 (MQ=255)
                                                          gTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCcgat                                                                                     <  7:163047/149‑1 (MQ=255)
                                                           tCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCcgatc                                                                                    <  3:309010/149‑1 (MQ=255)
                                                                  ccAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTcc                                                                             <  4:24146/149‑1 (MQ=255)
                                                                                tCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCcgcg                                                               >  8:8413/1‑149 (MQ=255)
                                                                                tCGTTGACAACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCAAGACCATGGCCCGATCGATCTCCACCTCCAGGCcgcg                                                               >  6:72511/1‑149 (MQ=255)
                                                                                 cGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGc                                                                                                     <  8:149801/110‑1 (MQ=255)
                                                                                 cGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGc                                                                                                     >  7:149801/1‑110 (MQ=255)
                                                                                  gTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCt                                                             >  5:287546/1‑149 (MQ=255)
                                                                                      accaccGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCgg                                                          <  5:368113/148‑1 (MQ=255)
                                                                                             tCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGtt                                                   >  4:60121/1‑148 (MQ=255)
                                                                                                acaACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATgg                                                                                           >  6:151578/1‑105 (MQ=255)
                                                                                                acaACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATgg                                                                                           <  5:151578/105‑1 (MQ=255)
                                                                                                 caACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGAtctc                                                                              <  2:264625/117‑1 (MQ=255)
                                                                                                 caACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGAtctc                                                                              >  1:264625/1‑117 (MQ=255)
                                                                                                 caACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGt                                              >  8:159160/1‑149 (MQ=255)
                                                                                                  aaCTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGtt                                             >  5:350629/1‑149 (MQ=255)
                                                                                                       cccGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGa                                          <  8:217876/147‑1 (MQ=255)
                                                                                                               tgagCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGGGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATAccc                                <  6:283300/146‑1 (MQ=255)
                                                                                                                   cTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCcaca                             >  7:269530/1‑148 (MQ=255)
                                                                                                                    tGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATa                           >  2:87824/1‑149 (MQ=255)
                                                                                                                      gATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATa                           <  8:269530/147‑1 (MQ=255)
                                                                                                                      gATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATAc                          <  4:344886/148‑1 (MQ=255)
                                                                                                                      gATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGAGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCcgatc                                                                                    <  4:563/90‑1 (MQ=255)
                                                                                                                      gATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGAGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCcgatc                                                                                    >  3:563/1‑90 (MQ=255)
                                                                                                                           tCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACaagaag                    >  3:93698/1‑149 (MQ=255)
                                                                                                                           tCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACaagaag                    >  4:332226/1‑149 (MQ=255)
                                                                                                                             aaTGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGc                   <  6:350629/148‑1 (MQ=255)
                                                                                                                              aTGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCg                                                           <  1:334592/107‑1 (MQ=255)
                                                                                                                              aTGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCg                                                           >  2:334592/1‑107 (MQ=255)
                                                                                                                              aTGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGa                 <  7:159160/149‑1 (MQ=255)
                                                                                                                              aTGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGa                 <  1:371782/149‑1 (MQ=255)
                                                                                                                              aTGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGa                 <  7:8413/149‑1 (MQ=255)
                                                                                                                              aTGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGa                 <  5:72511/149‑1 (MQ=255)
                                                                                                                                  cGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGATCCGATACCCACATACAAGAAGCGACCAt             >  4:249954/1‑149 (MQ=255)
                                                                                                                                           aCGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCTATACCcac                              <  6:160661/123‑1 (MQ=255)
                                                                                                                                           aCGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCcac                              >  5:160661/1‑123 (MQ=255)
                                                                                                                                             gCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTCTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATg  >  5:354287/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
TGTAAGGCCAACCATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATG  >  NC_002947/4106074‑4106363

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