Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,835,646 T→C 7.1% intergenic (‑88/+44) PP_1639 ← / ← yohC SprT family protein/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,835,6460TC 7.1% 203.9 / 3.6 70intergenic (‑88/+44)PP_1639/yohCSprT family protein/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (33/32);  new base C (2/3);  total (35/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.65e-01

GAAAAAGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAG  >  NC_002947/1835499‑1835792
                                                                                                                                                   |                                                                                                                                                  
gAAATAGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGGCAGCGCACAGCCAGAGGCCGCCtt                                                                                                                                                    <  4:117852/148‑1 (MQ=255)
gAAAAAGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCtt                                                                                                                                                    >  1:242151/1‑148 (MQ=255)
gAAAAAGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCtt                                                                                                                                                    <  2:242151/148‑1 (MQ=255)
  aaaaGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTgcg                                                                                                                                                 >  4:293599/1‑149 (MQ=255)
  aaaaGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTgcg                                                                                                                                                 >  2:57804/1‑149 (MQ=255)
              cTTGCTGGTAACAGGTTTCGACGCGTTGTGTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGt                                                                                                                                     <  3:224578/149‑1 (MQ=255)
               ttGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGtt                                                                                                                                    >  1:53311/1‑149 (MQ=255)
                 gCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGGTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTg                                                                                                                                   <  7:185594/148‑1 (MQ=255)
                  cTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGGGGCCGCCGCAAGGCGGCCTCTGGCTgt                                                                                                                                  <  7:200463/148‑2 (MQ=255)
                   tGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATg                                                                                                                                >  1:243169/1‑149 (MQ=255)
                      tAACCGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGccc                                                                                                                             <  6:216495/149‑1 (MQ=255)
                        aCAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCa                                                                                                                           <  5:60919/149‑1 (MQ=255)
                        aCAGGGTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCa                                                                                                                           <  7:297104/149‑1 (MQ=255)
                               tCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTg                                                                                                                     <  5:330425/148‑1 (MQ=255)
                               tCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTg                                                                                                                     <  2:184530/148‑1 (MQ=255)
                                         gttTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTg                                                                                                           <  7:346099/148‑1 (MQ=255)
                                            tGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTGCGGTCTGTTc                                                                                                        >  6:177682/1‑148 (MQ=255)
                                              aGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAg                                                                                                     <  5:210168/149‑1 (MQ=255)
                                               gctgctCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGt                                                                                                    >  5:107225/1‑149 (MQ=255)
                                               gctgctCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGt                                                                                                    <  6:244246/149‑1 (MQ=255)
                                               gctgctCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGt                                                                                                    >  5:244246/1‑149 (MQ=255)
                                               gctgctCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGt                                                                                                    >  8:107331/1‑149 (MQ=255)
                                                   ctCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGgt                                                                                                >  5:44695/1‑149 (MQ=255)
                                                        gCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGGGGCCGCCGCAAGGCGGCCTCTgg                                                                                                                                      <  1:359171/106‑1 (MQ=38)
                                                        gCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGGGGCCGCCGCAAGGCGGCCTCTgg                                                                                                                                      >  2:359171/1‑106 (MQ=38)
                                                        gCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGa                                                                                           <  5:24616/149‑1 (MQ=255)
                                                          aTGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGgt                                                                                                >  3:237529/1‑142 (MQ=255)
                                                          aTGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGgt                                                                                                <  4:237529/142‑1 (MQ=255)
                                                            ggcggggcggTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACggg                                                                                       >  5:306999/1‑149 (MQ=255)
                                                             gcggggcggTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGc                                                                                      <  3:293599/149‑1 (MQ=255)
                                                               ggggcggTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGccc                                                                                    >  6:57154/1‑149 (MQ=255)
                                                                        ctcCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGcc                                                                              >  8:260840/1‑146 (MQ=255)
                                                                        ctcCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCa                                                                            <  6:107225/148‑1 (MQ=255)
                                                                        ctcCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCa                                                                            <  2:53311/148‑1 (MQ=255)
                                                                          ccTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGACGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTc                                                                         >  8:160464/1‑149 (MQ=255)
                                                                           cTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTcc                                                                        >  8:321281/1‑149 (MQ=255)
                                                                           cTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTcc                                                                        >  3:133830/1‑149 (MQ=255)
                                                                                 aGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGa                                                                   >  4:296020/1‑148 (MQ=255)
                                                                                 aGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGa                                                                   >  1:239922/1‑148 (MQ=255)
                                                                                     aTTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTc                                                                         >  1:154426/1‑138 (MQ=255)
                                                                                     aTTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTc                                                                         <  2:154426/138‑1 (MQ=255)
                                                                                      ttATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCa                                                             >  4:199003/1‑149 (MQ=255)
                                                                                        aTGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCAc                                                            <  1:57804/148‑1 (MQ=255)
                                                                                         tGCCGCAAGCGCAGGTGGCTTGCCAAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACgg                                                          >  7:231107/1‑149 (MQ=255)
                                                                                           ccGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGTGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACgg                                                          >  5:280954/1‑147 (MQ=255)
                                                                                            cGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAg                                                       >  4:331733/1‑149 (MQ=255)
                                                                                            cGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAg                                                       <  5:57154/149‑1 (MQ=255)
                                                                                            cGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAg                                                       >  2:344873/1‑149 (MQ=255)
                                                                                                       gTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGAtt                                            <  7:107331/149‑1 (MQ=255)
                                                                                                         ggCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGc                                          <  7:260840/149‑1 (MQ=255)
                                                                                                         ggCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGc                                          <  8:188914/149‑1 (MQ=255)
                                                                                                                    gcACCGGTCAGCGCACAGCCAGATGCCGCCTTGGGGCGGCCCCTGGGTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGtacta                               <  5:177682/148‑1 (MQ=255)
                                                                                                                       ccGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTAcc                             >  6:137856/1‑148 (MQ=255)
                                                                                                                        cGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCaa                           <  7:160464/149‑1 (MQ=255)
                                                                                                                          gTCAGCGCACAGCCAGAGGCGGCATTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGc                         <  5:305794/149‑1 (MQ=255)
                                                                                                                           tCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGcc                        <  2:239922/149‑1 (MQ=255)
                                                                                                                           tCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGcc                        <  3:331733/149‑1 (MQ=255)
                                                                                                                               cgcACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACt                    >  6:217905/1‑149 (MQ=255)
                                                                                                                                 cacaGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGc                  >  1:231986/1‑149 (MQ=255)
                                                                                                                                  tcaaccAGGGGCCGCCGCAAGGCGGCCTCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGc                  <  7:109242/144‑1 (MQ=255)
                                                                                                                                  tcaaccAGGGGCCGCCGCAAGGCGGCCTCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGc                  >  8:109242/5‑148 (MQ=255)
                                                                                                                                         gagGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCAcc           >  6:295723/1‑148 (MQ=255)
                                                                                                                                         gagGCCGCCTTGCGGCGGCCCCTGGTAGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCCGCAcc           >  3:310430/1‑148 (MQ=255)
                                                                                                                                          agGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCAcc           >  6:64357/1‑147 (MQ=255)
                                                                                                                                           ggCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGc              >  1:21524/1‑143 (MQ=255)
                                                                                                                                           ggCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGc              <  2:21524/143‑1 (MQ=255)
                                                                                                                                           ggCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAg        <  5:64357/149‑1 (MQ=255)
                                                                                                                                            gccggcTTGCGGGGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGTTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAg        <  4:310430/148‑1 (MQ=255)
                                                                                                                                            gccgccTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGcc                                                  >  3:99085/1‑106 (MQ=255)
                                                                                                                                            gccgccTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGcc                                                  <  4:99085/106‑1 (MQ=255)
                                                                                                                                            gccgccTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGc       >  1:306271/1‑149 (MQ=255)
                                                                                                                                               gccTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTgg    <  3:199003/149‑1 (MQ=255)
                                                                                                                                                 cTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGgag  <  6:44695/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GAAAAAGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAG  >  NC_002947/1835499‑1835792

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: