Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,033,047 C→G 7.0% intergenic (+212/‑198) PP_4435 → / → PP_4437 hypothetical protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,033,0470CG 7.0% 219.0 / 6.5 71intergenic (+212/‑198)PP_4435/PP_4437hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base C (31/35);  new base G (3/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

ACCGCTAAATCGATCCAGGAACAAACTCATGAACCACGTGCTTGCCGACATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGG  >  NC_002947/5032901‑5033187
                                                                                                                                                   |                                                                                                                                            
aCCGCTAAATCGATCCAGGAACAAACTCATGAACCACGTGCTTGCCGACATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCg                                                                                                                                            <  7:125950/149‑1 (MQ=255)
  cGCTAAATCGATCCAGGAACAAACTCATGAACCACGTGCTTGCCGACATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCCGGTCGGCGAc                                                                                                                                          >  7:150127/1‑149 (MQ=255)
  cGCTAAATCGATCCAGGAACAAACTCATGAACCACGTGCTTGCCGACATTAAAGAGTTTCAACAAGAGATCATCGATATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGAc                                                                                                                                          >  5:297022/1‑149 (MQ=255)
           gATCCAGGAACAAACTCATGAACCACGTGCTTGCCGACATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCCGGTCGGCGACCAGctgtga                                                                                                                                 <  8:150127/149‑5 (MQ=255)
               cAGGAACAAACTCATGAACCACGTGCTTGCCGACATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGt                                                                                                                             >  4:3935/1‑149 (MQ=255)
                aGGAACAAACTCATGAACCACGTGCTTGCCGACATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTg                                                                                                                            >  3:49299/1‑149 (MQ=255)
                aGGAACAAACTCATGAACCACGTGCTTGCCGACATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTg                                                                                                                            <  7:214565/149‑1 (MQ=255)
                    aCAAACTCATGAACCACGTGCTTGCCGACATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCGCTCAGTGGGGc                                                                                                                        >  5:145396/1‑149 (MQ=255)
                      aaaCTCATGAACCACGTGCTTGCCGACATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCtt                                                                                                                      <  4:44378/149‑1 (MQ=255)
                       aaCTCATGAACCACGTGCTTGCCGACATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTg                                                                                                                     <  4:49299/149‑1 (MQ=255)
                          tCATGAACCACGTGCTTGCCGACATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATAGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGgag                                                                                                                  <  7:269263/149‑1 (MQ=255)
                            aTGAACCACGTGCTTGCCGACATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGgaga                                                                                                                 <  4:173403/148‑1 (MQ=255)
                            aTGAACCACGTGCTTGCCGACATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCCGGTCGGCGACCAGct                                                                                                                                     >  8:174131/1‑128 (MQ=255)
                            aTGAACCACGTGCTTGCCGACATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCCGGTCGGCGACCAGct                                                                                                                                     <  7:174131/128‑1 (MQ=255)
                                  cACGTGCTTGCCGACATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCg                                                                                                          >  2:136263/1‑149 (MQ=255)
                                    cGTGCTTGCCGACATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGt                                                                                                                             >  5:282641/1‑128 (MQ=255)
                                    cGTGCTTGCCGACATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGt                                                                                                                             <  6:282641/128‑1 (MQ=255)
                                           gCCGACATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTgg                                                                                                  >  7:160365/1‑148 (MQ=255)
                                               aCATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCgg                                                                                                                                       >  7:243058/1‑107 (MQ=255)
                                               aCATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCgg                                                                                                                                       <  8:243058/107‑1 (MQ=255)
                                               aCATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGt                                                                                                                             >  2:297495/1‑117 (MQ=255)
                                               aCATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGt                                                                                                                             <  1:297495/117‑1 (MQ=255)
                                                   tAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCgg                                                                                                                                       <  5:382658/103‑1 (MQ=255)
                                                   tAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCg                                                                                         >  8:284344/1‑149 (MQ=255)
                                                   tAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCg                                                                                         >  5:237582/1‑149 (MQ=255)
                                                   tAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACGCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCg                                                                                         >  5:112045/1‑149 (MQ=255)
                                                   tAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATAGCACCTCACAGCTGGTCGCCGACCgg                                                                                                                                       >  6:382658/1‑103 (MQ=255)
                                                    aaaGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCgg                                                                                        >  5:277654/1‑149 (MQ=255)
                                                     aaGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCggg                                                                                       <  6:185009/149‑1 (MQ=255)
                                                     aaGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCggg                                                                                       <  6:356228/149‑1 (MQ=255)
                                                          tttCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGt                                                                                  >  2:132649/1‑149 (MQ=255)
                                                          tttCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGt                                                                                  >  2:207273/1‑149 (MQ=255)
                                                            tCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCGCTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGtt                                                                                 <  6:145396/148‑1 (MQ=255)
                                                            tCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTg                                                                                >  1:217591/1‑149 (MQ=255)
                                                             caacaaGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGtt                                                                                 >  6:354011/1‑147 (MQ=255)
                                                             caacaaGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGtt                                                                                 <  5:354011/147‑1 (MQ=255)
                                                              aacaaGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGa                                                                              <  4:335705/149‑1 (MQ=255)
                                                                caaGAGCTCATCGCTATCAGGCGCGACGATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTg                                                                                                   >  1:209237/1‑127 (MQ=255)
                                                                caaGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTg                                                                                                   <  2:209237/126‑1 (MQ=255)
                                                                      cTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTc                                                                        <  2:13298/147‑1 (MQ=255)
                                                                      cTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTc                                                                        >  1:13298/1‑147 (MQ=255)
                                                                      cTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGt                                                                      <  6:237582/149‑1 (MQ=255)
                                                                         aTCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTgg                                                                   <  1:207273/149‑1 (MQ=255)
                                                                            gCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTggg                                                                                                                          >  4:281242/1‑91 (MQ=255)
                                                                            gCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTggg                                                                                                                          <  3:281242/91‑1 (MQ=255)
                                                                                tCAGGCGCGACATCCACAAACATCCCGACACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAACCCTGGGTCGTTGGAACGCTCGTTGGACGCtgg                                                            >  8:86995/1‑147 (MQ=255)
                                                                                      gcgACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGa                                                      <  3:3935/149‑1 (MQ=255)
                                                                                              ccACAAACATCCCGAAACCGCTTTCGAAGAACCTCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCa                                               <  7:86995/149‑1 (MQ=255)
                                                                                              ccACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCAcc                                                >  3:20669/1‑148 (MQ=255)
                                                                                               cacaAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCAt                                              >  4:252954/1‑149 (MQ=255)
                                                                                                acaAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATc                                             <  6:277654/149‑1 (MQ=255)
                                                                                                    aCATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTaa                                         <  3:133141/149‑1 (MQ=255)
                                                                                                          cGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTcc                                   >  5:139949/1‑149 (MQ=255)
                                                                                                          cGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTcc                                   >  3:91316/1‑149 (MQ=255)
                                                                                                           gAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATGACCATCGTAAGCGTCCg                                  <  6:297022/149‑1 (MQ=255)
                                                                                                              aCCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCg                               <  7:284344/149‑1 (MQ=255)
                                                                                                              aCCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCg                               <  1:136263/149‑1 (MQ=255)
                                                                                                                cGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGaa                             <  3:252954/149‑1 (MQ=255)
                                                                                                                        aagaaCATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCtt                     >  4:342599/1‑149 (MQ=255)
                                                                                                                        aagaaCATCGCACCTCACAGCTGGTCGCAGACCGGCTCACTCAGACGGGCTTGGAGATTCACCGTCGACTGGCAAAAACCGGGGTAGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACAATCGTAAGCGTCCGGCGAGCTCACCtt                     >  7:289490/1‑149 (MQ=255)
                                                                                                                                 gCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTcc            <  4:91316/149‑1 (MQ=255)
                                                                                                                                 gCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTc             <  8:160365/148‑1 (MQ=255)
                                                                                                                                 gCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTc             <  1:132649/148‑1 (MQ=255)
                                                                                                                                  cACCTCACAGCTGGTCGCCGACCGGCTCACTCAGGGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCGCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCa           <  6:112045/149‑1 (MQ=255)
                                                                                                                                      agtgaGCCGGTCGGCGACCAGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCAt       >  7:359907/5‑149 (MQ=255)
                                                                                                                                         caGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGTGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTaa    <  4:20669/149‑1 (MQ=255)
                                                                                                                                          aGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGAGTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGc         >  1:147872/1‑143 (MQ=255)
                                                                                                                                          aGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGaa                >  7:336528/1‑136 (MQ=255)
                                                                                                                                          aGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGaa                <  8:336528/136‑1 (MQ=255)
                                                                                                                                                  gCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAgg  <  3:154843/142‑1 (MQ=255)
                                                                                                                                                  gCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAgg  >  4:154843/1‑142 (MQ=255)
                                                                                                                                                   ccGACCTGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGcc                                                          <  2:46523/85‑1 (MQ=255)
                                                                                                                                                   ccGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGcc                                                          >  1:46523/1‑85 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
ACCGCTAAATCGATCCAGGAACAAACTCATGAACCACGTGCTTGCCGACATTAAAGAGTTTCAACAAGAGCTCATCGCTATCAGGCGCGACATCCACAAACATCCCGAAACCGCTTTCGAAGAACATCGCACCTCACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGG  >  NC_002947/5032901‑5033187

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