Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,011,124:1 +T 12.9% intergenic (+43/+39) 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,1241.T12.9% 279.3 / 45.9 85intergenic (+43/+39)norM/PP_5263multidrug resistance protein NorM/GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (34/40);  new base T (3/8);  total (37/48)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

GGGGCTGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCAGCC  >  NC_002947/6010976‑6011270
                                                                                                                                                       |                                                                                                                                                  
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:2556429/150‑1 (MQ=255)
            tgCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTgg                                                                                                                                       <  1:2680071/150‑1 (MQ=255)
            tgCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTgg                                                                                                                                       <  1:524266/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:988953‑M2/22‑1 (MQ=255)
                 caggataacgacagggtggtggggatactcggggataccgagctttaccacgcgctgctcggcaagaaccacggttgaagcaattggggccgcttcgcggcccttcgcgggcacgcccgctcccacagGTTTGCTGAGATCTCTGTGGGa                                                                                                                                     <  1:792974‑M2/22‑1 (MQ=255)
                  ggTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCCTGTGGGAGCgg                                                                                                                                 <  2:2936308/150‑1 (MQ=37)
                  ggTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCCTGTGGGAGCgg                                                                                                                                 <  2:2764103/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                                                                                                                                   >  1:4866498‑M2/86‑109 (MQ=255)
                                                            tttaccacgcgctgctcggcaagaaccacggttgaagcaattggggccgcttcgcggcccttcgcgggcacgcccgctcccacagGTTTGCTGAGATCTCTGTGGGAGc                                                                                                                                   <  2:4866498‑M2/24‑1 (MQ=255)
                                                               gggTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGt                                                                                                                             <  2:1814876/109‑1 (MQ=255)
                                                               gggTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGt                                                                                                                             >  1:1814876/1‑109 (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)
                                                                                                                                                       |                                                                                                                                                  
GGGGCTGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCAGCC  >  NC_002947/6010976‑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: