Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,725,040 G→A 5.5% intergenic (+25/+43) PP_5661 → / ← PP_4183 putative lipoprotein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,725,0400GA 5.5% 331.2 / 4.0 107intergenic (+25/+43)PP_5661/PP_4183putative lipoprotein/membrane protein
Reads supporting (aligned to +/- strand):  ref base G (50/51);  new base A (4/2);  total (54/53)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 110 total aligned reads displayed.

GCCGAGGAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAACTCTCGACGGTAC  >  NC_002947/4724899‑4725188
                                                                                                                                             |                                                                                                                                                    
gccgAGGAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCgg                                                                                                                                               <  2:76316/149‑1 (MQ=255)
gccgAGGAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCgg                                                                                                                                               <  2:47662/149‑1 (MQ=255)
     ggAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGtttt                                                                                                                                          <  2:359507/149‑1 (MQ=255)
        aaaCGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGttttttt                                                                                                                                       >  5:304264/1‑149 (MQ=255)
          agtccGAAGCTTCCAACGTTAAGGAGGAAGCCTAGCAGAAGTCTGCCTAAGAAGCCGCTCCGGCTACGCGGGCACCTACTACACCGGCATAACATGCCAAGCAGTAATTGGCTAAAACTAAAAAACCCGAGGATGTCGGGTTTTTTtgt                                                                                                                                     <  6:1401/146‑1 (MQ=255)
                aaGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGATTTATTGTCCCTGc                                                                                                                               >  8:34565/1‑149 (MQ=255)
                    tGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACCCCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAATAACCCGACGATGTCGGGTTTGTTTGTGCCTGGGacg                                                                                                                           >  2:302120/1‑147 (MQ=255)
                    tGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGAAGTAATTGGCTAAAACAAAAAAAACCGACGATGTCGGGTTTTTTTCTGCCTGCGAAc                                                                                                                           >  7:218005/1‑149 (MQ=255)
                      ccAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACtt                                                                                                                         >  4:193507/1‑149 (MQ=255)
                       cAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACttt                                                                                                                        >  3:15400/1‑149 (MQ=255)
                         aaGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTc                                                                                                                       >  6:326592/1‑148 (MQ=255)
                         aaGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCg                                                                                                                      <  3:36979/149‑1 (MQ=255)
                          aGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGa                                                                                                                     <  5:106205/149‑1 (MQ=255)
                               aGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACATCGTCGGGttttttt                                                                                                                                       <  3:367901/126‑1 (MQ=255)
                               aGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACATCGTCGGGttttttt                                                                                                                                       >  4:367901/1‑126 (MQ=255)
                               aGGCCGAACCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAAt                                                                                                                <  1:390026/149‑1 (MQ=255)
                                ggCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATc                                                                                                               >  4:325699/1‑149 (MQ=255)
                                    gAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGttt                                                                                                           >  7:313128/1‑149 (MQ=255)
                                       gCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAg                                                                                                        <  4:17438/149‑1 (MQ=255)
                                         cGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGtt                                                                                                      <  3:325699/149‑1 (MQ=255)
                                          gAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACttt                                                                                                                        <  5:164930/130‑1 (MQ=255)
                                          gAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACttt                                                                                                                        >  6:164930/1‑130 (MQ=255)
                                          caGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTg                                                                                                     >  7:346712/2‑149 (MQ=255)
                                           aGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTgg                                                                                                    <  1:44454/149‑1 (MQ=255)
                                             cAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGa                                                                                                   >  5:389649/1‑148 (MQ=255)
                                                   gCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAAcc                                                                                            >  2:403686/1‑149 (MQ=255)
                                                      gCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGttttttt                                                                                                                                       <  8:151674/103‑1 (MQ=255)
                                                      gCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGttttttt                                                                                                                                       >  7:47382/1‑103 (MQ=255)
                                                      gCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGttttttt                                                                                                                                       >  7:151674/1‑103 (MQ=255)
                                                       ccGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAg                                                                                         <  6:262843/148‑1 (MQ=255)
                                                       ccGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGc                                                                                        >  2:382235/1‑149 (MQ=255)
                                                          aagaagCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTgg                                                                                                    >  4:280657/1‑134 (MQ=255)
                                                          aagaagCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTgg                                                                                                    <  3:280657/134‑1 (MQ=255)
                                                          aagaagCCGCTCCGGCTACTCCGGCTCCTACTACACCGACAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCttt                                                                                     >  3:168896/1‑149 (MQ=255)
                                                            gaagCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTcc                                                                                   <  5:99056/149‑1 (MQ=255)
                                                               gCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTg                                                                                <  3:202181/149‑1 (MQ=255)
                                                               gCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTg                                                                                >  4:202181/1‑149 (MQ=255)
                                                                 cGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGaa                                                                                                  <  3:151251/129‑1 (MQ=255)
                                                                 cGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGaa                                                                                                  >  4:151251/1‑129 (MQ=255)
                                                                  gCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGcc                                                                             <  5:55330/149‑1 (MQ=255)
                                                                   cTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCg                                                                            <  1:180927/149‑1 (MQ=255)
                                                                      cGGCTACTCCCGCTCCAACTACCCCGCCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCt                                                                         <  4:261712/149‑1 (MQ=255)
                                                                            cTCCGGCTCCTACTACACCGGCAGAACCTGCCAATCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACa                                                                   <  4:226988/149‑1 (MQ=255)
                                                                            cTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGa                                                                     <  8:313128/147‑1 (MQ=255)
                                                                            cTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGAc                                                                    <  8:346712/148‑1 (MQ=255)
                                                                            cTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACa                                                                   <  4:328571/149‑1 (MQ=255)
                                                                              ccGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAg                                                                  >  3:187002/1‑148 (MQ=255)
                                                                               cGGCTCCTACTACACCGGCATAACCTGCCAAGCAGTAATTGTCTAAAACAAAAAAACCCTACGATGTCGGGTTTTTTTGTGCCTGAGAACTTTCGATTAAt                                                                                                                <  3:351118/101‑1 (MQ=255)
                                                                               cGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTTGCTAAAACAAAAAAACACGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAAt                                                                                                                >  4:351118/1‑101 (MQ=255)
                                                                               cGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGTCTCCGACAGt                                                                 >  4:4728/1‑148 (MQ=255)
                                                                               cGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGt                                                                 >  5:201497/1‑148 (MQ=255)
                                                                                      tactacACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTtc                                                         >  4:346632/1‑149 (MQ=255)
                                                                                      tactacACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACAATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTtc                                                         >  1:302889/1‑149 (MQ=255)
                                                                                         tacACGGGCAGAACCTGCCAAGCAGTAATTGGCTAAAAGAAAAAAACCCGACGATGTGGGGTTTTTTTGTGCCTGCGAACTTTTGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTtcggt                                                      <  8:218005/149‑1 (MQ=255)
                                                                                         tacACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTcc                                                                                   <  2:132197/120‑1 (MQ=255)
                                                                                         tacACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTcc                                                                                   >  1:132197/1‑120 (MQ=255)
                                                                                         tacACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTtcggt                                                      <  6:201497/149‑1 (MQ=255)
                                                                                            acCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCt                                                                         <  3:108888/127‑1 (MQ=255)
                                                                                            acCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCt                                                                         >  4:108888/1‑127 (MQ=255)
                                                                                            acCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTtcggtcgg                                                   >  5:197911/1‑149 (MQ=255)
                                                                                             ccGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTcggtcggt                                                  <  6:389649/149‑1 (MQ=255)
                                                                                              cGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCGGACAGTACACTTTCGGTCGgtg                                                 >  4:217948/1‑149 (MQ=255)
                                                                                                  aGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGaa                                             >  8:367366/1‑149 (MQ=255)
                                                                                                  aGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGaa                                             >  4:387164/1‑149 (MQ=255)
                                                                                                   gAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAac                                            <  5:75871/149‑1 (MQ=255)
                                                                                                      ccTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACAcc                                         <  5:326592/149‑1 (MQ=255)
                                                                                                       cTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCa                                        >  5:102264/1‑149 (MQ=255)
                                                                                                         gCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTc                                                                        <  1:303455/115‑1 (MQ=255)
                                                                                                         gCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTAGATTAATCGTTTCAGTTGGAAACAAACGAGCTTTCATTGGCCGGCTc                                                                        >  2:303455/1‑115 (MQ=255)
                                                                                                          ccAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATca                                     >  8:14341/1‑149 (MQ=255)
                                                                                                           cAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGa                                              >  2:183897/1‑139 (MQ=255)
                                                                                                           cAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGa                                              <  1:183897/139‑1 (MQ=255)
                                                                                                           cAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACAcc                                         <  4:154707/144‑1 (MQ=255)
                                                                                                           cAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACAcc                                         >  3:154707/1‑144 (MQ=255)
                                                                                                           cAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATca                                     >  4:47506/1‑148 (MQ=255)
                                                                                                            aaGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATcaca                                   >  5:2994/1‑149 (MQ=255)
                                                                                                              gCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATc                                 <  1:403686/149‑1 (MQ=255)
                                                                                                              gCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATc                                 <  3:47506/149‑1 (MQ=255)
                                                                                                               cAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCAGGCTCCGACAGTACACATTCGGTCGGTCTGAACACCATCACATc                                 >  2:59414/1‑148 (MQ=255)
                                                                                                                aGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGa                               <  6:197911/149‑1 (MQ=255)
                                                                                                                aGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGa                               <  3:193507/149‑1 (MQ=255)
                                                                                                                        gCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATcc                       >  5:143205/1‑149 (MQ=255)
                                                                                                                          tAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGgtgt                                                <  4:309779/122‑1 (MQ=255)
                                                                                                                          tAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGgtgt                                                >  3:309779/1‑122 (MQ=255)
                                                                                                                          tAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCAGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATc                                      >  8:164037/1‑132 (MQ=255)
                                                                                                                          tAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCAGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATc                                      <  7:164037/132‑1 (MQ=255)
                                                                                                                           aaaaCAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGaa                    >  1:102351/1‑149 (MQ=255)
                                                                                                                            aaaCAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAAt                   <  4:187002/149‑1 (MQ=255)
                                                                                                                            aaaCAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTATCGGTCGGTGTGAACACCATCACATCGAGTACATCCGaa                    <  7:34565/148‑1 (MQ=255)
                                                                                                                                aaaaaaaCCCGCCATCGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGCAACAACCGAGCTTTCCTTGGACGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTTCTTAGGAATGGaa               >  5:328616/1‑149 (MQ=255)
                                                                                                                                aaaaaaaCCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGaa               >  7:124197/1‑149 (MQ=255)
                                                                                                                                aaaaaaaCCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGaa               <  3:346632/149‑1 (MQ=255)
                                                                                                                                aaaaaaaCCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGaa               >  8:356165/1‑149 (MQ=255)
                                                                                                                                  aaaaaCCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAAct             >  2:196370/1‑149 (MQ=255)
                                                                                                                                  aaaaaCCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAAc              >  4:216311/1‑148 (MQ=255)
                                                                                                                                  aaaaaCCCGACATCGTCGGGGTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTAGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGTAAc              >  4:52368/1‑148 (MQ=255)
                                                                                                                                   aaaaCCCGACAATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAActc            <  2:302889/149‑1 (MQ=255)
                                                                                                                                           acgaTGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAACTCTCGACGGt    <  3:216311/149‑1 (MQ=255)
                                                                                                                                            cgaTGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAACTCTCGACGGTa   <  6:143205/149‑1 (MQ=255)
                                                                                                                                             gaTGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAACTCTCGACGGTAc  <  6:347219/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
GCCGAGGAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAACTCTCGACGGTAC  >  NC_002947/4724899‑4725188

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