Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,900,074 C→T 5.9% intergenic (+35/‑463) PP_4308 → / → PP_4309 AsnC family transcriptional regulator/NCS1 family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,900,0740CT 5.9% 251.2 / 3.7 84intergenic (+35/‑463)PP_4308/PP_4309AsnC family transcriptional regulator/NCS1 family transporter
Reads supporting (aligned to +/- strand):  ref base C (36/43);  new base T (2/3);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-01

GCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAACTAAT  >  NC_002947/4899927‑4900220
                                                                                                                                                      |                                                                                                                                                  
gcTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCcg                                                                                                                                                   >  6:285980/1‑149 (MQ=255)
 cTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCcgc                                                                                                                                                  >  7:281865/1‑149 (MQ=255)
  tCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCcgc                                                                                                                                                  >  2:350613/1‑148 (MQ=255)
    gATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGaa                                                                                                                                               <  6:249664/149‑1 (MQ=255)
    gATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTACTGGCCTCTTCGCGGGCACGCCTGCGaa                                                                                                                                               <  2:195195/149‑1 (MQ=39)
     aTGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAag                                                                                                                                              >  7:146102/1‑149 (MQ=255)
      tGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCGGGCACGCCTGCGAaga                                                                                                                                             <  3:96438/149‑1 (MQ=255)
      tGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAaga                                                                                                                                             <  5:340796/149‑1 (MQ=255)
        gAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAagag                                                                                                                                            <  6:180508/148‑1 (MQ=255)
                   cacCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGCACAGGa                                                                                                                                <  5:228243/149‑1 (MQ=255)
                      ccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGAAGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAgg                                                                                                                                 >  1:117820/1‑148 (MQ=255)
                      ccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGaa                                                                                                                               >  4:69384/1‑147 (MQ=255)
                      ccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGaa                                                                                                                               >  7:174705/1‑147 (MQ=255)
                      ccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GCTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAACAGGCCAGTACAGGaa                                                                                                                               >  6:158416/1‑147 (MQ=255)
                          cTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCTGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCt                                                                                                                          <  8:174705/148‑1 (MQ=255)
                          cTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCt                                                                                                                          <  4:115521/148‑1 (MQ=255)
                           tCGAGCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCt                                                                                                                          >  3:216857/1‑147 (MQ=255)
                             gACCTCATGAATGTGCTGTCCTAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTTATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGa                                                                                                                        <  8:98293/147‑1 (MQ=255)
                             gACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGa                                                                                                                        <  1:350613/147‑1 (MQ=255)
                             gACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAg                                                                                                                       >  2:171290/1‑148 (MQ=255)
                                 tCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTAttt                                                                                                                  >  2:326903/1‑149 (MQ=255)
                                 tCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTACTGGCCTCTTCGCGGGCACGCCTGCGAagag                                                                                                                                            >  1:70436/1‑123 (MQ=38)
                                 tCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTACTGGCCTCTTCGCGGGCACGCCTGCGAagag                                                                                                                                            <  2:70436/123‑1 (MQ=38)
                                  cATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCTGC‑‑‑GGTTGCAGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGAAtttt                                                                                                                 >  7:126906/1‑149 (MQ=255)
                                  cATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTAtttt                                                                                                                 >  4:211336/1‑149 (MQ=255)
                                   aTGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTAttttt                                                                                                                <  8:188466/149‑1 (MQ=255)
                                    tGAATGTGCTGGCAGAGACCATTACCGCGTTCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTTGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTAtttttt                                                                                                               <  7:72519/149‑1 (MQ=255)
                                     gAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTg                                                                                                              >  3:87752/1‑149 (MQ=255)
                                      aaTGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTg                                                                                                              <  7:244633/148‑1 (MQ=255)
                                               gCCGAGACCAAAACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACaaaa                                                                                                     >  8:328512/1‑148 (MQ=255)
                                                ccGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGTCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAAc                                                                                                   <  8:49697/149‑1 (MQ=255)
                                                   agaCCAAGACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTa                                                                                                <  7:328512/149‑1 (MQ=255)
                                                          tACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAg                                                                                         <  1:326903/149‑1 (MQ=255)
                                                          tACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAg                                                                                         <  2:117820/149‑1 (MQ=255)
                                                             cgcgATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTg                                                                                      >  8:34595/1‑149 (MQ=255)
                                                              gcgATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGt                                                                                                             <  5:76859/125‑1 (MQ=255)
                                                              gcgATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGt                                                                                                             >  6:76859/1‑125 (MQ=255)
                                                              gcgATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGc                                                                                     <  7:34595/149‑1 (MQ=255)
                                                                gATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGACTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCa                                                                                   >  1:135144/1‑149 (MQ=255)
                                                                     tcaAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCTAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCAttttt                                                                              <  7:62892/149‑1 (MQ=255)
                                                                     tcaAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAACGTAACCATAGGTGCCATTTTtg                                                                             >  6:204962/1‑149 (MQ=255)
                                                                       aagaagACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtgt                                                                            <  5:271118/149‑1 (MQ=255)
                                                                       aagaagACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtg                                                                             <  3:69384/148‑1 (MQ=255)
                                                                         gaagaCCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtgtgt                                                                          <  8:146102/149‑1 (MQ=255)
                                                                           agaCCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGc                                                                        <  4:77347/149‑1 (MQ=255)
                                                                            gaCCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCa                                                                       <  3:250808/149‑1 (MQ=255)
                                                                               cccGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGt                                                                                                             <  4:293915/108‑1 (MQ=255)
                                                                               cccGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGt                                                                                                             >  3:293915/1‑108 (MQ=255)
                                                                               cccGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACtt                                                                    >  4:231696/1‑149 (MQ=255)
                                                                                 cGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTttc                                                                  <  5:285980/149‑1 (MQ=255)
                                                                                  ggTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTttct                                                                 <  5:21608/149‑1 (MQ=255)
                                                                                   gTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCt                                                                                                                          >  3:178997/1‑91 (MQ=255)
                                                                                   gTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCt                                                                                                                          >  3:173160/1‑91 (MQ=255)
                                                                                   gTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCt                                                                                                                          <  4:173160/91‑1 (MQ=255)
                                                                                   gTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCt                                                                                                                          <  4:178997/91‑1 (MQ=255)
                                                                                   gTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTttctt                                                                <  2:135144/149‑1 (MQ=255)
                                                                                      aaGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCa                                                              <  2:274090/148‑1 (MQ=255)
                                                                                           gc‑‑‑ggTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGgcg                                                        >  1:260913/1‑149 (MQ=255)
                                                                                                     ccgccAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCaa                                                  <  8:8217/148‑1 (MQ=255)
                                                                                                      cgccAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAAtt                                                <  3:211336/149‑1 (MQ=255)
                                                                                                              gattgattGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTaca                                         >  4:146996/1‑148 (MQ=255)
                                                                                                                ttgattgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGaa                                                                                                                                               >  6:183261/1‑44 (MQ=255)
                                                                                                                ttgattgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGaa                                                                                                                                               <  5:183261/44‑1 (MQ=255)
                                                                                                                   attgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACttt                                   >  6:40512/1‑149 (MQ=255)
                                                                                                                       gCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCAc                               <  5:341953/149‑1 (MQ=255)
                                                                                                                       gCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCAc                               <  5:313065/149‑1 (MQ=255)
                                                                                                                       gCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCAc                               >  6:32194/1‑149 (MQ=255)
                                                                                                                         tCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCAccc                             >  3:23361/1‑149 (MQ=255)
                                                                                                                         tCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCAccc                             >  5:200558/1‑149 (MQ=255)
                                                                                                                         tCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGAAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCAccc                             >  1:279518/1‑149 (MQ=255)
                                                                                                                           cTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCctc                           >  1:183735/1‑149 (MQ=255)
                                                                                                                            tgtgCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCctc                           >  2:56429/1‑148 (MQ=255)
                                                                                                                                cTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTAttt                      <  5:32194/149‑1 (MQ=255)
                                                                                                                                 tGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACAATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTAtttt                     >  8:282339/1‑149 (MQ=255)
                                                                                                                                   gCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTAtttt                     <  8:281865/147‑1 (MQ=255)
                                                                                                                                   gCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGa                   <  5:204962/148‑1 (MQ=255)
                                                                                                                                      tctTCGCGGGCACGCCTGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGAcc                 >  1:363489/1‑148 (MQ=255)
                                                                                                                                      tctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGAcc                 >  7:253493/1‑148 (MQ=255)
                                                                                                                                      tctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGAcc                 >  8:70714/1‑148 (MQ=255)
                                                                                                                                            cAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGAt          <  4:216857/149‑1 (MQ=255)
                                                                                                                                             aGGCGTGCCCGCTAAGAGGCCAGTACATGAAAAGCTGAGTATTTTTTGTACAAAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGAtg         <  5:158416/149‑2 (MQ=255)
                                                                                                                                             aGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATc         <  1:56429/149‑1 (MQ=255)
                                                                                                                                              ggCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCCGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATc         >  1:217566/1‑148 (MQ=255)
                                                                                                                                                cGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAAc      <  1:171290/149‑1 (MQ=255)
                                                                                                                                                 gTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTAtttt                     >  5:279260/1‑133 (MQ=255)
                                                                                                                                                 gTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTAtttt                     <  6:279260/133‑1 (MQ=255)
                                                                                                                                                 gTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAAc      >  2:3779/1‑148 (MQ=255)
                                                                                                                                                    cccGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACGCTCTATTTTGACCAAATGATCAACTaat  <  1:212626/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGC‑‑‑GGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAACTAAT  >  NC_002947/4899927‑4900220

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: