Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,900,119 C→T 17.4% 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,1190CT17.4% 213.5 / 35.9 98intergenic (+35/‑463)PP_4308/PP_4309AsnC family transcriptional regulator/NCS1 family transporter
Reads supporting (aligned to +/- strand):  ref base C (36/44);  new base T (8/9);  total (44/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

Only 100 of 101 total aligned reads displayed.

GCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAACTAATCA  >  NC_002947/4899972‑4900267
                                                                                                                                                    |                                                                                                                                                    
gcTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCcgc                                                                                                                                                    >  1:913084/1‑150 (MQ=255)
 cTCGATGGAGCAGCTGGACACCCTGCTCGACCGCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGAGGTTGCCGCCAATGGATTAATTGGCTCCTGTGCTGGCCTCTTCGCAGGTGTGGCcgcg                                                                                                                                                   >  2:549992/1‑150 (MQ=255)
  tCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGa                                                                                                                                                  >  2:1285777/1‑150 (MQ=255)
   cGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGCCCCCGGTCAAGCGGCGGTTTCCGCCATTGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCTTGCCCGCGaa                                                                                                                                                 <  1:549992/150‑1 (MQ=255)
   cGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGaa                                                                                                                                                 <  2:2003915/150‑1 (MQ=255)
     aTGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGTGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAaga                                                                                                                                               <  2:34418/150‑1 (MQ=255)
      tGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAATAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCGGGCACGCCTGCGAagag                                                                                                                                              <  1:1749689/150‑1 (MQ=255)
      tGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAagag                                                                                                                                              <  2:2577977/150‑1 (MQ=255)
      tGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTACTGGCCGATACCAATACCGCGATCATCAAGATGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAagag                                                                                                                                              <  1:1869123/150‑1 (MQ=255)
        gAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGATCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGc                                                                                                                                            <  1:1098142/150‑1 (MQ=255)
        gAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGc                                                                                                                                            <  2:1212064/150‑1 (MQ=255)
         agcagcTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGCCCCATGAGAGCCATAGCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGcc                                                                                                                                           >  2:1518756/1‑150 (MQ=255)
          gcagcTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCa                                                                                                                                          <  2:1227086/150‑1 (MQ=255)
                     cccTGCTCGACCGCAGGAATAGGCTGGACGAGACCAATAGCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGACGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGaaaa                                                                                                                               <  2:682701/150‑1 (MQ=255)
                      ccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCTCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAg                                                                                                                              >  1:1814113/1‑150 (MQ=255)
                      ccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAg                                                                                                                              >  2:2620095/1‑150 (MQ=255)
                       cTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGc                                                                                                                             >  1:2281900/1‑150 (MQ=255)
                        tGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCt                                                                                                                            >  1:1763790/1‑150 (MQ=255)
                          cTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGa                                                                                                                          <  2:2695494/150‑1 (MQ=255)
                          cTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGa                                                                                                                          <  1:2586574/150‑1 (MQ=255)
                                cTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTAttt                                                                                                                    >  2:390926/1‑150 (MQ=255)
                                cTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTAttt                                                                                                                    <  1:1285777/150‑1 (MQ=255)
                                 tCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTAtttt                                                                                                                   <  2:1142703/150‑1 (MQ=255)
                                  cATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTAttttt                                                                                                                  >  2:1836537/1‑150 (MQ=255)
                                      aaTGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTa                                                                                                              >  1:166392/1‑150 (MQ=255)
                                             tGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACaaaa                                                                                                       >  1:942993/1‑150 (MQ=255)
                                             tGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACaaaa                                                                                                       <  1:390926/150‑1 (MQ=255)
                                             tGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACaaaa                                                                                                       >  1:2926139/1‑150 (MQ=255)
                                                  gagaCCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCCCTTCGCGGGCACGCCTGCGAagag                                                                                                                                              <  1:368058/106‑1 (MQ=38)
                                                  gagaCCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCCCTTCGCGGGCACGCCTGCGAagag                                                                                                                                              >  2:368058/1‑106 (MQ=38)
                                                       cAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGaa                                                                                                                                                 <  1:502616/98‑1 (MQ=255)
                                                       cAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGaa                                                                                                                                                 >  2:502616/1‑98 (MQ=255)
                                                       cAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAaga                                                                                                                                               >  2:1516428/1‑100 (MQ=255)
                                                       cAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAaga                                                                                                                                               <  1:1516428/100‑1 (MQ=255)
                                                        aaTACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTACTGGCCTCTTCGTGGGCACGCCTGCGAAGAGGCCAGCACAGGa                                                                                                                                  <  2:1320735/112‑1 (MQ=21)
                                                        aaTACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTACTGGCCTCTTCGTGGGCACGCCTGCGAAGAGGCCAGCACAGGa                                                                                                                                  >  1:1320735/1‑112 (MQ=21)
                                                        aaTACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTACTGGCCTCTTCGCGGGCACGCCTGCGAagag                                                                                                                                              >  1:1886502/1‑100 (MQ=25)
                                                        aaTACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTACTGGCCTCTTCGCGGGCACGCCTGCGAagag                                                                                                                                              <  2:1886502/100‑1 (MQ=25)
                                                         aTACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAg                                                                                           <  1:1836537/150‑1 (MQ=255)
                                                          tACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAgg                                                                                          >  2:2263365/1‑150 (MQ=255)
                                                               gcgATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCGGGCACGCCTGCGAaga                                                                                                                                               <  2:1280607/93‑1 (MQ=37)
                                                               gcgATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCGGGCACGCCTGCGAaga                                                                                                                                               >  1:1280607/1‑93 (MQ=37)
                                                                     atcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCAttttt                                                                                <  2:1814113/150‑1 (MQ=255)
                                                                      tcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCACGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCACAGGTGCGATTTTtg                                                                               >  2:1158185/1‑150 (MQ=255)
                                                                       caAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCt                                                                                                                            >  2:3047598/1‑104 (MQ=255)
                                                                       caAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCt                                                                                                                            <  1:3047598/104‑1 (MQ=255)
                                                                       caAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtgt                                                                              <  1:1518756/150‑1 (MQ=255)
                                                                         agaagaCCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGACAGCGAAGAGGCCAGTATAGGAAAAGCTGAGTATTAAAAAAAAAAAAAAACGTAACCATAGGTGCCATTTTtgtgt                                                                            <  2:1435856/150‑1 (MQ=255)
                                                                            agaCCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCa                                                                         <  1:2964027/150‑1 (MQ=255)
                                                                              aCCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACt                                                                       >  1:1759608/1‑150 (MQ=255)
                                                                              aCCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACt                                                                       <  2:2165494/150‑1 (MQ=255)
                                                                                cccGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtgt                                                                              <  2:493815/141‑1 (MQ=255)
                                                                                cccGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtgt                                                                              >  1:493815/1‑141 (MQ=255)
                                                                                    gTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTttcttc                                                                 >  1:354397/1‑150 (MQ=255)
                                                                                       aaGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACa                                                                                                                                     <  1:1843931/79‑1 (MQ=255)
                                                                                       aaGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACa                                                                                                                                     >  2:1843931/1‑79 (MQ=255)
                                                                                       aaGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGt                                                              <  1:2620095/150‑1 (MQ=255)
                                                                                         gcggcgGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCGGGCACGCCTGCGAAGAGGc                                                                                                                                            >  1:2490534/1‑70 (MQ=18)
                                                                                         gcggcgGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCGGGCACGCCTGCGAAGAGGc                                                                                                                                            <  2:2490534/70‑1 (MQ=21)
                                                                                         gcggcgGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGa                                                                                                                          >  1:1809804/1‑88 (MQ=255)
                                                                                         gcggcgGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGa                                                                                                                          <  2:1809804/88‑1 (MQ=255)
                                                                                              ggTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGccc                                                       <  2:1207890/150‑1 (MQ=255)
                                                                                                  gccgccAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAAt                                                   <  2:1763790/150‑1 (MQ=255)
                                                                                                   ccgccAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAAtt                                                  <  2:2453909/150‑1 (MQ=255)
                                                                                                     cccAATGGATTGATTGGCTCCTGTCCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAATCGTACCCAGAGGTGCCATTTTTGTGGGCACTTTCTTCAGTGGCGCCCCAATTCg                                                <  1:1158185/149‑1 (MQ=39)
                                                                                                        aaTGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAACAACCTGGCCAGAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTa                                             <  2:2281900/150‑1 (MQ=255)
                                                                                                        aaTGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTa                                             <  2:2926139/150‑1 (MQ=255)
                                                                                                        aaTGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTa                                             <  1:2263365/150‑1 (MQ=255)
                                                                                                           ggATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGGGGAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTacac                                          >  2:2124041/1‑150 (MQ=255)
                                                                                                                 attgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGccc                                                       <  1:2612644/131‑1 (MQ=255)
                                                                                                                 attgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGccc                                                       >  2:2612644/1‑131 (MQ=255)
                                                                                                                     gCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCAcc                                <  2:354397/150‑1 (MQ=255)
                                                                                                                     gCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCAcc                                >  1:828975/1‑150 (MQ=255)
                                                                                                                     gCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCAcc                                >  2:2780858/1‑150 (MQ=255)
                                                                                                                       tCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTTTTCTTTTTCTTCAGTTGTGCCCCAATTCGTTACACCCTTTCCACCct                              >  2:390778/1‑150 (MQ=255)
                                                                                                                         cTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCt                                                                                                                            <  1:3045215/54‑1 (MQ=255)
                                                                                                                         cTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCt                                                                                                                            >  2:3045215/1‑54 (MQ=255)
                                                                                                                         cTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCctct                            >  1:248783/1‑150 (MQ=255)
                                                                                                                          tgtgCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCAttt                                                                                  >  1:115171/1‑95 (MQ=255)
                                                                                                                          tgtgCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCAttt                                                                                  <  2:115171/95‑1 (MQ=255)
                                                                                                                          tgtgCTGGCCTCTTCGCAAGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTa                           >  2:1722442/1‑150 (MQ=255)
                                                                                                                            tgCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTttc                                                                    >  1:1599647/1‑107 (MQ=255)
                                                                                                                            tgCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTttc                                                                    <  2:1599647/107‑1 (MQ=255)
                                                                                                                               tGGCCTCTTCGCGGGCACGCCTGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtgtgt                                                                            <  1:312160/96‑1 (MQ=37)
                                                                                                                               tGGCCTCTTCGCGGGCACGCCTGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtgtgt                                                                            >  2:312160/1‑96 (MQ=25)
                                                                                                                                ggCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGa                     >  1:1855127/1‑150 (MQ=255)
                                                                                                                                 gCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGAc                    <  2:913084/150‑1 (MQ=255)
                                                                                                                                  ccTCTTCGCGGGCACGCCTGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCa                                                                <  1:1024222/105‑1 (MQ=37)
                                                                                                                                  ccTCTTCGCGGGCACGCCTGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCa                                                                >  2:1024222/1‑105 (MQ=37)
                                                                                                                                   ctctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAAc                                                                                                     <  1:669808/67‑1 (MQ=255)
                                                                                                                                   ctctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAAc                                                                                                     >  2:669808/1‑67 (MQ=255)
                                                                                                                                    tctTCGCGGGCACGCCTGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGt                                                                                         <  2:922656/78‑1 (MQ=25)
                                                                                                                                    tctTCGCGGGCACGCCTGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGt                                                                                         >  1:922656/1‑78 (MQ=25)
                                                                                                                                      ttCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAAt               >  1:415364/1‑150 (MQ=255)
                                                                                                                                      ttCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCAATTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAAt               >  2:2717830/1‑150 (MQ=255)
                                                                                                                                          cAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATc           <  2:1513492/150‑1 (MQ=255)
                                                                                                                                          cAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTAACCAAAAATGTAACTTTAAGGTCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATc           <  1:1898271/150‑1 (MQ=255)
                                                                                                                                               gTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAACTa      >  2:252110/1‑150 (MQ=255)
                                                                                                                                                tGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAACTaa     <  1:2780858/150‑1 (MQ=255)
                                                                                                                                                   ccGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTTTTTTTTTACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAACTaatca  >  2:46780/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAACTAATCA  >  NC_002947/4899972‑4900267

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