Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,148,799 G→T 5.8% T151K (ACG→AAG)  PP_1007 ← transmembrane sensor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,148,7990GT 5.8% 219.2 / 2.4 70T151K (ACG→AAG) PP_1007transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base G (29/36);  new base T (2/2);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

CCTCGGGGACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTGGCCA  >  NC_002947/1148653‑1148946
                                                                                                                                                    |                                                                                                                                                   
ccTCGGGGACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTaa                                                                                                                                                  >  1:1422874/1‑150 (MQ=255)
       gACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTc                                                                                                                                           <  1:841228/150‑1 (MQ=255)
          cTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGTCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCTTACGGTGGTCGGc                                                                                                                                        <  1:860472/150‑1 (MQ=255)
          cTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGc                                                                                                                                        >  2:347971/1‑150 (MQ=255)
          cTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGc                                                                                                                                        >  1:998323/1‑150 (MQ=255)
           tGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGcc                                                                                                                                       >  1:793611/1‑150 (MQ=255)
            ggAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCt                                                                                                                                      >  1:484102/1‑150 (MQ=255)
            ggAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAACAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCt                                                                                                                                      >  1:169480/1‑150 (MQ=255)
             gAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTg                                                                                                                                     >  1:881580/1‑150 (MQ=255)
              aaGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGc                                                                                                                                    <  2:1410960/150‑1 (MQ=255)
               aGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTTCTGGCGTTCGCCGGCGAATGCCGAATTGGTGTTGAACAGCACTTTTGCGCCCTCTTCCAGTTGCAT‑TCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCa                                                                                                                                   <  1:330521/150‑1 (MQ=255)
                gTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAg                                                                                                                                  >  2:1098395/1‑150 (MQ=255)
                  aaGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGgc                                                                                                                                <  1:1239562/150‑1 (MQ=255)
                  aaGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTCTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGgc                                                                                                                                >  1:1336030/1‑150 (MQ=255)
                     gcTCGCTCTGTAGCAGGCGGGCGCCATGGCGGTCGTCCTTGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAC‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCa                                                                                                                             >  2:173637/2‑150 (MQ=255)
                     ccTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCa                                                                                                                             >  2:203882/1‑150 (MQ=255)
                      cTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTCAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCAt                                                                                                                            <  1:882679/150‑1 (MQ=255)
                       tCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATg                                                                                                                           <  2:992143/150‑1 (MQ=255)
                            cTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAg                                                                                                                      >  1:925988/1‑150 (MQ=255)
                                   aGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCg                                                                                                               >  2:248683/1‑150 (MQ=255)
                                     gcgggcgACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAg                                                                                                                      <  1:465684/141‑1 (MQ=255)
                                     gcgggcgACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAg                                                                                                                      >  2:465684/1‑141 (MQ=255)
                                      cgggcgACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCCCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACg                                                                                                            <  2:695536/150‑1 (MQ=255)
                                       gggcgACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCGTCGCGT‑‑CAGCGTACTGGTTGTTGGAGATCGGCCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGAggg                                                                                                           >  1:168422/1‑150 (MQ=255)
                                        ggcgACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTCAGCAGCACTTTTGGGTCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGt                                                                                                          <  2:436123/150‑1 (MQ=255)
                                           gACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGaa                                                                                                       <  2:396173/150‑1 (MQ=255)
                                             tcTGGCGGCCGTCCTGGCGTTTGCTGGCGAATGCCTAATTGGTGTTGAGGAGCACTTTTGCGCCGTCTTCGAGTAGCGC‑GCGCGTGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATa                                                                                                     >  2:474521/2‑150 (MQ=255)
                                               tGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GAGCTGGCGCTCGCCGC‑CATGCATAATGTGGTCGGTCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGc                                                                                                   <  1:248683/150‑1 (MQ=255)
                                                    gcCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGAc                                                                                              <  1:347971/150‑1 (MQ=255)
                                                    gcCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGAc                                                                                              <  1:693395/150‑1 (MQ=255)
                                                            ggCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACcagcagca                                                                                      <  1:1052874/150‑1 (MQ=255)
                                                                   gCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGc                                                                               <  2:1360397/150‑1 (MQ=255)
                                                                     tGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCgg                                                                             >  1:795704/1‑150 (MQ=255)
                                                                      ggCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGt                                                                            <  1:726148/150‑1 (MQ=255)
                                                                       gCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTg                                                                           <  2:496709/150‑1 (MQ=255)
                                                                               cGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAgggg                                                                   >  1:185843/1‑150 (MQ=255)
                                                                                  aTTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGttt                                                                <  1:420518/150‑1 (MQ=255)
                                                                                      gtgtTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAg                                                            >  1:33203/1‑150 (MQ=255)
                                                                                       tgtTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGt                                                           <  2:1080691/150‑1 (MQ=255)
                                                                                       tgtTGAGCAGCACTTTTGCGCCGGCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGt                                                           <  1:1310157/150‑1 (MQ=255)
                                                                                             gcagcaCTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTgcgcgt                                                     <  2:1422874/150‑6 (MQ=255)
                                                                                             gcagcaCTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTgc                                                     <  2:333152/150‑1 (MQ=255)
                                                                                              cagcaCTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCa                                                             <  1:447539/141‑1 (MQ=255)
                                                                                              cagcaCTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCa                                                             >  2:447539/1‑141 (MQ=255)
                                                                                               agcaCTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTgcgc                                                   >  1:1252309/1‑150 (MQ=255)
                                                                                                   cTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGgc                                               <  2:795704/150‑1 (MQ=255)
                                                                                                    ttttGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGgcg                                              >  1:1222720/1‑150 (MQ=255)
                                                                                                             gTCTTCCAGTTGCAG‑GCACTGGCGCTCGCCGA‑CCACCGTAAGGTGGAACAACTCAATGACGACGGGCAGGTTACCGACGGTGAATAGCCCGACCCGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCa                                     >  2:703051/1‑150 (MQ=255)
                                                                                                              tCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCAt                                    <  2:925988/150‑1 (MQ=255)
                                                                                                                  ccAGTTGCAG‑GCACTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGaa                                                                                                       >  2:643270/1‑79 (MQ=255)
                                                                                                                  ccAGTTGCAG‑GCACTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGaa                                                                                                       <  1:643270/79‑1 (MQ=255)
                                                                                                                    aGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTg                              <  1:1098395/150‑1 (MQ=255)
                                                                                                                     gTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGc                             <  1:652572/150‑1 (MQ=255)
                                                                                                                          aG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAgg                                                                     <  2:1270316/105‑1 (MQ=255)
                                                                                                                          aG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAgg                                                                     >  1:1270316/1‑105 (MQ=255)
                                                                                                                             gcgctggcgctCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTtgctg                      >  1:781329/1‑150 (MQ=255)
                                                                                                                              cgctggcgctCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGTGCAGGCGCCCCAGCCATGACCTGCGTtgctgc                     <  1:474521/150‑1 (MQ=255)
                                                                                                                              cgctggcgctCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTtgctgc                     <  2:1252309/150‑1 (MQ=255)
                                                                                                                               gctggcgctCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTgctgct                    >  2:813610/1‑150 (MQ=255)
                                                                                                                                       tCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCAtg            >  2:1285394/1‑150 (MQ=255)
                                                                                                                                         gCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCAtgtg          >  2:198491/1‑150 (MQ=255)
                                                                                                                                           cGA‑CCACCTTACGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCg        >  2:1236705/1‑150 (MQ=255)
                                                                                                                                            gA‑CCACCTTACGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTg                                                          <  1:987738/99‑1 (MQ=255)
                                                                                                                                            gA‑CCACCTTACGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTg                                                          >  2:987738/1‑99 (MQ=255)
                                                                                                                                            gA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAACGTCCGGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGt       <  1:15001/150‑1 (MQ=255)
                                                                                                                                            gA‑CCACCGTAAGGTGGTCGGCATGCAAGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGt       <  2:793611/150‑1 (MQ=255)
                                                                                                                                             a‑ccaccGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTg      <  1:813610/150‑1 (MQ=255)
                                                                                                                                             a‑ccaccGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTg      <  2:185843/150‑1 (MQ=255)
                                                                                                                                                caccGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTGGc    <  1:1373545/150‑1 (MQ=255)
                                                                                                                                                  ccGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTGGCCa  <  1:173637/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCTCGGGGACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAG‑GCGCTGGCGCTCGCCGA‑CCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTGGCCA  >  NC_002947/1148653‑1148946

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: ‑