Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,707,382 T→A 17.1% intergenic (‑152/+25) kdpF ← / ← PP_4162 potassium‑transporting ATPase subunit F/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,707,3820TA17.1% 386.8 / 67.7 158intergenic (‑152/+25)kdpF/PP_4162potassium‑transporting ATPase subunit F/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (71/60);  new base A (15/12);  total (86/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

Only 100 of 162 total aligned reads displayed.

GTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAAC  >  NC_002947/4707233‑4707530
                                                                                                                                                     |                                                                                                                                                    
gtgGAGCGCCCCTTGTTACGGCTTGGGCCGGTTTGGGGAACTTGTCCGAGGGCGCGGCGTACCGGGGCGCGTGCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCtt                                                                                                                                                      <  2:3553783/150‑1 (MQ=255)
    tcctcCCCGTGTTCGGGCTTGGGCGGGTTGGGGGAATTGGTCGGCGGGCGCGGCGTCCGGGGGCGTGCGGGGGGGGGGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCTTTCGCGGGCAAGccc                                                                                                                                                  <  1:702992/148‑1 (MQ=255)
    aTCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGccc                                                                                                                                                  >  1:2998434/1‑150 (MQ=255)
      cgcccgGTGTTCGGGCTTGGGCCGCTTGGGGGAAATGGGCCGCGGGCGCGGCGTAAGGGGCCGCGCGCGGGGGGGGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCAAGCCcgc                                                                                                                                                <  1:3279044/148‑1 (MQ=255)
       tcCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCcgcg                                                                                                                                               >  1:1539325/1‑150 (MQ=255)
         cccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                             >  2:4828526/1‑150 (MQ=255)
          ccGTGTTACGGCTTGGGCGGGTTTGGGGAATTTGTCCGAGGGTGAGGGGTTATGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaa                                                                                                                                            <  2:91/150‑1 (MQ=255)
          ccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaa                                                                                                                                            >  1:3291691/1‑150 (MQ=255)
          ccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaa                                                                                                                                            >  2:4395612/1‑150 (MQ=255)
          ccGTGTTACGGCCGGGGCCGCTTTGGGGAATTTGTCCGCGGGAGCGGCGTCCCGGGGGGAGGGGGGGGGGGGGGGGGGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaa                                                                                                                                            <  2:1216409/150‑1 (MQ=255)
           ggtgtTCGGGGTTGGGCCGGTTTGGGGAATTGGTCCGCGGGAGCGGCGTACGGGGGCGCGCGGGGGGGGGGGGGGGGGGCATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaga                                                                                                                                           <  2:238484/149‑3 (MQ=255)
           cGTTTTCCGGCTGGGGCGTCTTGGGGGACTGTGTCCGCGGGTGCGGCGTACGGGGGGGCGTGGGGGGGGTGGGGGGGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAg                                                                                                                                           <  2:4558508/150‑1 (MQ=255)
           cGTGTTAGGGGTTGGGCCGGTTTGGGGAATGTGGCCGTGGGTGCGGGGTACCGGGGCGCGCTCGGGGGGGGGGGGGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAgg                                                                                                                                          <  2:2927186/150‑1 (MQ=255)
           cGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAg                                                                                                                                           >  1:4625151/1‑150 (MQ=255)
            gggtcccGGCTTGGGCCGCTTTGGGCAACTGGCCCGCGGGCGCGGCGCCCCGGGGCGAGCGCGGGGGGTGGGGGGGGCCATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAgg                                                                                                                                          <  1:389721/144‑1 (MQ=255)
            ggggtccGGCTGGGGCCGCTTGGGGGTATTTGTCCGCGGGCGAGGCGTAAGGGGGCGAGCTCGGGGGGTGGGGGTGGGCATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAgg                                                                                                                                          <  2:4194344/144‑1 (MQ=255)
               ttcgggCTTGGGCGGGTTGGGGTAATGTGCCGGTGGGCGCGGCGTACCGGGGCGCCCTGGGGGGGGGGGGGTGGCTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGcc                                                                                                                                       <  1:709235/146‑1 (MQ=12)
               ttACGGCTCGGGCCGCTTTGGGGAATGTGTCCGCGGGCGCGGCTTCCGGGGGCGAGCGCGGGGGGGGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGcc                                                                                                                                       <  2:242927/150‑1 (MQ=21)
                 agggCTGGGGCCGGTTTGGGTATATGGTCGGAGGGCGAGGCGTAACGGGGCGTGTTGGGGGGGTGGGGGGGCGAATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGt                                                                                                                                     <  2:4625151/148‑1 (MQ=255)
                     ctggggCCGCTGTGGGCACCTGGCCCGCGGGCGCGGCGTCACGGGGCGCGCTCGGGGGGGGCGGGCGCCTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACa                                                                                                                                 <  1:1928668/147‑1 (MQ=255)
                             gCTTTCGCGAATTTGCCCGAGGCCGTTGCGTACCGGGGCGAGATCGGGGGTGGGGGGCGCCTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTAtt                                                                                                                         <  2:716974/150‑1 (MQ=255)
                              cTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTAttg                                                                                                                        >  1:881965/1‑150 (MQ=255)
                               cttGGGTAACTGGGCCGCGGGCGCGGCGCAACGGGGCGCGAGCGGGGGGGGCGGGCCGGTATAAAGAAGGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTAttgt                                                                                                                       <  1:3996248/149‑1 (MQ=255)
                                  ggggAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATtgttgt                                                                                                                    >  2:1006460/1‑150 (MQ=255)
                                       aaTTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTGGCGGGCTTGGCCGCGGAAGGGCCGTTTGAGGGTTTTTGTTGTggggg                                                                                                               >  1:250171/1‑150 (MQ=255)
                                       aaTTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTggcgg                                                                                                               >  2:3241302/1‑150 (MQ=255)
                                            gccGAGGGAGAGGCGTAAGGGGGCGCGCGCGGGGGGGGGGGGTGGGAATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGgctgct                                                                                                          <  2:2032966/149‑1 (MQ=255)
                                               gAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAAAGGCCACATCGCGGGCTTTCCCCCGAAAGGGCCGTTACAGGCATATTTTTGTTGCGGCTGCTgcg                                                                                                       >  1:1071483/1‑150 (MQ=255)
                                               gAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGGATAAAGAATGCGTCAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTCTTGTTGTGGCGGCTGCTgcg                                                                                                       >  2:1437474/1‑150 (MQ=255)
                                                   gtgagGCGTTCCGGGGGGAGTTGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGcgccgc                                                                                                   <  2:1297472/148‑1 (MQ=255)
                                                     gagGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATCCGTAAAAATCAGCGGTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCATGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCAGCGCCGCGt                                                                                                 >  2:2811537/1‑150 (MQ=255)
                                                        gCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTTTGGCGGCTGCTGCGCCGCGTTAc                                                                                              >  1:2972467/1‑150 (MQ=255)
                                                          gtatgggggCGAGTTCGGGGGGGGGGGGGGGGTATAAAGAATGCGTAAAACTCAGAGTTTTGGTCTGCCTGTATCGGCCTCTTCGCGTTCTTGCCCGCGAATGTGCCTTTACAGGATTTTTTTTGTGGCGTCTGCTGCGACGCGTTActg                                                                                            <  2:2105295/145‑3 (MQ=255)
                                                          gcaccggggCGTGCGCGGGGGGGGGGGGTGGCAATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACgg                                                                                            <  1:346385/145‑1 (MQ=255)
                                                          gTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTTACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTAAAGGCTTTTTGTTTTTGCGGCTGCTTGGCCGCGTTACgg                                                                                            >  1:2101030/1‑150 (MQ=255)
                                                          gTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTAGCGGGCTTGCCCGCGAAAGGGCCGTTTCAGGCTTTTTGTTGTGGCGGCTGCTGCGCCGCGTTACgg                                                                                            >  1:2101133/1‑150 (MQ=255)
                                                           tAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCGTATTTTTGTGGCGGCTGCTGCGCCCCGTTACGGc                                                                                           >  1:5030263/1‑150 (MQ=255)
                                                              cggggCGTGCGCGGGGGGGGGGGGTGGCTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGACCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtcc                                                                                        <  1:4736490/149‑1 (MQ=255)
                                                               ggggCGCGCGCGGGGGGGGGGGGTGGGAATAAAGAAGGCGGAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtcca                                                                                       <  2:726133/150‑1 (MQ=255)
                                                                   cgcccccGGGCGGCGGGGGTCGCTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCGGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcc                                                                                   <  1:1053105/144‑1 (MQ=255)
                                                                   cGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACAGCtccagtcc                                                                                   >  2:1357984/1‑150 (MQ=255)
                                                                    gagcgcGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcca                                                                                  <  2:3680294/145‑1 (MQ=255)
                                                                    gagaTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCGTCGCGGGCTTGCCCCCGCAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcca                                                                                  >  1:88038/1‑150 (MQ=255)
                                                                      gaTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaa                                                                                                                                            <  2:271861/90‑1 (MQ=255)
                                                                      gaTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaa                                                                                                                                            >  1:271861/1‑90 (MQ=255)
                                                                       ctgggggggTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCcagca                                                                               <  2:3935678/147‑1 (MQ=255)
                                                                          ggggggTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCcagcagca                                                                            >  1:3816832/1‑150 (MQ=255)
                                                                              ggtgggggtgggTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCt                                                                        >  2:4075098/1‑150 (MQ=255)
                                                                                     gggggaATAAAGAATGCGTAAAAATCACCGTTGGGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCt                                                                 <  2:357094/148‑1 (MQ=255)
                                                                                       gggtATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGACGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTAc                                                               <  2:4350664/150‑1 (MQ=255)
                                                                                        ggtATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTAcc                                                              <  2:3816832/150‑1 (MQ=255)
                                                                                            taAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCcagca                                                                               <  2:333475/129‑1 (MQ=255)
                                                                                             aaaGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCGTATTGTTTTTGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCAACCCCCac                                                         >  1:2550650/1‑149 (MQ=39)
                                                                                                gAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCaa                                                      >  2:4794396/1‑150 (MQ=255)
                                                                                                 aaTGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcca                                                                                  >  2:2648501/1‑121 (MQ=255)
                                                                                                 aaTGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcca                                                                                  <  1:2648501/121‑1 (MQ=255)
                                                                                                   tGCGTAAAAATCAGCGGTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGt                                                   <  1:2951810/150‑1 (MQ=255)
                                                                                                     cGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCg                                                 <  1:4075098/150‑1 (MQ=255)
                                                                                                       tAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGc                                               >  2:365110/1‑150 (MQ=255)
                                                                                                           aaTCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGa                                           <  2:635000/150‑1 (MQ=255)
                                                                                                                    tgtgGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGc                                  <  1:1439307/150‑1 (MQ=255)
                                                                                                                    tgtgGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGc                                  <  1:1439315/150‑1 (MQ=255)
                                                                                                                      caataagCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGc                                               >  1:2076558/7‑135 (MQ=25)
                                                                                                                      caataagCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGc                                               <  2:2076558/129‑1 (MQ=25)
                                                                                                                       ggTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCaa                               <  2:1566949/150‑1 (MQ=255)
                                                                                                                           tGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGaa                           <  2:1861264/150‑1 (MQ=255)
                                                                                                                           tGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGaa                           >  1:5219837/1‑150 (MQ=255)
                                                                                                                            gCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAGGGCCAATGaa                           >  2:1677148/1‑149 (MQ=38)
                                                                                                                            gCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGaa                           <  1:1677148/149‑1 (MQ=255)
                                                                                                                              cTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGcc                        >  1:3359690/1‑150 (MQ=255)
                                                                                                                              cTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGcc                        <  1:2844300/150‑1 (MQ=255)
                                                                                                                               tGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCt                                                                 >  1:4396554/1‑108 (MQ=255)
                                                                                                                               tGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGaaaa                                     >  2:5376641/1‑136 (MQ=255)
                                                                                                                               tGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGaaaa                                     <  1:5376641/136‑1 (MQ=255)
                                                                                                                                 tAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCaaa                     <  2:1420930/150‑1 (MQ=255)
                                                                                                                                  aaCGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAg                    <  2:4657178/150‑1 (MQ=255)
                                                                                                                                   aCGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAg                                                     <  1:2193448/116‑1 (MQ=37)
                                                                                                                                     ggccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGAt                                          <  2:3008254/119‑1 (MQ=38)
                                                                                                                                     ggccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGAt                                          >  1:3008254/7‑125 (MQ=38)
                                                                                                                                     ggccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgca                 >  1:1196236/7‑150 (MQ=38)
                                                                                                                                      gccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaa                >  1:768710/6‑150 (MQ=39)
                                                                                                                                      gccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaa                >  1:1978324/6‑150 (MQ=39)
                                                                                                                                       ccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcca                                                                                  >  1:5402756/5‑83 (MQ=21)
                                                                                                                                       ccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcca                                                                                  <  2:5402756/79‑1 (MQ=21)
                                                                                                                                       ccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaag               >  1:4722851/5‑150 (MQ=39)
                                                                                                                                        ctctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaagc              >  1:4951097/1‑150 (MQ=255)
                                                                                                                                        ctctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaagc              >  1:1598578/1‑150 (MQ=255)
                                                                                                                                        ctctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaagc              >  1:4344186/1‑150 (MQ=255)
                                                                                                                                        cctttCGCGGGCAAGCCCGCGAAGGGGCCGTTACAGGCTTATTGTTGTGTCGGCTGCTGCGCCGCGTTACGGCtccagtcca                                                                                  <  1:4165099/79‑1 (MQ=25)
                                                                                                                                        cctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaagc              >  1:844691/4‑150 (MQ=39)
                                                                                                                                        cctttCGCGGGCAAGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGaaa                                      <  1:1257470/123‑1 (MQ=255)
                                                                                                                                         tctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaagca             <  1:1344052/150‑1 (MQ=255)
                                                                                                                                             cgcgGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAgcc         <  2:4897856/150‑1 (MQ=255)
                                                                                                                                              gcgGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAgccg        >  2:2569560/1‑150 (MQ=255)
                                                                                                                                                gggCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAGGCCAATGAAGCCAAAGGCAAGCAAGccgccg     <  1:388225/150‑1 (MQ=255)
                                                                                                                                                 ggCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGccgccg     <  2:881965/150‑1 (MQ=255)
                                                                                                                                                  gCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGAt                                          <  2:3203310/112‑1 (MQ=255)
                                                                                                                                                  gCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGAt                                          >  1:3203310/1‑112 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAgg                                                                          <  2:612224/78‑1 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAgg                                                                          >  1:612224/1‑78 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAAC  >  NC_002947/4707233‑4707530

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