Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,011,122 T→C 12.6% intergenic (+41/+41) 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,1220TC12.6% 226.4 / 18.4 88intergenic (+41/+41)norM/PP_5263multidrug resistance protein NorM/GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (34/41);  new base C (3/8);  total (37/51)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGGGCTGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCAGCC  >  NC_002947/6010973‑6011270
                                                                                                                                                       |                                                                                                                                                     
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                                                                                                                                         <  2:2713185/150‑1 (MQ=255)
             tgtgCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCtgt                                                                                                                                         <  1:1317588/150‑1 (MQ=255)
             tgtgCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCtgt                                                                                                                                         <  1:3937275/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:4443332‑M2/22‑1 (MQ=255)
                    caggataacgacagggtggtggggatactcggggataccgagctttaccacgcgctgctcggcaagaaccacggttgaagcaattggggccgcttcgcggcccttcgcgggcacgcccgctcccacagGTTTGCTGAGATCTCTGTGGGa                                                                                                                                     <  1:792974‑M2/22‑1 (MQ=255)
                    caggataacgacagggtggtggggatactcggggataccgagctttaccacgcgctgctcggcaagaaccacggttgaagcaattggggccgcttcgcggcccttcgcgggcacgcccgctcccacagGTTTGCTGAGATCTCTGTGGGa                                                                                                                                     <  1:988953‑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)
                                        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                                                                                                                                   >  1:4866498‑M2/86‑109 (MQ=255)
                                                               tttaccacgcgctgctcggcaagaaccacggttgaagcaattggggccgcttcgcggcccttcgcgggcacgcccgctcccacagGTTTGCTGAGATCTCTGTGGGAGc                                                                                                                                   <  2:4866498‑M2/24‑1 (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)
                                                                                                                                                       |                                                                                                                                                     
GCTGGGGCTGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCAGCC  >  NC_002947/6010973‑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: