Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,106,219 G→C 7.1% G251G (GGC→GGG)  PP_3613 ← L‑sorbosone dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,106,2190GC 7.1% 215.6 / 6.9 70G251G (GGC→GGG) PP_3613L‑sorbosone dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (32/33);  new base C (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.37e-01

TGTAAGGCCAACCATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCA  >  NC_002947/4106074‑4106367
                                                                                                                                                 |                                                                                                                                                    
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)
              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                                                                                                                               >  1:112184/1‑126 (MQ=255)
                                           gTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTgc                                                                                                                               <  2:112184/126‑1 (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)
                                                                                                                                                 gCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCa  <  6:64680/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TGTAAGGCCAACCATAGAAGCCGCCTTCACGTACCGAAGTGAGGTAGTCTGGCACCAGGTCGGGGCCCAGCTCATCCCGCTCGTTGACCACCGTCCACAACTGCCCGGTCTGCGGCTGGATGGTCAATGCCGTCGGGTTACGCAAGCCGGTGGCATAGGGCCGGTGTGCGCCGGTCTCGGCATCAACCTGCCAGACCATGGCCCGATCGATCTCCACCTCCAGGCCGCGCTCGGTGATGTTGCTGTTCGAGCCGATACCCACATACAAGAAGCGACCATCCGGGCTGATGGTCA  >  NC_002947/4106074‑4106367

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