Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,342,055 G→C 32.9% intergenic (+74/‑257) ribAB‑II → / → PP_3814 3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,342,0550GC32.9% 140.1 / 85.2 91intergenic (+74/‑257)ribAB‑II/PP_38143,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (36/25);  new base C (20/10);  total (56/35)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGC  >  NC_002947/4341906‑4342198
                                                                                                                                                         |                                                                                                                                                
gAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGTGACACCACCg                                                                                                                                                  <  2:185837/150‑1 (MQ=255)
gAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGTGACACCACCg                                                                                                                                                  >  2:4978416/1‑150 (MQ=255)
  aaTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGAc                                                                                                                                                <  1:4146922/150‑1 (MQ=255)
       gcgGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCAtt                                                                                                                                           <  1:2172756/150‑1 (MQ=255)
          gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTcgc                                                                                                                                        <  1:2857142/150‑1 (MQ=255)
          gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTcgc                                                                                                                                        <  2:2688199/150‑1 (MQ=255)
           ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑TCGGTGGTGTCACCCCCGCATTcgcg                                                                                                                                      <  1:3631985/150‑1 (MQ=37)
           ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTcgcg                                                                                                                                      <  1:1376718/150‑1 (MQ=255)
            gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCACTTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCgg                                                                                                                                     <  2:4358302/150‑1 (MQ=255)
            gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                     <  1:3724075/150‑1 (MQ=25)
            gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCgg                                                                                                                                     <  1:3536372/150‑1 (MQ=255)
            gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCgg                                                                                                                                     <  1:3536475/150‑1 (MQ=255)
            gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCgg                                                                                                                                     <  2:4099128/150‑1 (MQ=255)
            gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCgg                                                                                                                                     <  2:1426633/150‑1 (MQ=255)
             cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCggg                                                                                                                                    <  1:973143/150‑1 (MQ=25)
                  gggCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAAc                                                                                                                               <  1:3635714/150‑1 (MQ=255)
                  gggCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAAc                                                                                                                               <  2:2136519/150‑1 (MQ=255)
                   ggCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAAcc                                                                                                                              <  2:3959904/150‑1 (MQ=25)
                     cTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGGTTAACCCCCCAATCG‑‑CGGGTGACACCACCCCCCCAATCGCGGGTAAAaccc                                                                                                                            >  2:4571179/1‑149 (MQ=255)
                        tGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTTTACCCGcga                                                                                                                         <  2:3961197/150‑3 (MQ=255)
                        tGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCCCCAATTCG‑‑GGGGGGACAACACCCCCCCCATCGCGGGGTAAAccccgc                                                                                                                          >  1:3636943/1‑147 (MQ=255)
                              ggtggtggtGCCGAGCACCCCGTGCGACCGGGGGTGGGTTGGCTGAGAACACACAAGCAATTTTTGTCGCGGGGCCCTGTGGGCGCGGTTTCCCCGGCTATTGGG‑‑TGGGTGACACCACCGACGAATTCGGGGTAAAACCCGCTCCTACAg                                                                                                                   <  2:3925790/147‑1 (MQ=255)
                                aGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTTCCCGCCAATGGG‑‑GGGGTGGCACCACCGACGCATTAGCGGGTAAAACCGCTCCGACAggg                                                                                                                 >  2:1232560/1‑150 (MQ=255)
                                       gCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGGTTTACCCCCGATGCG‑‑AGGGGGGCACCAACCACGCAAAAGCGGGTCAAACCCCACCAACACGGccccggg                                                                                                          >  2:199312/1‑146 (MQ=255)
                                       gCAGAGCATCCCGGTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                     >  2:89442/1‑123 (MQ=18)
                                       gCAGAGCATCCCGGTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCGGTAGGAGCGGGTTTACCCGCGAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                     <  1:89442/123‑1 (MQ=18)
                                                          ccGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGGGGCACCACCGACGCATTCGCGGGTAAACCCGCACCGACAGGGCCCGGGGTGCTTTTGTCAGGGCGGAg                                                                                       >  2:3588345/1‑150 (MQ=37)
                                                                        cTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCgg                                                                                                                                     >  2:761864/1‑90 (MQ=255)
                                                                        cTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCgg                                                                                                                                     <  1:761864/90‑1 (MQ=255)
                                                                          gacaaCACACCCGAACTTGATTTCGGGGGCCCCTTTTGGAGGGGGTTTCCCCGCGAAGGGG‑‑GGGGTGCCACCCCCGACGATTTGGGGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGccc                                                                       <  2:190057/147‑1 (MQ=255)
                                                                                             aTTTAGCGGGGCCCTGTAGGAGCGGGTTTAACCCCGAAAGAG‑‑AGGGTGAAACCACCCACGCATTCGCGGGTAAAACCCCACCAACAGGGCCCCTGGTGGATTTGACAGGGCGGAGTGATGGCCCGGAGCCCCAAACCGCTTTCCcaaaac                                                    >  2:2699114/1‑148 (MQ=19)
                                                                                                                      ggTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGata                               >  1:5418654/1‑150 (MQ=25)
                                                                                                                      ggTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATcc                            >  1:2806267/1‑150 (MQ=255)
                                                                                                                      ggTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATcc                            >  2:3926899/1‑150 (MQ=255)
                                                                                                                       gTTTACCCGCGAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATccc                           >  1:140123/1‑150 (MQ=25)
                                                                                                                       gTTTACCCGCGAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATacc                           >  1:3974091/1‑150 (MQ=21)
                                                                                                                       gTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCa                                                                     >  1:3451849/1‑108 (MQ=255)
                                                                                                                       gTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCa                                                                     <  2:3451849/108‑1 (MQ=255)
                                                                                                                        tttACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAGTCCCATAATATGATATCCCg                          >  2:3169737/1‑150 (MQ=255)
                                                                                                                           aCCCGCGAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAg                       >  1:565383/1‑150 (MQ=25)
                                                                                                                           aCCCGCGAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAg                       >  1:565431/1‑150 (MQ=25)
                                                                                                                            cccGCGATTGCG‑‑TGGGTGACACCACCGACGCTTTGGGGTTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa                      <  1:3588345/150‑1 (MQ=255)
                                                                                                                            cccGCGACTGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACCGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa                      >  2:2496703/1‑150 (MQ=255)
                                                                                                                            cccGCGAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa                      >  1:4224442/1‑150 (MQ=25)
                                                                                                                            cccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa                      >  1:1725614/1‑150 (MQ=255)
                                                                                                                            cccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa                      >  2:2496795/1‑150 (MQ=255)
                                                                                                                             ccGCGAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                                 >  2:1341186/1‑122 (MQ=21)
                                                                                                                             ccGCGAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                                 <  1:1341186/122‑1 (MQ=21)
                                                                                                                             ccGCGAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGata                               >  2:4624337/1‑140 (MQ=25)
                                                                                                                             ccGCGAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGata                               <  1:4624337/140‑1 (MQ=25)
                                                                                                                             ccGCGAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  2:25202/1‑150 (MQ=25)
                                                                                                                             ccGCGAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  1:4782699/1‑150 (MQ=25)
                                                                                                                             ccGCGAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATACCGTAGAc                     >  2:1724384/1‑150 (MQ=25)
                                                                                                                             ccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                                 >  1:1202751/1‑122 (MQ=255)
                                                                                                                             ccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                                 <  2:5075315/122‑1 (MQ=255)
                                                                                                                             ccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                                 <  2:3961119/122‑1 (MQ=255)
                                                                                                                             ccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                                 <  2:1202751/122‑1 (MQ=255)
                                                                                                                             ccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                                 <  2:1202551/122‑1 (MQ=255)
                                                                                                                             ccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                                 >  1:1202551/1‑122 (MQ=255)
                                                                                                                             ccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                                 >  1:5075315/1‑122 (MQ=255)
                                                                                                                             ccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                                 >  1:3961119/1‑122 (MQ=255)
                                                                                                                             ccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  2:3992653/1‑150 (MQ=255)
                                                                                                                             ccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  2:991272/1‑150 (MQ=255)
                                                                                                                             ccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  1:1523056/1‑150 (MQ=255)
                                                                                                                             ccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  2:4785952/1‑150 (MQ=255)
                                                                                                                             ccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  1:2242073/1‑150 (MQ=255)
                                                                                                                             ccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  1:2473374/1‑150 (MQ=255)
                                                                                                                             ccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  2:2014869/1‑150 (MQ=255)
                                                                                                                             ccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  1:4149317/1‑150 (MQ=255)
                                                                                                                             ccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  1:4456445/1‑150 (MQ=255)
                                                                                                                             ccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  2:2844371/1‑150 (MQ=255)
                                                                                                                             ccGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  1:5201030/1‑150 (MQ=255)
                                                                                                                              cgcgAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc                    >  2:726084/1‑150 (MQ=25)
                                                                                                                               gcgAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCtt                                                  >  1:3336702/1‑119 (MQ=18)
                                                                                                                               gcgAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCtt                                                  <  2:3336702/119‑1 (MQ=18)
                                                                                                                               gcgAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCtt                                                  >  1:3336701/1‑119 (MQ=18)
                                                                                                                               gcgAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCtt                                                  <  2:3336701/119‑1 (MQ=18)
                                                                                                                               gcgAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTa                        >  1:4972035/1‑145 (MQ=25)
                                                                                                                               gcgAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTa                        <  2:4972035/145‑1 (MQ=25)
                                                                                                                               gcgAATGCG‑‑TCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt                   >  1:540276/1‑150 (MQ=25)
                                                                                                                               gcgAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCTGTAGACCt                   >  2:4688124/1‑150 (MQ=255)
                                                                                                                               gcgAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt                   >  2:3784634/1‑150 (MQ=255)
                                                                                                                               gcgAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt                   >  1:1964256/1‑150 (MQ=255)
                                                                                                                               gcgAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt                   >  2:5345712/1‑150 (MQ=255)
                                                                                                                                 gAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTg                 >  2:4633580/1‑150 (MQ=255)
                                                                                                                                 gAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTg                 >  2:5416776/1‑150 (MQ=255)
                                                                                                                                    tGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGt              <  2:4149317/150‑1 (MQ=255)
                                                                                                                                       g‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                                 <  1:438965/112‑1 (MQ=255)
                                                                                                                                       g‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                                 >  2:438965/1‑112 (MQ=255)
                                                                                                                                              gTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAgg       >  1:5375539/1‑150 (MQ=255)
                                                                                                                                               tGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGa                <  1:1780916/140‑1 (MQ=255)
                                                                                                                                               tGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGa                >  2:1780916/1‑140 (MQ=255)
                                                                                                                                                   accaccGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGc  >  1:262439/1‑150 (MQ=255)
                                                                                                                                                         |                                                                                                                                                
GAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCG‑‑AGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGC  >  NC_002947/4341906‑4342198

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