Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,915,872 A→T 7.8% intergenic (+258/‑31) fklB‑II → / → PP_1715 FKBP‑type peptidyl‑prolyl cis‑trans isomerase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,915,8720AT 7.8% 183.1 / 4.4 62intergenic (+258/‑31)fklB‑II/PP_1715FKBP‑type peptidyl‑prolyl cis‑trans isomerase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base T (3/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGGATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCGCGCCTCGGC  >  NC_002947/1915724‑1916016
                                                                                                                                                      |                                                                                                                                                
tAGGATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAAt                                                                                                                                                 >  2:1648317/1‑150 (MQ=255)
  ggATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAAtgt                                                                                                                                               <  2:1888702/150‑1 (MQ=255)
  ggATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAAtgt                                                                                                                                               <  1:754237/150‑1 (MQ=255)
        ttCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATaa                                                                                                                                         <  2:1701570/150‑1 (MQ=255)
          cTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAAcc                                                                                                                                       <  2:1488412/150‑1 (MQ=255)
             aCGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCttt                                                                                                                                    >  1:1172762/1‑150 (MQ=255)
                   aTCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCctct                                                                                                                              <  2:1344461/150‑1 (MQ=255)
                    tCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCtctc                                                                                                                             <  1:1252470/150‑1 (MQ=255)
                       tttGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGa                                                                                                                          <  1:1351927/150‑1 (MQ=255)
                        ttGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAg                                                                                                                         <  2:368149/150‑1 (MQ=255)
                            gccgcAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTc                                                                                                                     <  1:716427/150‑1 (MQ=255)
                            gccgcAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTc                                                                                                                     <  1:657299/150‑1 (MQ=255)
                               gcAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCAt                                                                                                                  >  2:1721454/1‑150 (MQ=255)
                                 aGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGACCATGGGGTCTGAAACGTCGCCCTGTTCACAAGGGTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGa                                                                                                                <  2:913641/150‑1 (MQ=255)
                                   gTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGTTGAAACGTCGACCTGTTCAACAGAAGTGAT‑‑ATTGAAAACGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAgc                                                                                                              >  1:1868439/1‑150 (MQ=255)
                                   gTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAgc                                                                                                              <  1:1180333/150‑1 (MQ=255)
                                                  aaCGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACATTTATTGTGGAt                                                                                                                                           >  1:356727/1‑106 (MQ=255)
                                                  aaCGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACATTTATTGTGGAt                                                                                                                                           <  2:356727/106‑1 (MQ=255)
                                                  aaCGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAAtttt                                                                                               >  2:785361/1‑150 (MQ=255)
                                                  aaCGACTGTGTAAAAAACGCCCGATCTGACCGCGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAAtttt                                                                                               <  1:1372697/150‑1 (MQ=255)
                                                     gACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAACCGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCGTGGAGGCTCCATGAGCGCACCCTGGAATTTTGcc                                                                                            <  1:1095914/150‑1 (MQ=255)
                                                             aaaaaaCGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGAATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCt                                                                                    >  1:690043/1‑150 (MQ=255)
                                                                   cGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGCATTTTGCCCGTTTCCTGCCAtt                                                                              >  2:282498/1‑150 (MQ=255)
                                                                    gCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTg                                                                             >  2:1640506/1‑150 (MQ=255)
                                                                      ccGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGc                                                                           >  2:1575281/1‑150 (MQ=255)
                                                                       cGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGcc                                                                          <  1:1213149/150‑1 (MQ=255)
                                                                             gACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAgcgc                                                                    <  2:1172762/150‑1 (MQ=255)
                                                                              aCCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGcc                                                                   <  1:25079/150‑1 (MQ=255)
                                                                               ccGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAACGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCt                                                                  >  2:706842/1‑150 (MQ=255)
                                                                                cGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCAc                                                                                                          >  2:535500/1‑109 (MQ=255)
                                                                                cGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCAc                                                                                                          <  1:535500/109‑1 (MQ=255)
                                                                                  ccccaggtgcctgccAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCt                                                               >  1:1018282/8‑150 (MQ=255)
                                                                                      ccTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGc                                                           >  1:325596/1‑150 (MQ=255)
                                                                                       cTTGCCTGGCAGGCACCTGGGGGCGAAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGcc                                                          >  1:1684218/1‑150 (MQ=255)
                                                                                           ccTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCgg                                                      <  2:325596/150‑1 (MQ=255)
                                                                                            cTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACATTTATTGTGGAt                                                                                                                                           >  2:1515906/1‑64 (MQ=25)
                                                                                            cTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACATTTATTGTGGAt                                                                                                                                           <  1:1515906/64‑1 (MQ=37)
                                                                                            cTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTCCCTTTCTAATGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGt                                                     >  2:65047/1‑150 (MQ=255)
                                                                                             tGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGATATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCAATGAGCGCACACTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTc                                                    >  1:894052/1‑150 (MQ=255)
                                                                                                 aggcaCCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGttt                                                >  2:44654/1‑150 (MQ=255)
                                                                                                       cTGGGCGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGc                                          >  2:355446/1‑150 (MQ=255)
                                                                                                           gggCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCttg                                      <  2:1018282/150‑1 (MQ=255)
                                                                                                             gCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCttgtt                                    <  1:44654/150‑1 (MQ=255)
                                                                                                               ggAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATATCCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTtgttgt                                  <  1:1648317/150‑1 (MQ=255)
                                                                                                               ggAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTtgttgt                                  >  1:560698/1‑150 (MQ=255)
                                                                                                                gAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCATCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTgttgtt                                 >  1:238867/1‑150 (MQ=255)
                                                                                                                    cGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGATCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTcgcg                             <  1:626434/150‑1 (MQ=255)
                                                                                                                     gTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCgg                            <  1:1412235/150‑1 (MQ=255)
                                                                                                                            tGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGgcgcg                     >  1:657374/1‑150 (MQ=255)
                                                                                                                                  acaAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTcgcg                             <  1:1557500/137‑1 (MQ=255)
                                                                                                                                  acaAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTcgcg                             >  2:1557500/1‑137 (MQ=255)
                                                                                                                                    aaGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCttg                                      <  1:517260/126‑1 (MQ=255)
                                                                                                                                    aaGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCttg                                      >  2:517260/1‑126 (MQ=255)
                                                                                                                                      aGGTTATCCACAATAGATTTGAATACCCTTTTCCTCTCGGAGGCTCCAGGAGCGCACACTGGAATTTTGCACGTTCCCTGCCATCGCCCGATTGTCTGCGCATCCGCGGTCTTTTCCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGcc             <  2:690043/150‑1 (MQ=255)
                                                                                                                                         ttATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCcgc          <  1:1575281/150‑1 (MQ=255)
                                                                                                                                            tCCACATTTATTGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTTGCCGAGCGCCGGCTCAGCCGCGGTCGGTTGCCGGCCACGTCGTTCGCGGTGGCGCACAAAGCCCCGCGcc       >  2:1595473/1‑150 (MQ=255)
                                                                                                                                            tCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCGCGcc       >  1:1730928/1‑150 (MQ=255)
                                                                                                                                              cacaATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCGCGCCTc     >  1:322395/1‑150 (MQ=255)
                                                                                                                                              cacaATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCGCGCCTc     >  2:1056025/1‑150 (MQ=255)
                                                                                                                                               acaATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCGCGCCTCg    >  1:395959/1‑150 (MQ=255)
                                                                                                                                               acaATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCGCGCCTCg    >  1:437537/1‑150 (MQ=255)
                                                                                                                                                 aaTAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATTAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCGCGCCTCGGc  >  2:596463/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
TAGGATGTTTCTGACGCAAATCGTTTGCGCCGCAGGTAGTCATGGCGTGAAACGACTGTGTAAAAAACGCCCGATCTGACCGGGCCCCTTGCCTGGCAGGCACCTGGGGGCGGAAACGTCGCCCTGTTCACAAGGTTATCCACAATAAATGTGGATAACCTTTTCCTCTCGGAGGCTCCATGAGCGCACCCTGGAATTTTGCCCGTTTCCTGCCATTGGCCGAGCGCCTGCTCAGCCGCGGTCGTTTGCCGGCCTTGTTGTTCGCGGTGGCGCGCAAAGGCCCGCGCCTCGGC  >  NC_002947/1915724‑1916016

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: