Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,011,120 A→T 15.0% intergenic (+39/+43) norM → / ← PP_5263 multidrug resistance protein NorM/GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,011,1200AT15.0% 219.2 / 26.9 90intergenic (+39/+43)norM/PP_5263multidrug resistance protein NorM/GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (32/41);  new base T (3/10);  total (38/52)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTGGGGCTGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCAG  >  NC_002947/6010971‑6011268
                                                                                                                                                       |                                                                                                                                                     
gcgcTGGGGCTGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGa                                                                                                                                                       <  1:24124/150‑1 (MQ=255)
gcgcTGGGGCTGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGa                                                                                                                                                       <  1:24107/150‑1 (MQ=255)
 cgcTGGGGCTGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGAt                                                                                                                                                      <  2:709731/150‑1 (MQ=255)
  gcTGGGGCTGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGAtt                                                                                                                                                     >  1:2713185/1‑150 (MQ=255)
     ggggCTGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTgcg                                                                                                                                                 <  2:1551452/150‑1 (MQ=255)
      gggCTGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTgcgc                                                                                                                                                <  2:591633/150‑1 (MQ=255)
          tGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTc                                                                                                                                            <  1:2870097/150‑1 (MQ=255)
          tGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTc                                                                                                                                            <  2:2090640/150‑1 (MQ=255)
           ggCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTcc                                                                                                                                           <  1:18496/150‑1 (MQ=255)
            gCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGAGTGCGCGGTccc                                                                                                                                          >  2:3959556/1‑150 (MQ=255)
               tgtgCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCtgt                                                                                                                                       <  1:3937275/150‑1 (MQ=255)
               tgtgCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCtgt                                                                                                                                       <  2:2713185/150‑1 (MQ=255)
               tgtgCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCtgt                                                                                                                                       <  1:1317588/150‑1 (MQ=255)
                gtgCGGCTGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCtgtg                                                                                                                                      <  1:3729020/150‑1 (MQ=255)
                gtgCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCtgtg                                                                                                                                      <  1:3355709/150‑1 (MQ=255)
                 tgatcgggTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTgg                                                                                                                                     <  1:1954315/144‑1 (MQ=255)
                 tgCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTgg                                                                                                                                     <  1:524266/150‑1 (MQ=255)
                 tgCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTgg                                                                                                                                     <  1:2556429/150‑1 (MQ=255)
                 tgCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTgg                                                                                                                                     <  1:2680071/150‑1 (MQ=255)
                 tgCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTgg                                                                                                                                     <  2:4329804/150‑1 (MQ=255)
                   cggcggTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGa                                                                                                                                   <  1:2024784/150‑1 (MQ=255)
                    ggcggTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAg                                                                                                                                  <  1:4123082/150‑1 (MQ=255)
                     gcggTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGc                                                                                                                                 <  1:3343825/150‑1 (MQ=255)
                      caggataacgacagggtggtggggatactcggggataccgagctttaccacgcgctgctcggcaagaaccacggttgaagcaattggggccgcttcgcggcccttcgcgggcacgcccgctcccacagGTTTGCTGAGATCTCTGTGGGa                                                                                                                                   <  1:988953‑M2/22‑1 (MQ=255)
                      caggataacgacagggtggtggggatactcggggataccgagctttaccacgcgctgctcggcaagaaccacggttgaagcaattggggccgcttcgcggcccttcgcgggcacgcccgctcccacagGTTTGCTGAGATCTCTGTGGGa                                                                                                                                   <  1:792974‑M2/22‑1 (MQ=255)
                      caggataacgacagggtggtggggatactcggggataccgagctttaccacgcgctgctcggcaagaaccacggttgaagcaattggggccgcttcgcggcccttcgcgggcacgcccgctcccacagGTTTGCTGAGATCTCTGTGGGa                                                                                                                                   <  1:4443332‑M2/22‑1 (MQ=255)
                       ggTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCCTGTGGGAGCgg                                                                                                                               <  2:2464843/150‑1 (MQ=37)
                       ggTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCCTGTGGGAGCgg                                                                                                                               <  2:2764103/150‑1 (MQ=37)
                       ggTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCCTGTGGGAGCgg                                                                                                                               <  2:2936308/150‑1 (MQ=37)
                       ggTGGCGCTGACCTTTGGCTGTGAGTGGCGGATGAAACGGTTGTAGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCAATTCGCGGGCAACCCCGCCCCCACCGGCATGCCCGGTACCGGGTGGAGCgg                                                                                                                               >  1:825210/1‑150 (MQ=255)
                        ggataacgacagggtggtggggatactcggggataccgagctttaccacgcgctgctcggcaagaaccacggttgaagcaattggggccgcttcgcggcccttcgcgggcacgcccgctcccacagGTTTGCTGAGATCTCTGTGGGAGc                                                                                                                                 <  2:3856533‑M2/24‑1 (MQ=255)
                         tGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGc                                                                                                                             <  1:5401921/150‑1 (MQ=255)
                         tGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGc                                                                                                                             <  1:3861861/150‑1 (MQ=255)
                         gataacgacagggtggtggggatactcggggataccgagctttaccacgcgctgctcggcaagaaccacggttgaagcaattggggccgcttcgcggcccttcgcgggcacgcccgctcccacagGTTTGCTGAGATCTCTGTGGGAGCg                                                                                                                                <  1:5417779‑M2/25‑1 (MQ=255)
                          ataacgacagggtggtggggatactcggggataccgagctttaccacgcgctgctcggcaagaaccacggttgaagcaattggggccgcttcgcggcccttcgcgggcacgcccgctcccacagGTTTGCTGAGATCTCTGTGGGAGCgg                                                                                                                               <  1:1234430‑M2/26‑1 (MQ=255)
                           gcgcTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCCTGTGGGAGCGGGCGt                                                                                                                           <  1:1480252/150‑1 (MQ=25)
                           gcgcTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGt                                                                                                                           <  1:309501/150‑1 (MQ=255)
                            cgcTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCCTGTGGGAGAGGGCGTg                                                                                                                          <  1:4336560/150‑1 (MQ=25)
                             gcTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGc                                                                                                                         <  1:4330529/150‑1 (MQ=255)
                               gacagggtggtggggatactcggggataccgagctttaccacgcgctgctcggcaagaaccacggttgaagcaattggggccgcttcgcggcccttcgcgggcacgcccgctcccacagGTTTGCTGAGATCTCTGTGGGAGCGGgcgcg                                                                                                                          >  2:1853807‑M2/120‑148 (MQ=255)
                                    tttGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTGGGCATGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                                  <  2:577815/150‑1 (MQ=21)
                                    tttGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCtg                                                                                                                                        <  1:289442/128‑1 (MQ=255)
                                    tttGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCtg                                                                                                                                        >  2:289442/1‑128 (MQ=255)
                                      gggCTTTGCGGGGCGGATGAAACGGTGGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTGCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATGGCGCGGTCCCTGTGGGTGCGGGCGTGCCCGCGAaga                                                                                                                <  2:3671973/149‑1 (MQ=38)
                                          tttGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGa                                                                                                                                                       >  1:3169514/1‑108 (MQ=255)
                                          tttGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGa                                                                                                                                                       <  2:3169514/108‑1 (MQ=255)
                                          tttGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGcc                                                                                                            >  2:4908501/1‑150 (MQ=255)
                                               gTGGCGGATGAACCGGTTGTTGGGTACGGCGGGGGTTCTTACATAAGCTGCTGTTTCGGCCTGAGTTCCGGGCCTCTTGGTGGGCACGCCCGCTCCCACTGGGACCGCGCAATCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGcccggaa                                                                                                       <  1:3114679/150‑5 (MQ=17)
                                                tGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGTGGCCCTAAAg                                                                                                      >  2:5121936/1‑150 (MQ=255)
                                                   cGGATGAAACGGTTGTTGGTTTCGGCGGGGGTTCTTACAGATGCGGCTGTTTCGGCCTGACTCCGGGGCCTCTTCGGGGCCACGCCCGCTCCCTTAGTTTTGCGCGTTCCCTGTGGGTGGGGGCGGGCCCGCGAAGAGGCCCTAAAGGac                                                                                                   <  2:1348515/150‑1 (MQ=255)
                                                          aaCGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGGACGCCCCCACCCACAGGATTTCGCGGGCCCTGTGGGAGCGGGCGGGCCCGCCTACACGGCCGAAAGGACCCCCAcc                                                                                            >  1:388767/1‑150 (MQ=255)
                                                          aaCGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTATTAGCGGGCAAGCCCGCTCCCACAGGCTTGGGCGGTCCCAGTGGGAGCCGGCGTGCCCCCGAATAGGCCCAAAATGACACCCGcc                                                                                            >  1:3085861/1‑150 (MQ=255)
                                                                 tttaccacgcgctgctcggcaagaaccacggttgaagcaattggggccgcttcgcggcccttcgcgggcacgcccgctcccacagGTTTGCTGAGATCTCTGTGGGAGc                                                                                                                                 <  2:4866498‑M2/24‑1 (MQ=255)
                                                                 tttaccacgcgctgctcggcaagaaccacggttgaagcaattggggccgcttcgcggcccttcgcgggcacgcccgctcccacagGTTTGCTGAGATCTCTGTGGGAGc                                                                                                                                 >  1:4866498‑M2/86‑109 (MQ=255)
                                                                    gggTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGt                                                                                                                           >  1:1814876/1‑109 (MQ=255)
                                                                    gggTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGt                                                                                                                           <  2:1814876/109‑1 (MQ=255)
                                                                              ggggTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGCCACCCTCCAACTGGCGGCCCGGCGCCTg                                                                        >  2:3994383/1‑150 (MQ=255)
                                                                                    caagaaccacggttgaagcaattggggccgcttcgcggcccttcgcgggcacgcccgctcccacagGTTTGCTGAGATCTCTGTGGGa                                                                                                                                   >  1:965546‑M2/67‑88 (MQ=255)
                                                                                    caagaaccacggttgaagcaattggggccgcttcgcggcccttcgcgggcacgcccgctcccacagGTTTGCTGAGATCTCTGTGGGa                                                                                                                                   <  2:965546‑M2/22‑1 (MQ=255)
                                                                                          gAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTgcg                                                                                                                                                 <  1:3676190/65‑1 (MQ=255)
                                                                                          gAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTgcg                                                                                                                                                 >  2:3676190/1‑65 (MQ=255)
                                                                                                  tgTTTCGGCCTGACTTCCGGGCCTCTTGGCGGGCACGCCCCCTCCCAAAGGATTGCGCGGTCCCGGTGGGCGCCGGCGTGCCCCCCAAGAGGCCCGAAAGGCCCCCCGCCAACAGGCGGCCCGGCGCCTGCAGCAAATACGCCATCaaag                                                    >  2:3691999/1‑148 (MQ=255)
                                                                                                                  ccGGGCCTCTTCGCGGGCACCCCCCCACCCACAGGATTGCGCCGTCCCGGTGGGAGCGGGCGGGCCCGCGCACGGGCCCGAAAGGGCACCCGCCCACTGGCGGCCCGGCGCCAGCCGCACATACGCCCGCAATTCCCCAATCGGGAGCgg                                    >  2:2882921/1‑150 (MQ=255)
                                                                                                                  ccGGGCCTCGTCGCGGGCACGCCCGCTCCCACAGGGTTGCGCGGTCCCGGTGGGAGCGGGGGTGCCCGCGAAGAGGCCCGAAAGGACACCCACCCACGGGCGGCCCGGCGCCGGCCGCCACTACGCCCGCACTTCCCCGATGGGGCGCgg                                    >  2:2882964/1‑150 (MQ=16)
                                                                                                                   ggggCCTCTTCGGGGGCACGCCCGCTCCCACTTTTTTGCGCGGTCCGTGTGGGTGGGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGt                                   <  1:3273141/149‑1 (MQ=255)
                                                                                                                            ttcgcgggcacgcccgctcccacaggGACCGCGCAATCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGAc                           >  2:1305100‑M2/27‑150 (MQ=255)
                                                                                                                                     cacgcccgctcccacagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTCTCGGTAGCGGTCTGCTGACCAGGAAGc                   >  2:3271008‑M2/18‑150 (MQ=255)
                                                                                                                                     cacgcccgctcccacagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGc                   >  1:4675073‑M2/18‑150 (MQ=255)
                                                                                                                                     acgcccgctcccacaggGACCGCGCAATCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCCGGCGCCTGCAGCAAATACGCCATCAATTCCCCTCTCGGTAGCGGTCTGCTGACCAGGAAGcc                  >  1:3339311‑M2/18‑150 (MQ=255)
                                                                                                                                      cgcccgctcccacaggGACCGCGCAATCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCt                 <  1:1305100‑M2/134‑1 (MQ=255)
                                                                                                                                          cgctcccacaggGACCGCGCAATCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCCGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCa             >  1:4544623‑M2/13‑150 (MQ=255)
                                                                                                                                            gctcccacagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCAc            >  2:1037342‑M2/11‑150 (MQ=255)
                                                                                                                                              tcccactgGCTTGCGCGGTCCGTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCt          <  1:3959556‑M2/142‑1 (MQ=255)
                                                                                                                                              tcccacagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGcc                                                           >  2:2185352‑M2/9‑101 (MQ=255)
                                                                                                                                              tcccacagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGa                      <  1:5193358‑M2/130‑1 (MQ=255)
                                                                                                                                              tcccacagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGa                      >  2:5193358‑M2/9‑138 (MQ=255)
                                                                                                                                              tcccacagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAAGACGCCAGCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCt          >  1:3916690‑M2/9‑150 (MQ=255)
                                                                                                                                              tcccacagGATTGCGCGGTCCCGGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGcc                                                           <  1:2185352‑M2/93‑1 (MQ=255)
                                                                                                                                                 cacagGATTGCGCGGTGCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCCTCAATTCCCCTATCGGTAGCGGTCTGCTGACCCGGAAGCCTTGCACCTGGt       >  1:2417620‑M2/6‑150 (MQ=255)
                                                                                                                                                 cacagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTc                                                                                             >  2:3701831‑M2/6‑64 (MQ=255)
                                                                                                                                                 cacagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTc                                                                                             <  1:3701831‑M2/59‑1 (MQ=255)
                                                                                                                                                 cacagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGTCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGt       >  1:4943299‑M2/6‑150 (MQ=255)
                                                                                                                                                 cacagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCCGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGt       >  2:2093955‑M2/6‑150 (MQ=255)
                                                                                                                                                 cacagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGt       >  1:2051357‑M2/6‑150 (MQ=255)
                                                                                                                                                 cacagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGt       >  2:5279106‑M2/6‑150 (MQ=255)
                                                                                                                                                 cacagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGt       >  2:4627226‑M2/6‑150 (MQ=255)
                                                                                                                                                 cacagGATTGCGCGGTCCCTGGGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGt       >  2:522881‑M2/6‑150 (MQ=255)
                                                                                                                                                  acagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTc      >  1:838570‑M2/5‑150 (MQ=255)
                                                                                                                                                  acagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTc      >  1:636689‑M2/5‑150 (MQ=255)
                                                                                                                                                     gGATGGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCa   <  1:4908501‑M2/149‑1 (MQ=255)
                                                                                                                                                      gATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTgcag                                                                     >  2:656469/1‑83 (MQ=255)
                                                                                                                                                      gATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTgcag                                                                     <  1:656469/83‑1 (MQ=255)
                                                                                                                                                      gATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCAg  <  1:1037342/150‑1 (MQ=255)
                                                                                                                                                      gATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCCGCTGACCAGGAAGCCTTGCACCTGGTCGCAg  >  2:2542403/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GCGCTGGGGCTGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCAG  >  NC_002947/6010971‑6011268

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: