Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,705,690 C→A 6.4% E268* (GAA→TAA)  PP_4162 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,705,6900CA 6.4% 227.0 / 4.0 71E268* (GAA→TAA) PP_4162membrane protein
Reads supporting (aligned to +/- strand):  ref base C (31/35);  new base A (3/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GCAAGCCGCCGAACACGCCCAGCAGCACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAGTGCCAGCAGGGCCT  >  NC_002947/4705549‑4705836
                                                                                                                                             |                                                                                                                                                  
gcaagcCGCCGAACACGCCCAGCAGCACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTc                                                                                                                                             <  7:150145/149‑1 (MQ=255)
    gccgccGAACACGCCCAGCAGCACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCcaca                                                                                                                                         <  6:106347/149‑1 (MQ=255)
     ccgccgAACACGCCCAGCAGCACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCacac                                                                                                                                        >  4:344185/1‑149 (MQ=255)
       gccgAACACGCCCAGCAGCACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCa                                                                                                                                      <  7:84872/149‑1 (MQ=255)
          gAACACGCCCAGCAGCACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGcc                                                                                                                                   <  3:344185/149‑1 (MQ=255)
             cacGCCCAGCAGCACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAg                                                                                                                                <  8:157708/149‑1 (MQ=255)
              acGCCCAGCAGCACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGc                                                                                                                               >  3:3914/1‑149 (MQ=255)
                    agcagcACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGc                                                                                                                         <  4:55770/149‑1 (MQ=255)
                        gcaCGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCggg                                                                                                                     <  4:3914/149‑1 (MQ=255)
                         caCGCTCACCAGCGGCAGGTGGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGGTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCCGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGa                                                                                                                    <  2:199589/149‑1 (MQ=255)
                          aCGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGa                                                                                                                    <  1:138128/148‑1 (MQ=255)
                          aCGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGAt                                                                                                                   >  6:381308/1‑149 (MQ=255)
                          aCGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTAGCCCTTCCACAACAGCCAAGCGGTGGCCGGGAt                                                                                                                   >  1:249632/1‑149 (MQ=255)
                            gATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATc                                                                                                                  >  4:157649/1‑148 (MQ=255)
                                aCCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAgg                                                                                                              <  2:169952/148‑1 (MQ=255)
                                 ccAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGc                                                                                                             >  8:167666/1‑148 (MQ=255)
                                    gCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCTGGATCCCGGCCacc                                                                                                         <  8:85951/149‑3 (MQ=255)
                                     cGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTg                                                                                                        <  3:157649/149‑1 (MQ=255)
                                      ggCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTgg                                                                                                       <  6:186345/149‑1 (MQ=255)
                                       gCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTggc                                                                                                      >  8:115097/1‑149 (MQ=255)
                                        cAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTggcg                                                                                                     <  1:18293/149‑1 (MQ=255)
                                        cAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTggcg                                                                                                     <  7:167666/149‑1 (MQ=255)
                                            ttGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGccc                                                                                                 >  2:11568/1‑149 (MQ=255)
                                            ttGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGccc                                                                                                 >  8:47499/1‑149 (MQ=255)
                                                 gccACGGCTGAGCAGGTTCGGTGTGAGCACGTTGTCCATACCGCAGATGACGAACGTGCCCCACAAACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACCCCAGCCAAGCGGTGGCCGGGAGCCAGGCCAGTGGCGGCCCCATgg                                                                                            <  8:262835/149‑1 (MQ=255)
                                                  ccACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCCCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATggg                                                                                           <  5:333720/149‑1 (MQ=255)
                                                      ggCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATc                                                                                       <  1:31970/149‑1 (MQ=255)
                                                        cTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAg                                                                                     >  5:179618/1‑149 (MQ=255)
                                                           aTCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCt                                                                                  >  8:371385/1‑149 (MQ=255)
                                                              aGGTACGGTTTGATCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGa                                                                                <  5:82268/148‑1 (MQ=255)
                                                               ggTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGc                                                                              <  4:332777/149‑1 (MQ=255)
                                                                gTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCCTGCCGTATTCGCCCTTCCACACCAGCCACGCGGTGGCCGGGATCCAGGCCAGGGGCGGCCCCATGGGGATCAGGCTGAGCa                                                                             >  2:205613/1‑149 (MQ=255)
                                                                gTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCa                                                                             >  4:361654/1‑149 (MQ=255)
                                                                 tACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCAt                                                                            >  3:48950/1‑149 (MQ=255)
                                                                              cGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAAc                                                               <  6:179618/149‑1 (MQ=255)
                                                                                   tCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCa                                                           <  1:11568/148‑1 (MQ=255)
                                                                                     cacaCCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGc                                                         <  7:371385/148‑1 (MQ=255)
                                                                                      acacCGCTGTTGCCGAACGTGCCCCACACTCCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCAc                                                       <  8:307221/149‑1 (MQ=255)
                                                                                        acCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATgg                                                                                            >  4:371391/1‑110 (MQ=255)
                                                                                        acCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATgg                                                                                            <  3:371391/110‑1 (MQ=255)
                                                                                         ccGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGAt                                                    <  7:47499/149‑1 (MQ=255)
                                                                                                           ccccACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCgg                                   <  5:381308/148‑1 (MQ=255)
                                                                                                              cacacaCCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAAcccc                                                            <  1:45063/120‑1 (MQ=255)
                                                                                                              cacacaCCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAAcccc                                                            >  2:45063/1‑120 (MQ=255)
                                                                                                              cacacaCCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGAt                               >  5:381687/1‑149 (MQ=255)
                                                                                                               acacacCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCCCCCCGGCGAt                               >  2:227364/1‑148 (MQ=255)
                                                                                                               acacacCCAGGAACACCGCCATGCCGTATTAGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGAt                               >  6:57800/1‑148 (MQ=255)
                                                                                                                cacacCCAGGAACACCGCCATGCCGTATTAGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCa                             <  5:57800/149‑1 (MQ=255)
                                                                                                                   acCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAgg                           >  2:21427/1‑148 (MQ=255)
                                                                                                                     ccAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGAt                               <  5:300236/142‑1 (MQ=255)
                                                                                                                     ccAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGAt                               >  6:300236/1‑142 (MQ=255)
                                                                                                                     ccAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAg                        >  2:146193/1‑149 (MQ=255)
                                                                                                                      cAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCCCCCCGGCGATCAGGAAGc                       >  8:205847/1‑149 (MQ=255)
                                                                                                                       aGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGcc                      <  6:381687/149‑1 (MQ=255)
                                                                                                                       aGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGcc                      >  2:314421/1‑149 (MQ=255)
                                                                                                                        ggAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGcc                      >  5:294799/1‑148 (MQ=255)
                                                                                                                         gAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGa                    <  1:146193/149‑1 (MQ=255)
                                                                                                                          aaCACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGAt                   <  3:361654/149‑1 (MQ=255)
                                                                                                                           acacCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGAt                   >  7:34213/1‑148 (MQ=255)
                                                                                                                           acacCGCCATGCCGTATTAGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGAt                   >  4:60888/1‑148 (MQ=255)
                                                                                                                           acacCGCCATGCCGTATTAGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATc                  <  3:60888/149‑1 (MQ=255)
                                                                                                                             acCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAg                >  4:205969/1‑149 (MQ=255)
                                                                                                                               cGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAGTg              >  1:244206/1‑149 (MQ=255)
                                                                                                                                gCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGAGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAGTGc             <  7:183531/149‑1 (MQ=255)
                                                                                                                                  cATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAGTGCca           <  7:115097/149‑1 (MQ=255)
                                                                                                                                    tGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAGTGCcag          >  1:181777/1‑148 (MQ=255)
                                                                                                                                     gCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCg                                 >  1:168464/1‑124 (MQ=255)
                                                                                                                                     gCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCg                                 <  2:168464/124‑1 (MQ=255)
                                                                                                                                      ccGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAGTGCcagcag       >  6:253531/1‑149 (MQ=255)
                                                                                                                                           ttCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAGTGCCAGCAGGGCCt  >  4:384813/1‑149 (MQ=255)
                                                                                                                                           ttCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAGTGCCAGCAGGGCCt  >  6:125638/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GCAAGCCGCCGAACACGCCCAGCAGCACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACGAACGTGCCCCACACACCCAGGAACACCGCCATGCCGTATTCGCCCTTCCACACCAGCCAAGCGGTGGCCGGGATCCAGGCCAGTGGCGGCCCCATGGGGATCAGGCTGAGCATGAAGGTGACCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAGTGCCAGCAGGGCCT  >  NC_002947/4705549‑4705836

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