Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,246,194 C→T 6.3% F183F (TTC→TTT)  PP_1087 → OmpA family outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,246,1940CT 6.3% 190.7 / 2.1 63F183F (TTC→TTT) PP_1087OmpA family outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (32/27);  new base T (2/2);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

GCGCAACTACACCGGTGAGGCGGGCAGGCCGATGGAGGTGCGTGTGCTGCCGGTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCATTCCGACAACAGTGGCAACCGCCTTACCAACCG  >  NC_002947/1246047‑1246334
                                                                                                                                                   |                                                                                                                                            
gcgcAACTACACCGGTGAGGCGGGCAGGCCGATGGAGGTGCGTGTCCTGCCGGTCAGCTTCGCCAAAGCCTACAAAGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGCTGCGCCAGACACAGGTGGGTTTcc                                                                                                                                             >  7:369683/1‑149 (MQ=255)
  gcAACTACACCGGTGAGGCGGGCAGGCCGATGGAGGTGCGTGTGCTGCCGGTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTcccc                                                                                                                                           >  2:18500/1‑149 (MQ=255)
        acacCGGTGAGGCGGGCAGGCCGATGGAGGTGCGTGTGCTGCCGGTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGc                                                                                                                                     >  2:218008/1‑149 (MQ=255)
         cacCGGTGAGGCGGGCAGGCCGATGGAGGTGCGTGTGCTGCCGGTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTTCCCGGGGGCg                                                                                                                                    >  8:334410/1‑149 (MQ=255)
         cacCGGTGAGGCGGGCAGGCCGATGGAGGTGCGTGTGCTGCCGGTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCg                                                                                                                                    <  3:257998/149‑1 (MQ=255)
          acCGGTGAGGCGGGCAGGCCGATGGAGGTGCGTGTGCTGCCGGTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCg                                                                                                                                    >  4:257998/1‑148 (MQ=255)
              gTGAGGCGGGCAGGCCGATGGAGGTGCGTTTGCTGCCGGTCATCTTCGGCATAGCCTACAATTAATATCAGGTCTGTGCTGGCATGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATc                                                                                                                               <  8:187096/149‑1 (MQ=255)
              gTGAGGCGGGCAGGCCGATGGAGGTGCGTGTGCTGCCGGTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATc                                                                                                                               <  4:287782/149‑1 (MQ=255)
              gTGAGGCGGGCAGGCCGATGGAGGTGCGTGTGCTGCCGGTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATc                                                                                                                               >  3:287782/1‑149 (MQ=255)
                 aGGCGGGCAGGCCGATGGAGGTGCGTGTGCTGCCGGTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGa                                                                                                                             <  1:305306/148‑1 (MQ=255)
                   gCGGGCAGGCCGATGGAGGTGCGTGTGCTGCAGGTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTTCGCCAGACACAGGTGGGTTTCCCCGTGGGCGGTATCGAGCt                                                                                                                          >  5:350834/1‑149 (MQ=255)
                   gCGGGCAGGCCGAGGGAGGTGCGGGTGCTGCCGGTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGAGGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCt                                                                                                                          <  7:326069/149‑1 (MQ=255)
                        cAGGCCGATGGAGGTGCGTGTGCTGCCGGTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACg                                                                                                                     <  5:186712/149‑1 (MQ=255)
                                tGGAGGTGCGTGTGCTGCCGGTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCggg                                                                                                                                        >  7:753/1‑122 (MQ=255)
                                tGGAGGTGCGTGTGCTGCCGGTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCggg                                                                                                                                        <  8:753/122‑1 (MQ=255)
                                             gctgcCGGTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGcagcag                                                                                                               >  3:201231/1‑134 (MQ=255)
                                             gctgcCGGTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGcagcag                                                                                                               <  4:201231/134‑1 (MQ=255)
                                               tgcCGGTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGAc                                                                                              >  5:87101/1‑149 (MQ=255)
                                                   ggTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCgg                                                                                                    <  5:267823/139‑1 (MQ=255)
                                                   ggTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCgg                                                                                                    >  6:267823/1‑139 (MQ=255)
                                                       aGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTTCCCGGGGGc                                                                                                                                     >  3:273591/1‑102 (MQ=255)
                                                       aGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTTCCCGGGGGc                                                                                                                                     <  4:273591/102‑1 (MQ=255)
                                                       aGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGa                                                                                                                             >  5:27451/1‑110 (MQ=255)
                                                       aGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGa                                                                                                                             <  6:27451/110‑1 (MQ=255)
                                                               aaaaGCCTACAATGACTATCAGGGCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACAGAGGTGGGTTTCCCCGGGGGCGGTAGCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACa                                                                              <  5:120078/148‑1 (MQ=255)
                                                                    ccTACAATGACTATGAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTTCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGaaa                                                                         <  7:334410/149‑1 (MQ=255)
                                                                    ccTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACAt                                                                             >  1:179325/1‑145 (MQ=255)
                                                                    ccTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACAt                                                                             <  2:179325/145‑1 (MQ=255)
                                                                    ccTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGaaa                                                                         <  5:121365/149‑1 (MQ=255)
                                                                     cTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCC‑GGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCg                                                                                    >  5:119585/1‑136 (MQ=255)
                                                                     cTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCC‑GGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCg                                                                                    <  6:119585/136‑1 (MQ=255)
                                                                     cTACAATGACTATAAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAg                                                                        >  7:63727/1‑149 (MQ=255)
                                                                         aaTGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGAt                                                                                  >  1:103387/1‑135 (MQ=255)
                                                                         aaTGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGAt                                                                                  <  2:103387/135‑1 (MQ=255)
                                                                                   aGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTg                                                          <  6:87101/149‑1 (MQ=255)
                                                                                        tgtgCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCa                                                     >  1:25167/1‑149 (MQ=255)
                                                                                          tgCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATg                                                   >  3:306312/1‑149 (MQ=255)
                                                                                          tgCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATg                                                   >  7:293008/1‑149 (MQ=255)
                                                                                           gCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGt                                                  >  7:72910/1‑149 (MQ=255)
                                                                                            cTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTc                                                 <  4:306312/149‑1 (MQ=255)
                                                                                                 aaGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCt                                            <  8:293008/149‑1 (MQ=255)
                                                                                                 aaGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCt                                            >  5:53666/1‑149 (MQ=255)
                                                                                                      gTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCt                                            >  6:13172/1‑144 (MQ=255)
                                                                                                      gTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCt                                            <  5:13172/144‑1 (MQ=255)
                                                                                                       ttGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATg                                       <  1:218008/148‑1 (MQ=255)
                                                                                                         gCCGTTGAACTACGACCAGGTTCGCGAGACTGAGGTGGGTTTCCCGGGTGGCGGTATCTAGGTGGACGCAGCATCCAAGGCGCGGTTGGACGTGATCCTTGATTACATTAAAGCCGTCCCTACGGTGAACCATGTCGAGCTCAATggg                                     <  8:369683/148‑1 (MQ=255)
                                                                                                           cGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCAt                                  >  1:45653/1‑149 (MQ=255)
                                                                                                           cGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCAt                                  >  2:306896/1‑149 (MQ=255)
                                                                                                            gATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCAtt                                 >  6:194485/1‑149 (MQ=255)
                                                                                                                  cTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCATTCCGac                            >  3:145410/1‑148 (MQ=255)
                                                                                                                  cTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCATTCCGac                            >  5:135008/1‑148 (MQ=255)
                                                                                                                        ccAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCATTCCGACAACAGTg                     >  3:372295/1‑149 (MQ=255)
                                                                                                                        ccAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCATTCCGACAACAGTg                     <  8:63727/149‑1 (MQ=255)
                                                                                                                        ccAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCATTCCGACAACAGTg                     >  2:306285/1‑149 (MQ=255)
                                                                                                                           ggTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCATTCCGa                             >  4:342727/1‑138 (MQ=255)
                                                                                                                           ggTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCATTCCGa                             <  3:342727/138‑1 (MQ=255)
                                                                                                                             tGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCATTCCGACAACAGTGGCAAc                >  2:295010/1‑149 (MQ=255)
                                                                                                                                gcCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCATTCCGACAACAGTGGCAACCGc             >  8:269319/1‑149 (MQ=255)
                                                                                                                                  cAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCATTCCGACAACAGTGGCAACCGCCt           <  7:269319/149‑1 (MQ=255)
                                                                                                                                  cAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCATTCCGACAACAGTGGCAACCGCCt           >  2:247230/1‑149 (MQ=255)
                                                                                                                                     acacAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCATTCCGACAACAGTGGCAACCGCCTTAc        <  1:18500/149‑1 (MQ=255)
                                                                                                                                     acacAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCATTCCGACAACAGTGGCAACCGCCTTAc        <  2:25167/149‑1 (MQ=255)
                                                                                                                                      cacaGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCATTCCGACAACAGTGGCAACCGCCTTAc        >  5:384836/1‑148 (MQ=255)
                                                                                                                                         aGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAATCATGTCGAGCTCAATGGGCATTCCGACAACAGTGGCAACCGCCTTACCAAc    <  1:306896/149‑1 (MQ=255)
                                                                                                                                           gTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCATTCCGACAACAGTGGCAACCGCCTTACCAACcg  >  4:184918/1‑149 (MQ=255)
                                                                                                                                                   ccccGGGGGCAGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCATTCCGACAACAGTGGCAACCGCCTTAc        >  5:186877/1‑135 (MQ=255)
                                                                                                                                                   ccccGGGGGCAGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCATTCCGACAACAGTGGCAACCGCCTTAc        <  6:186877/135‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
GCGCAACTACACCGGTGAGGCGGGCAGGCCGATGGAGGTGCGTGTGCTGCCGGTCAGCTTCGCCAAAGCCTACAATGACTATCAGGTCTGTGCTGGCAAGCTGTTGCCGATGAACTACGACCAGGTGCGCCAGACACAGGTGGGTTTCCCCGGGGGCGGTATCGAGCTGGACGCAGCAGCCAAGGCGCGGCTGGACGTGATCCTCGATTACATGAAAGCCGACCCGACGGTGAACCATGTCGAGCTCAATGGGCATTCCGACAACAGTGGCAACCGCCTTACCAACCG  >  NC_002947/1246047‑1246334

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