Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,011,121 T→C 12.7% intergenic (+40/+42) 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,1210TC12.7% 224.2 / 18.5 88intergenic (+40/+42)norM/PP_5263multidrug resistance protein NorM/GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (32/43);  new base C (3/8);  total (37/51)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

CGCTGGGGCTGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCAGCC  >  NC_002947/6010972‑6011270
                                                                                                                                                       |                                                                                                                                                      
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:1317588/150‑1 (MQ=255)
              tgtgCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCtgt                                                                                                                                         <  1:3937275/150‑1 (MQ=255)
              tgtgCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCtgt                                                                                                                                         <  2:2713185/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:2680071/150‑1 (MQ=255)
                tgCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTgg                                                                                                                                       <  1:2556429/150‑1 (MQ=255)
                tgCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTgg                                                                                                                                       <  2:4329804/150‑1 (MQ=255)
                tgCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTgg                                                                                                                                       <  1:524266/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:4443332‑M2/22‑1 (MQ=255)
                     caggataacgacagggtggtggggatactcggggataccgagctttaccacgcgctgctcggcaagaaccacggttgaagcaattggggccgcttcgcggcccttcgcgggcacgcccgctcccacagGTTTGCTGAGATCTCTGTGGGa                                                                                                                                     <  1:988953‑M2/22‑1 (MQ=255)
                     caggataacgacagggtggtggggatactcggggataccgagctttaccacgcgctgctcggcaagaaccacggttgaagcaattggggccgcttcgcggcccttcgcgggcacgcccgctcccacagGTTTGCTGAGATCTCTGTGGGa                                                                                                                                     <  1:792974‑M2/22‑1 (MQ=255)
                      ggTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCCTGTGGGAGCgg                                                                                                                                 <  2:2764103/150‑1 (MQ=37)
                      ggTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCCTGTGGGAGCgg                                                                                                                                 <  2:2936308/150‑1 (MQ=37)
                      ggTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCCTGTGGGAGCgg                                                                                                                                 <  2:2464843/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                                                                                                                                          >  2:289442/1‑128 (MQ=255)
                                   tttGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCtg                                                                                                                                          <  1:289442/128‑1 (MQ=255)
                                     gggCTTTGCGGGGCGGATGAAACGGTGGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTGCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATGGCGCGGTCCCTGTGGGTGCGGGCGTGCCCGCGAaga                                                                                                                  <  2:3671973/149‑1 (MQ=38)
                                         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                                                                                                                                     <  2:965546‑M2/22‑1 (MQ=255)
                                                                                   caagaaccacggttgaagcaattggggccgcttcgcggcccttcgcgggcacgcccgctcccacagGTTTGCTGAGATCTCTGTGGGa                                                                                                                                     >  1:965546‑M2/67‑88 (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         >  2:4627226‑M2/6‑150 (MQ=255)
                                                                                                                                                cacagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGt         >  1:2051357‑M2/6‑150 (MQ=255)
                                                                                                                                                cacagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGt         >  2:5279106‑M2/6‑150 (MQ=255)
                                                                                                                                                cacagGATTGCGCGGTCCCTGGGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGt         >  2:522881‑M2/6‑150 (MQ=255)
                                                                                                                                                 acagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTc        >  1:636689‑M2/5‑150 (MQ=255)
                                                                                                                                                 acagGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTc        >  1:838570‑M2/5‑150 (MQ=255)
                                                                                                                                                    gGATGGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCa     <  1:4908501‑M2/149‑1 (MQ=255)
                                                                                                                                                     gATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTgcag                                                                       <  1:656469/83‑1 (MQ=255)
                                                                                                                                                     gATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTgcag                                                                       >  2:656469/1‑83 (MQ=255)
                                                                                                                                                     gATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCAg    <  1:1037342/150‑1 (MQ=255)
                                                                                                                                                     gATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCCGCTGACCAGGAAGCCTTGCACCTGGTCGCAg    >  2:2542403/1‑150 (MQ=255)
                                                                                                                                                       ttGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCAGcc  >  1:5435891/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
CGCTGGGGCTGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCAGCC  >  NC_002947/6010972‑6011270

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: