Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,632,499 A→C 5.5% intergenic (+104/+55) PP_4096 → / ← PP_t59 hypothetical protein/tRNA‑Gly

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,632,4990AC 5.5% 233.0 / 2.3 72intergenic (+104/+55)PP_4096/PP_t59hypothetical protein/tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base A (35/33);  new base C (2/2);  total (37/35)
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.93e-01

TTGAAATACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGC  >  NC_002947/4632355‑4632643
                                                                                                                                                |                                                                                                                                                 
ttGAAATACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAggg                                                                                                                                               >  2:41024/1‑149 (MQ=255)
 ttaaaTACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGaa                                                                                                                                                  <  6:348213/143‑1 (MQ=255)
       aCCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCTTCGGGGCCtttt                                                                                                                                         <  4:309781/148‑1 (MQ=255)
       aCCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCTTCGGGGCCtttt                                                                                                                                         >  3:309781/1‑148 (MQ=255)
        cccGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCttttt                                                                                                                                        >  3:271875/1‑148 (MQ=255)
         ccGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCtttttt                                                                                                                                       >  7:172458/1‑148 (MQ=255)
         ccGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTc                                                                                                                                      <  2:376921/149‑1 (MQ=255)
                    aTCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTa                                                                                                                             >  7:113038/1‑147 (MQ=255)
                    aTCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACa                                                                                                                           >  1:136230/1‑149 (MQ=255)
                    aTCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACa                                                                                                                           <  8:113038/149‑1 (MQ=255)
                       aGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGa                                                                                                                         <  1:242732/148‑1 (MQ=255)
                       aGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGa                                                                                                                         <  8:178267/148‑1 (MQ=255)
                          gATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTc                                                                                                                     <  2:331072/149‑1 (MQ=255)
                              tGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGAc                                                                                                                  <  3:307747/148‑1 (MQ=255)
                               gTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACt                                                                                                                 >  2:121966/1‑148 (MQ=255)
                                  cAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGaaa                                                                                                             >  1:358750/1‑149 (MQ=255)
                                   atgTATCCCACGGCGTCACGGCCAGGGGGCCGGGTTGTGGTGTTTGTACTTTAATTTCGTAGAGGGCGGTGTATTGAGGCAGAATGCATTCTAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAAc                                                                                                            <  8:174037/147‑1 (MQ=255)
                                                cGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTTAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTa                                                                                               <  3:116244/149‑1 (MQ=255)
                                                  tCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTAtt                                                                                             <  1:41024/149‑1 (MQ=255)
                                                    aCGGCCAGGGCGCCGCGTTGTTGTGTTTGTACTTGAATTTCGTAGAGGTCGGGGTATTGAGGCAGAATTCATTCGAACGAAAAAAGGCCCCGAAGTGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTg                                                                                           <  3:271204/149‑1 (MQ=255)
                                                             gcgcCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGa                                                                                   >  7:112983/1‑148 (MQ=255)
                                                                ccgcGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACg                                                                               <  2:358750/149‑1 (MQ=255)
                                                                                   cTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGAccc                                                            >  2:121057/1‑149 (MQ=255)
                                                                                    ttGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGAcccc                                                           >  4:360703/1‑149 (MQ=255)
                                                                                    ttGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGAcccc                                                           >  1:88399/1‑149 (MQ=255)
                                                                                     tGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAACCCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGCACTCGCGAcccc                                                           >  2:249692/1‑148 (MQ=255)
                                                                                      gAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTc                                                                         >  2:377860/1‑133 (MQ=255)
                                                                                      gAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTc                                                                         <  1:377860/133‑1 (MQ=255)
                                                                                      gAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGAcccc                                                           <  2:136230/147‑1 (MQ=255)
                                                                                       aaTTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCg                                                                                      <  3:27721/119‑1 (MQ=255)
                                                                                       aaTTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCg                                                                                      >  4:27721/1‑119 (MQ=255)
                                                                                         tttCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCt                                                      >  1:36885/1‑149 (MQ=255)
                                                                                             gTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTa                                                                                                                             >  5:70454/1‑74 (MQ=255)
                                                                                             gTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTa                                                                                                                             <  6:70454/74‑1 (MQ=255)
                                                                                             gTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTgg                                                   >  2:46177/1‑148 (MQ=255)
                                                                                             gTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTgg                                                   <  8:172458/148‑1 (MQ=255)
                                                                                              tAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGc                                                  >  8:108492/1‑148 (MQ=255)
                                                                                                   ggCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTc                                           <  2:88399/149‑1 (MQ=255)
                                                                                                   ggCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTc                                           <  7:108492/149‑1 (MQ=255)
                                                                                                     cGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGt                                         >  7:320529/1‑149 (MQ=255)
                                                                                                      ggggTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTg                                        <  8:112983/149‑1 (MQ=255)
                                                                                                      ggggTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTg                                        <  4:271875/149‑1 (MQ=255)
                                                                                                        gctATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCggggt                                      <  1:121057/147‑5 (MQ=255)
                                                                                                           aTTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGt                                                                                                        >  5:178972/1‑81 (MQ=255)
                                                                                                           aTTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGt                                                                                                        <  6:178972/81‑1 (MQ=255)
                                                                                                            ttGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTa                                   <  2:265066/148‑1 (MQ=255)
                                                                                                               aGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCaa                               >  4:290472/1‑149 (MQ=255)
                                                                                                                 gCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGctc                                     <  8:348358/141‑1 (MQ=255)
                                                                                                                 gCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGctc                                     >  7:348358/1‑141 (MQ=255)
                                                                                                                   aGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACt                             >  3:112453/1‑147 (MQ=255)
                                                                                                                    gAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAg                          >  2:187408/1‑149 (MQ=255)
                                                                                                                     aaTGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGAc                                                              <  1:95769/113‑1 (MQ=255)
                                                                                                                     aaTGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGAc                                                              >  2:95769/1‑113 (MQ=255)
                                                                                                                      aTGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTc                                                                         >  5:308323/1‑101 (MQ=255)
                                                                                                                      aTGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTc                                                                         <  6:308323/101‑1 (MQ=255)
                                                                                                                         cATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTAtt                     >  7:166433/1‑149 (MQ=255)
                                                                                                                          aTTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTc                    >  5:128132/1‑149 (MQ=255)
                                                                                                                          aTTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTc                    <  4:112453/149‑1 (MQ=255)
                                                                                                                             cgaacgaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTccc                  <  8:166433/148‑1 (MQ=255)
                                                                                                                              gaacgaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTgg                                                   <  6:287864/115‑1 (MQ=255)
                                                                                                                              gaacgaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTgg                                                   >  5:287864/1‑115 (MQ=255)
                                                                                                                               aacgaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCCACTGAAACCTGGAGATCTTTATTTGGAGCGCGAAACGAGACACGAACTCGCGACCACGACCTTGGCAAGGACGTGCTCTACCAACTGAGCTATTccccc                >  3:85933/1‑146 (MQ=21)
                                                                                                                                ttgaaAAAAGGCCCTTCGGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCcgcg               >  6:25749/3‑148 (MQ=255)
                                                                                                                                ttgaaAAAAGGCCCTTCGGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCcgcg               <  5:25749/146‑1 (MQ=255)
                                                                                                                                  gaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGcgtcg            >  1:172224/1‑149 (MQ=255)
                                                                                                                                   aaaaaaGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTTCGTGCTCTACCAACTGAGCTATTCCCGcgtcg            <  1:46177/149‑1 (MQ=255)
                                                                                                                                     aaaaGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTCCCAACTGAGCTATTCGCGCGGCGTGa         >  4:31225/1‑149 (MQ=255)
                                                                                                                                          gCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGt    <  2:36885/149‑1 (MQ=255)
                                                                                                                                           ccccGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTg                                        >  3:78607/1‑112 (MQ=255)
                                                                                                                                           ccccGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTg                                        <  4:78607/112‑1 (MQ=255)
                                                                                                                                            cccGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACgaga                                                                            <  7:33447/76‑1 (MQ=255)
                                                                                                                                            cccGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACgaga                                                                            >  8:33447/1‑76 (MQ=255)
                                                                                                                                            cccGAAGGGACTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTAACAACTGTGCTATTGCCGCGTCGTGATGGGTGc  >  8:236113/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TTGAAATACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGC  >  NC_002947/4632355‑4632643

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: