Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,106,210 G→C 6.8% N254K (AAC→AAG PP_3613 ← L‑sorbosone dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,106,2100GC 6.8% 233.0 / 6.6 73N254K (AAC→AAGPP_3613L‑sorbosone dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (35/32);  new base C (2/3);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

CCAGTAGCTGTAAGGCCAACCATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCAT  >  NC_002947/4106066‑4106352
                                                                                                                                                |                                                                                                                                              
ccAGTAGCTGTAAGGCCAACCATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTAc                                                                                                                                            >  6:350018/1‑149 (MQ=255)
ccAGTAGCTGTAAGGCCAACCATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTAc                                                                                                                                            <  2:320492/149‑1 (MQ=255)
  aGTAGCTGTAAGGCCAACCATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACg                                                                                                                                           >  8:301288/1‑148 (MQ=255)
  aGTAGCTGTAAGGCCAACCATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACg                                                                                                                                           >  3:181400/1‑148 (MQ=255)
        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                                                                                                                         >  1:281761/1‑149 (MQ=255)
                   ccATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAggg                                                                                                                         >  4:309010/1‑149 (MQ=255)
                   acATAGAAGACGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGAGGCTTGATGGTCAATTCAGTCTGGTTACGCAAGCCGGTGGCATattg                                                                                                                         >  1:286528/2‑146 (MQ=255)
                    cATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCACCGGCTTGCGTAACCCGACGGCAttg                                                                                                                           <  6:306912/146‑3 (MQ=255)
                    cATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCACCGGCTTGCGTAACCCGACGGCAttg                                                                                                                           >  5:306912/1‑144 (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                                                                                 >  3:344886/1‑149 (MQ=255)
                                                           gCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATg                                                                                 >  6:368113/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                                                                         >  3:563/1‑90 (MQ=255)
                                                                                                                              gATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGAGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCcgatc                                                                         <  4:563/90‑1 (MQ=255)
                                                                                                                                   tCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACaagaag         >  4:332226/1‑149 (MQ=255)
                                                                                                                                   tCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACaagaag         >  3:93698/1‑149 (MQ=255)
                                                                                                                                     aaTGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGc        <  6:350629/148‑1 (MQ=255)
                                                                                                                                      aTGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCg                                                >  2:334592/1‑107 (MQ=255)
                                                                                                                                      aTGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCg                                                <  1:334592/107‑1 (MQ=255)
                                                                                                                                      aTGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGa      <  7:8413/149‑1 (MQ=255)
                                                                                                                                      aTGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGa      <  7:159160/149‑1 (MQ=255)
                                                                                                                                      aTGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGa      <  1:371782/149‑1 (MQ=255)
                                                                                                                                      aTGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGa      <  5:72511/149‑1 (MQ=255)
                                                                                                                                          cGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGATCCGATACCCACATACAAGAAGCGACCAt  >  4:249954/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                              
CCAGTAGCTGTAAGGCCAACCATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCAT  >  NC_002947/4106066‑4106352

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: