Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 875,954 A→G 9.7% intergenic (‑63/+24) PP_0760 ← / ← PP_mr65 membrane protein/YybP‑YkoY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947875,9540AG 9.7% 168.1 / 7.6 62intergenic (‑63/+24)PP_0760/PP_mr65membrane protein/YybP‑YkoY
Reads supporting (aligned to +/- strand):  ref base A (28/28);  new base G (4/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CGAGAATGAGTGCGAGCAGTTGCGTCTTGTCGCCGATCTCGGCCAGGGCAACGATAGCGGTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCTGTGGACCAAGTATGTTGACGCGTACCGGGCGAGCTGGGTGCTCGCGGGAGA  >  NC_002947/875807‑876098
                                                                                                                                                   |                                                                                                                                                  
cGAGAATGAGTGCGAGCAGTTGCGTCTTGTCGCCGATCTCGGCCAGGGCAACGATAGCGGTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAg                                                                                                                                                   <  8:118533/149‑1 (MQ=255)
      tGAGTGCGAGCAGTTGCGTCTTGTCGCCGATCTCGGCCAGGGCAACGATAGCGGTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGc                                                                                                                                             <  3:995/149‑1 (MQ=255)
         gTGCGAGCAGTTGCGTCTTGTCGCCGATCTCGGCCAGGGCAACGATAGCGGTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCgtg                                                                                                                                          >  5:312059/1‑149 (MQ=255)
          ggCGAGCAGCTGCGTCTTGTCGCCGATCTCGGCCAGGGCAACGATAGCGGTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCgtgt                                                                                                                                         <  8:313052/148‑1 (MQ=255)
                  gTTGCGTCTTGTCGCCGATCTCGGCCAGGGCAACGATAGCGGTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTc                                                                                                                                 <  1:251445/149‑1 (MQ=255)
                   ttGCGTCGTGGCGCCGATCTCGGCCAGGGCAACGATAGCGGTGGGGACCACCAGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGCACACCTTACCCGGGGCGGGGAGGCTGTACGTGTCATAgccgt                                                                                                                                <  5:197058/149‑5 (MQ=21)
                   ttGCGTCATGTCGCCGATCTCGGCCAGGGCAACGATAGCGGTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGGGGGCCGGTCGAGAAGGCTATGACACGTACAGCCTCACCGCCCCGGGTAAGGTGTGCCTGTCAAAGGTCt                                                                                                                                >  1:277558/1‑149 (MQ=255)
                            gTCGCCGATCTCGGCCAGGGCAACGATAGCGGTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCaa                                                                                                                          <  4:197751/146‑1 (MQ=255)
                            gTCGCCGATCTCGGCCAGGGCAACGATAGCGGTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑cc                                                                                                                     >  5:415733/1‑149 (MQ=255)
                             tCGCCGATCTCGGCCAGGGCAACGATAGCGGTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAg                                                                                                                                                   <  2:391470/120‑1 (MQ=255)
                             tCGCCGATCTCGGCCAGGGCAACGATAGCGGTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAg                                                                                                                                                   >  1:391470/1‑120 (MQ=255)
                               gCCGATCTCGGCCAGGGCAACGATAGCGGTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAGGCTGTACGTGTCa                                                                                                                                       >  4:77900/1‑130 (MQ=255)
                               gCCGATCTCGGCCAGGGCAACGATAGCGGTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGCGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAGGCTGTACGTGTCa                                                                                                                                       <  3:77900/130‑1 (MQ=255)
                                ccGATCTCGGCCAGGGCAACGATAGCGGTGGTGACCAACAGTGATTCCAGCATCAGGATTCCTACGGTGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCTGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCg                                                                                                                 <  3:272661/149‑1 (MQ=255)
                                        ggCCAGGGCAACGATAGCGGTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATc                                                                                                              <  3:34782/144‑1 (MQ=255)
                                        ggCCAGGGCAACGATAGCGGTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGt                                                                                                           <  1:179862/147‑1 (MQ=255)
                                            aGGGCAACGATAGCGGTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGc                                                                                                     <  6:312059/149‑1 (MQ=255)
                                                cAACGATAGCGGTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCtt                                                                                                                               <  8:166958/121‑1 (MQ=255)
                                                cAACGATAGCGGTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCtt                                                                                                                               >  7:166958/1‑121 (MQ=255)
                                                           gTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCa                                                                                      <  8:192200/149‑1 (MQ=255)
                                                                       aGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCAtggt                                                                          >  4:150810/1‑149 (MQ=255)
                                                                         tGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGt                                                                                 >  8:236619/1‑140 (MQ=255)
                                                                         tGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGt                                                                                 <  7:236619/140‑1 (MQ=255)
                                                                              ccAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAATGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCACGCCTGCAtggagagatc                                                                    <  8:160958/148‑8 (MQ=255)
                                                                                 ggttCAGGATTCCTTCGGGGGCGGGTGGACACGGCTATGACACGTACTGCGTCGCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCGATGCAGGCGTGCATGGTGGGTCGCAcg                                                                <  2:277558/146‑1 (MQ=255)
                                                                                  cATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCAcgc                                                               >  7:375140/1‑149 (MQ=255)
                                                                                    tCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACgcgc                                                             >  5:315751/1‑149 (MQ=255)
                                                                                     cAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACgcgc                                                             >  1:365470/1‑148 (MQ=255)
                                                                                      aGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAg                                                                                                                                    <  8:335143/78‑1 (MQ=255)
                                                                                      aGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAg                                                                                                                                    >  7:335143/1‑78 (MQ=255)
                                                                                          ttCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCAcgcg                                                              >  6:19765/1‑142 (MQ=255)
                                                                                          ttCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCAcgcg                                                              <  5:19765/142‑1 (MQ=255)
                                                                                             cTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCtg                                                    >  7:128901/1‑149 (MQ=255)
                                                                                             cTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCtg                                                    >  1:311729/1‑149 (MQ=255)
                                                                                             cTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCtg                                                    >  1:38439/1‑149 (MQ=255)
                                                                                              taagggggCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCtgt                                                   >  5:65438/4‑149 (MQ=255)
                                                                                              tACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTTTGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCtgt                                                   >  1:331381/1‑149 (MQ=255)
                                                                                              tACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCtgt                                                   >  1:382397/1‑149 (MQ=255)
                                                                                              tACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCtgt                                                   >  3:46392/1‑149 (MQ=255)
                                                                                                cGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAATGGTCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCtgt                                                   <  1:126031/149‑1 (MQ=255)
                                                                                                cGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCtgtg                                                  >  6:369366/1‑148 (MQ=255)
                                                                                                     gCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCTGTGGACCaa                                            >  6:205411/1‑149 (MQ=255)
                                                                                                         gTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCtgt                                                   <  1:404107/138‑1 (MQ=255)
                                                                                                         gTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCtgt                                                   >  2:404107/1‑138 (MQ=255)
                                                                                                         gTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCTGTGGACCAAGTAt                                        <  6:315751/149‑1 (MQ=255)
                                                                                                          tCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCTGTGGACCAAGTATg                                       <  8:128901/149‑1 (MQ=255)
                                                                                                            gACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCTGTGGACCAAGTATGt                                      <  8:375140/148‑1 (MQ=255)
                                                                                                             acacGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCAcgc                                                               >  6:314883/1‑122 (MQ=255)
                                                                                                             acacGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCAcgc                                                               <  5:314883/122‑1 (MQ=255)
                                                                                                             acacGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCTGTGGACCAAGTATGtt                                     >  5:299086/1‑148 (MQ=255)
                                                                                                                   cTATGACACGCACACCTTACCCGGGGCGGGGAGGCTGTACGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCTGTGGACCAAGTATGTTGAcgcg                               >  4:220466/1‑148 (MQ=25)
                                                                                                                   cTATGACACGCACACCTTACCCGGGGCGGGGAGGCTGTACGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCTGTGGACCAAGTATGTTGAcgcg                               >  1:362745/1‑148 (MQ=25)
                                                                                                                    tATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGCTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTATCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCTGTGGAGCAAGTATGTTGACGCGTa                             >  5:252540/1‑149 (MQ=255)
                                                                                                                    tATGACACGCACACCTTACCCGGGGCGGGGAGGCTGTACGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCTGTGGACCAAGTATGTTGACGCGTa                             >  4:375061/1‑149 (MQ=25)
                                                                                                                     aTGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATACGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGTTCGCACGCGCCATGGTCTGTGGACCAAGTATGTTGACGCGTAc                            >  5:236596/1‑149 (MQ=255)
                                                                                                                       gACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCTGTGGACCAAGTATGTTGACGCGTACCg                          >  7:283271/1‑149 (MQ=255)
                                                                                                                          acGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCTGTGGACCAAGTATGTTGACGCGTACCGGGc                       <  6:65438/149‑1 (MQ=255)
                                                                                                                                 gCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCTGTGGACCAAGTATGTTGACGCGTACCGGGCGAGCTgg                <  2:365470/149‑1 (MQ=255)
                                                                                                                                        cgccccgGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCTGTGGACCAAGTATGTTGACGCGTACCGGGCGAGCTGGGTGCTcg         >  6:14026/1‑149 (MQ=255)
                                                                                                                                            ccgGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCTGTGGACCAAGTATGTTGACGCGTACCGGGCGAGCTGGGTGCTCGAggg     <  2:331381/149‑1 (MQ=255)
                                                                                                                                             cgGGTTAGGTGTGCGTGTGATAGGTCTTGTCATA‑‑CCCTCTATCCGTCTATGCGGATCTCTCCATGCAGGCGTGGATGGTGGGTCGCACGCGCCATGGTCTGTGGTCCAAGTATGTTGACGCGTACGGGGCGAGCTTGGTGCTGGCGGga    <  4:46392/149‑1 (MQ=255)
                                                                                                                                               ggTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCTGTGGACCAAGTATGTTGACGCGTACCGGGCGAGCTGGGTGCTCGCGGgaga  <  2:382397/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGAGAATGAGTGCGAGCAGTTGCGTCTTGTCGCCGATCTCGGCCAGGGCAACGATAGCGGTGGGGACCAACAGTGATTCCAGCATCAGGATTCCTACGGGGGCGGGTCGACACGGCTATGACACGTACAGCCTCCCCGCCCCGGGTAAGGTGTGCGTGTCATAGGTCTTGTCAAA‑‑CCCTCGATCCGTCTATGCGGATCTCTCCATGCAGGCGTGCATGGTGGGTCGCACGCGCCATGGTCTGTGGACCAAGTATGTTGACGCGTACCGGGCGAGCTGGGTGCTCGCGGGAGA  >  NC_002947/875807‑876098

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: