Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,725,042 T→C 7.3% intergenic (+27/+41) 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,0420TC 7.3% 193.2 / 4.4 68intergenic (+27/+41)PP_5661/PP_4183putative lipoprotein/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (34/29);  new base C (3/2);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CCGCCGAGGAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAACTCTCGACGG  >  NC_002947/4724897‑4725185
                                                                                                                                                   |                                                                                                                                                
ccgccgAGGAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTc                                                                                                                                               <  5:294961/149‑1 (MQ=255)
  gccgAGGAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCgg                                                                                                                                             >  5:207252/1‑149 (MQ=255)
          aaaCGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGttttttt                                                                                                                                     <  7:233789/149‑1 (MQ=255)
                          aaaCATCAGTCCGAAGCCGATCAGATGGCTGCCGAACAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGCCAGAA‑‑CCTGCCAAGCAGTAATTGCCTAAAACAAAAAAACCCCACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTc                                                                                                                     <  4:232449/149‑1 (MQ=255)
                            aGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGa                                                                                                                   >  4:14098/1‑149 (MQ=255)
                            aGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGa                                                                                                                   <  6:207252/149‑1 (MQ=255)
                                  ggCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAAGGAGTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTGGGGTTTTTTTGTGCCTGCGAACTTTCGATTaa                                                                                                               >  1:282556/1‑149 (MQ=255)
                                   gCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCg                                                                                                            >  8:92739/1‑149 (MQ=255)
                                          ccGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGAt                                                                                                                  <  8:255257/136‑1 (MQ=255)
                                          ccGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGAt                                                                                                                  >  7:255257/1‑136 (MQ=255)
                                           cGAGCAGAAGGCTGCCGAATCAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCATAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGtt                                                                                                    <  6:97359/149‑1 (MQ=255)
                                            gAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTg                                                                                                   >  1:52869/1‑149 (MQ=255)
                                             aGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGACAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACATCGTCGGGttttttt                                                                                                                                     <  5:348371/114‑1 (MQ=39)
                                             aGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGACAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACATCGTCGGGttttttt                                                                                                                                     >  6:348371/1‑114 (MQ=39)
                                                aGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGtttt                                                                                                                                        <  4:321938/108‑1 (MQ=255)
                                                aGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGtttt                                                                                                                                        >  3:321938/1‑108 (MQ=255)
                                                  gtaggagCCGGAGTAGCCGGAGCGGCTTCTTCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAaac                                                                                              >  1:249607/7‑148 (MQ=38)
                                                    ggCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGttt                                                                                                         <  4:2654/135‑1 (MQ=255)
                                                    ggCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGttt                                                                                                         >  3:2654/1‑135 (MQ=255)
                                                       tGCCGAAGAAGCCGCTCCGGCTACTCAGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATAGTTTCAGTTGGAAACAAcc                                                                                          >  3:254739/1‑147 (MQ=255)
                                                       agCCGGAGTAGCCGGAGCGGCTTCTTCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGa                                                                                        <  2:249607/148‑1 (MQ=255)
                                                              gaagCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTcc                                                                                 >  7:74135/1‑149 (MQ=255)
                                                                 gCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAg                                                                                       >  2:83723/1‑140 (MQ=255)
                                                                 gCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAg                                                                                       <  1:83723/140‑1 (MQ=255)
                                                                        cGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTACGa                                                                                                                   >  4:160271/1‑105 (MQ=255)
                                                                        cGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTACGa                                                                                                                   <  3:160271/105‑1 (MQ=255)
                                                                         ggCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTc                                                                     <  2:52869/149‑1 (MQ=255)
                                                                          gCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCCGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGTGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTcc                                                                    <  2:36545/149‑1 (MQ=255)
                                                                                    cTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTacac                                                          <  3:14098/149‑1 (MQ=255)
                                                                                      ccTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACtt                                                        >  7:224766/1‑149 (MQ=255)
                                                                                       ctactaCACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACttt                                                       <  3:129303/149‑1 (MQ=255)
                                                                                        tactacACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTACCTTGGCCGGCTCCGACAGTACACTTtcg                                                     >  3:264167/1‑149 (MQ=255)
                                                                                          ctacACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTtcgg                                                    >  3:276878/1‑149 (MQ=255)
                                                                                           tacACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTtcggt                                                   <  7:92739/149‑1 (MQ=255)
                                                                                           tacACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGAATGCGAAATTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCATTGGCCGCATCCTACAGAACACTTtcggt                                                   >  8:224060/1‑149 (MQ=255)
                                                                                            acacCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCg                                                                                         >  6:17615/1‑111 (MQ=255)
                                                                                            acacCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCg                                                                                         >  8:302592/1‑111 (MQ=255)
                                                                                                cGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGgtg                                              <  8:224766/149‑1 (MQ=255)
                                                                                                 ggCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGgtgt                                             <  2:282556/149‑1 (MQ=255)
                                                                                                 ggCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGgtgt                                             >  3:103380/1‑149 (MQ=255)
                                                                                                   cAGAA‑‑CCTGTCAAGCAGTAATTGGATAAAACTAAAAAACCCGACTATGTCGGTTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGTTCGgtgt                                             <  4:254739/147‑1 (MQ=255)
                                                                                                      aa‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTAAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAaca                                        >  3:65009/1‑149 (MQ=255)
                                                                                                           cTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCa                                     >  6:1699/1‑149 (MQ=255)
                                                                                                              ccAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATc                                   <  4:133476/148‑1 (MQ=255)
                                                                                                               cAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATca                                  >  2:291835/1‑148 (MQ=255)
                                                                                                                aaGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATcaca                                >  3:291535/1‑149 (MQ=255)
                                                                                                                     gTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAg                           >  6:269549/1‑149 (MQ=255)
                                                                                                                      tAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGt                          <  5:1699/149‑1 (MQ=255)
                                                                                                                            gCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTtc                                                      <  2:140122/115‑1 (MQ=255)
                                                                                                                            gCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTtc                                                      >  1:140122/1‑115 (MQ=255)
                                                                                                                            gCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTcggtcggt                                               >  8:284421/1‑122 (MQ=255)
                                                                                                                            gCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTcggtcggt                                               <  7:284421/122‑1 (MQ=255)
                                                                                                                              tAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCg                   >  3:242589/1‑148 (MQ=255)
                                                                                                                                 aaCAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGt                          <  4:70698/138‑1 (MQ=255)
                                                                                                                                 aaCAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGt                          >  3:70698/1‑138 (MQ=255)
                                                                                                                                    aaaaaaaCCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTtcggt                                                   >  1:141409/1‑110 (MQ=255)
                                                                                                                                    aaaaaaaCCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTtcggt                                                   <  2:141409/110‑1 (MQ=255)
                                                                                                                                    aaaaaaaCCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGaa            <  1:291835/149‑1 (MQ=255)
                                                                                                                                    aaaaaaaCCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGaa            <  8:74135/149‑1 (MQ=255)
                                                                                                                                    aaaaaaaCCCGACATCGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGt                          <  2:224288/135‑1 (MQ=255)
                                                                                                                                    aaaaaaaCCCGACATCGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGt                          >  1:224288/1‑135 (MQ=255)
                                                                                                                                      aaaaaCCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAAct          >  4:27909/1‑149 (MQ=255)
                                                                                                                                        aaaCCCGACATCGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGATTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATAGATctg         >  3:281829/1‑147 (MQ=255)
                                                                                                                                           cccGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAACTCTCGa     <  4:276878/149‑1 (MQ=255)
                                                                                                                                            ccGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAACTCTCGAc    >  4:165846/1‑149 (MQ=255)
                                                                                                                                             cgacgaTGTCGGGTTTTTTTGTTCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTCGGCCGGCTCCGACAGTACACTTTcggtcggt                                               <  1:213314/106‑1 (MQ=255)
                                                                                                                                             cgacgaTGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTcggtcggt                                               >  2:213314/1‑105 (MQ=255)
                                                                                                                                              gacgaTGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAACTCTCGACgg  >  2:141540/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCGCCGAGGAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAA‑‑CCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAACTCTCGACGG  >  NC_002947/4724897‑4725185

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: