Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,779,105 G→A 6.0% L241L (CTG→CTA)  PP_1587 → NhaP‑type Na+(K+)/H+ antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,779,1050GA 6.0% 195.1 / 2.2 67L241L (CTG→CTA) PP_1587NhaP‑type Na+(K+)/H+ antiporter
Reads supporting (aligned to +/- strand):  ref base G (32/30);  new base A (2/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

GCTGACGCCAGCTTCAGCTTCGTGTTGGTCGCACTGGGCGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTCAGACCAGTACCCGCCTGCTCAACCGCAGTGTCTGGTCGTTGC  >  NC_002947/1778959‑1779253
                                                                                                                                                  |                                                                                                                                                    
gCTGACGCCAGCTTCAGCTTCGTGTTGGTCGCACTGGGCGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTtgctgctg                                                                                                                                                   <  2:102074/150‑1 (MQ=255)
    aCGCCAGCTTCAGCTTCGTGTTGGTCGCACTGGGCGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGt                                                                                                                                               <  2:1618493/150‑1 (MQ=255)
     cGCCAGCTTCAGGTTCGTGTTGGTCGAACTGGGCGGCCTTGCGGTGGGCGTGGTGGTGAGTTGGTTGATCGGTGGCCTGCGTAGGTGGGTGATCACCGGGGGCTGGGACGATCCGGCTGCCCATGTGGTGTTCATGTTGCTGCTGCCGtt                                                                                                                                              <  2:746194/150‑1 (MQ=255)
      gCCAGCTTCAGCTTCGTGTTGGTCGCACTGGGCGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTc                                                                                                                                             >  2:774613/1‑150 (MQ=255)
        cagcttcagcttCGTGTTGGTCGCACTGGGCGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGc                                                                                                                                           >  1:295884/1‑150 (MQ=255)
           cttcagcttcGTGTTGGTCGCACTGGGCGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCCACGAGGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGc                                                                                                                                        <  1:1019620/150‑1 (MQ=255)
                 cttcGTGTTGGTCGCACTGGGCGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGt                                                                                                                                  <  2:1275560/150‑1 (MQ=255)
                 cttcGTGTTGGTCGCACTGGGCGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACACATGTGGTGTTCGTTTTTTTTCGGCCTTGCTCAGCCTACGt                                                                                                                                  >  2:1870672/1‑150 (MQ=255)
                     gtgtTGGTCGCACTGGGCGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGcc                                                                                                                                                 <  1:225384/131‑1 (MQ=255)
                     gtgtTGGTCGCACTGGGCGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTg                                                                                                                              >  2:1161347/1‑150 (MQ=255)
                     gtgtTGGTCGCACTGGACGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGcc                                                                                                                                                 >  2:225384/1‑131 (MQ=255)
                           gTCGCACTGGGCGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAg                                                                                                                        >  1:1712546/1‑150 (MQ=255)
                           gTCGCACTGGGCGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAg                                                                                                                        >  2:1592227/1‑150 (MQ=255)
                           gTCGCACTGGGCGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCAGCTGCCCTTCGCAGCCTACGTGCTGGCCGAg                                                                                                                        <  2:1193824/150‑1 (MQ=255)
                            tCGCACTGGGCGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGc                                                                                                                       <  2:295884/150‑1 (MQ=255)
                              gCACTGGGCGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGt                                                                                                                     >  2:289614/1‑150 (MQ=255)
                                aCTGGGCGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGttt                                                                                                                   <  2:1046281/150‑1 (MQ=255)
                                 cTGGGCGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTg                                                                                                                  <  2:23362/150‑1 (MQ=255)
                                   gggCGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTggg                                                                                                                <  1:165682/150‑1 (MQ=255)
                                    ggcggcCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGc                                                                                                               <  2:1712546/150‑1 (MQ=255)
                                          cTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCg                                                                                                         >  2:245591/1‑150 (MQ=255)
                                           ttGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCgg                                                                                                        <  2:1165662/150‑1 (MQ=255)
                                               ggTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTCACGTTCCCGTGATAGGTGCTGTCCGATCGTTTGGGCGCTTCGGGTAt                                                                                                    >  2:1772561/1‑150 (MQ=255)
                                                gTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTAtt                                                                                                   <  1:1870672/150‑1 (MQ=255)
                                                      gTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATtctctct                                                                                             >  1:668983/1‑149 (MQ=255)
                                                       tGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCgg                                                                                            >  2:636155/1‑150 (MQ=255)
                                                            cTGAGTGGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTg                                                                                       <  1:275047/150‑1 (MQ=255)
                                                                gttggttgATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTggcgg                                                                                   <  1:1772561/150‑1 (MQ=255)
                                                                    gttgATCGGTCCCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTACTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGcggccgg                                                                               >  2:1443456/1‑150 (MQ=255)
                                                                        aTCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCAGGCatg                                                                           >  2:108/1‑150 (MQ=255)
                                                                         tCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCatga                                                                          >  2:296931/1‑150 (MQ=255)
                                                                            gTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGc                                                                       >  1:695798/1‑150 (MQ=255)
                                                                             tCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCa                                                                      <  1:296931/150‑1 (MQ=255)
                                                                                 cTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGc                                                                  <  1:245591/150‑1 (MQ=255)
                                                                                   gcgcATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTg                                                                <  1:1161347/150‑1 (MQ=255)
                                                                                   gcgcATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTg                                                                >  1:1176022/1‑150 (MQ=255)
                                                                                    cgcATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCTTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTgg                                                               <  1:930910/150‑1 (MQ=255)
                                                                                          tgGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGc                                                                       <  2:949936/136‑1 (MQ=255)
                                                                                          tgGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGc                                                                       >  1:949936/1‑136 (MQ=255)
                                                                                          tgGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGAc                                                         <  1:186045/150‑1 (MQ=255)
                                                                                          tgGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTACTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCgg                                                                                                        >  2:1056594/1‑103 (MQ=255)
                                                                                          tgGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTACTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCgg                                                                                                        <  1:1056594/103‑1 (MQ=255)
                                                                                          tgGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTACTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGNGTATTCNCTCGGCAGTGGCNGCCGGCATGATGCAGAGCTGGCTCGAc                                                         <  1:1443456/150‑1 (MQ=255)
                                                                                           ggATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGc                                                                                                                       >  1:407441/1‑87 (MQ=255)
                                                                                           ggATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGc                                                                                                                       <  2:407441/87‑1 (MQ=255)
                                                                                            gatgatCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACct                                                       >  1:970507/1‑150 (MQ=255)
                                                                                                       ggggCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGCTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTc                                            >  2:376493/1‑150 (MQ=255)
                                                                                                       ggggCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGCTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTc                                            <  1:376493/150‑1 (MQ=255)
                                                                                                        gggCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTCa                                           <  1:774613/150‑1 (MQ=255)
                                                                                                         ggCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCgg                                                          >  2:257054/1‑133 (MQ=255)
                                                                                                         ggCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGa                                                          <  1:257054/134‑1 (MQ=255)
                                                                                                           cTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTCAGAc                                        <  2:1055452/150‑1 (MQ=255)
                                                                                                            tGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTCAGAcc                                       >  2:986718/1‑150 (MQ=255)
                                                                                                            tGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTCAGAcc                                       >  2:1055189/1‑150 (MQ=255)
                                                                                                            tGGGACGATCCGGATACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTCCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCCGAGCTGGCTCGACCTGCTGCCGCGCCAGAcc                                       >  2:562206/1‑150 (MQ=255)
                                                                                                                aCGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTCAGACCAGTa                                   <  1:1055189/150‑1 (MQ=255)
                                                                                                                aCGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTCAGACCAGTa                                   <  2:1176022/150‑1 (MQ=255)
                                                                                                                  gATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGCGCGTTTGGGGGTTTGAGCGCCTCATGTATCACTGTCGGGGGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTCAGACCAGTAcc                                 >  1:1406377/1‑150 (MQ=255)
                                                                                                                  gATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTCAGACCAGTAcc                                 >  2:1657599/1‑150 (MQ=255)
                                                                                                                        gCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCatgatg                                                                        >  2:996472/1‑105 (MQ=255)
                                                                                                                        gCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCatgatg                                                                        <  1:996472/105‑1 (MQ=255)
                                                                                                                           aCCCATGTGGTGTTCATGTTGCTGCTGGCGTTCGAAGAATAAGTGCAGGGCGAGTGGTTGGGCGTTTCGGGTGTTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCAACCTGCGGCCGGGTCAGACCAGTACCCGCCTGCTc                        >  1:1306622/1‑150 (MQ=255)
                                                                                                                             ccATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTCAGACCAGTACCCGCCTGCTCaa                      >  1:1598097/1‑150 (MQ=255)
                                                                                                                                tgtgGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTCAGACCAGTACCCGCCTGCTCAACCg                   >  2:1646988/1‑150 (MQ=255)
                                                                                                                                    gtgtTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTCAGACCAGTACCCGCCTGCTCAACCGCAgt               <  1:636155/150‑1 (MQ=255)
                                                                                                                                    gtgtTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTCGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTCAGACCAGTACCCGCCTGCTCAACCGCAgt               >  1:630274/1‑150 (MQ=255)
                                                                                                                                       ttCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGAGGTTCTGGCTCGACCTGCTGCCGCGTCAGACCAGTACCCGCCTGCTCAACCGCAGTGTc            >  1:1477985/1‑150 (MQ=255)
                                                                                                                                          aTGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTCAGACCAGTACCCGCCTGCTCAACCGCAGTGTCTgg         >  1:275808/1‑150 (MQ=255)
                                                                                                                                           tgttgctgttgcGGCTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTTTCGGCAGTGGCGGCCGGCATCATGCAGAGCTGGCTCGCCCTGCTGCGGAGTCACTCCAGTACCCTCCTGCTCAACCGCAGTGTCTGGt        <  1:1592227/150‑1 (MQ=255)
                                                                                                                                               gctgctgcCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTCAGACCAGTACCCGCCTGCTCAACCGCAGTGTCTGGTCGtt    <  2:1306622/150‑1 (MQ=255)
                                                                                                                                                ctgctgcCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTCAGACCAGTACCCGCCTGCTCAACCGCAGTGTCTGGTCGtt    <  2:679700/149‑1 (MQ=255)
                                                                                                                                                ctgctgcCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTCAGACCAGTACCCGCCTGCTCAACCGCAGTGTCTGGTCGtt    >  1:679700/1‑149 (MQ=255)
                                                                                                                                                 tgctgcCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTCAGACCAGTACCCGCCTGCTCAACCGCAGTGTCTGGTCGTTGc  >  2:1148857/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCTGACGCCAGCTTCAGCTTCGTGTTGGTCGCACTGGGCGGCCTTGCGGTGGGCGTGGTACTGAGTTGGTTGATCGGTCGCCTGCGCATGTGGATGATCACCCGGGGCTGGGACGATCCGGCTACCCATGTGGTGTTCATGTTGCTGCTGCCGTTCGCAGCCTACGTGCTGGCCGAGCGTTTGGGCGTTTCGGGTATTCTCTCGGCAGTGGCGGCCGGCATGATGCAGAGCTGGCTCGACCTGCTGCCGCGTCAGACCAGTACCCGCCTGCTCAACCGCAGTGTCTGGTCGTTGC  >  NC_002947/1778959‑1779253

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