Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,486,766 G→C 42.3% L180L (CTG→CTC yhdZ → amino acid ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,486,7660GC42.3% 686.2 / inf 537L180L (CTG→CTCyhdZamino acid ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (151/159);  new base C (106/121);  total (257/280)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

Only 100 of 544 total aligned reads displayed.

TGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCA  >  NC_002947/1486618‑1486915
                                                                                                                                                    |                                                                                                                                                     
   aCTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTACCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGAcc                                                                                                                                                   <  1:3749009/150‑1 (MQ=255)
           gAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGACATg                                                                                                                                           >  1:4787243/1‑150 (MQ=255)
            aGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATgg                                                                                                                                          <  1:7972358/150‑1 (MQ=255)
                 gtgaaaaTTCCGGAACAGGCCAACAAATATCCGGGCCAGCGTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAAGGGTCCAg                                                                                                                                     >  2:7481591/7‑150 (MQ=9)
                 gtgCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAg                                                                                                                                     <  2:7946708/150‑1 (MQ=255)
                  tgaaaaTTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAgg                                                                                                                                    <  1:6796648/145‑1 (MQ=9)
                   gaaaaTTCCGGAACAGGCCCACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGa                                                                                                                                   >  2:7288364/5‑150 (MQ=12)
                    cGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGaa                                                                                                                                  <  1:9015991/150‑1 (MQ=255)
                    cGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGGTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGaa                                                                                                                                  <  1:6418773/150‑1 (MQ=255)
                     ggATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAg                                                                                                                                 >  2:6949527/1‑150 (MQ=255)
                     ggATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAg                                                                                                                                 >  2:1233496/1‑150 (MQ=255)
                      aaTTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGt                                                                                                                                <  1:5852903/149‑1 (MQ=12)
                          ccGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCtg                                                                                                                            <  2:2255167/150‑2 (MQ=12)
                           cGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCtgg                                                                                                                           <  1:1658513/150‑3 (MQ=9)
                              aacAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATa                                                                                                                        <  2:3604595/148‑1 (MQ=9)
                                  ggCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCAt                                                                                                                    >  1:2305410/1‑150 (MQ=255)
                                     caacaAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGAGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGAtggt                                                                                                                 <  2:8096017/148‑1 (MQ=12)
                                      cacaAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCAtggtg                                                                                                                <  2:646792/150‑1 (MQ=255)
                                      aacaAATATCCGGGCCAGCTTTCGGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGAtggtg                                                                                                                <  2:274915/149‑1 (MQ=17)
                                      aacaAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGAtggtg                                                                                                                >  1:3546266/2‑150 (MQ=17)
                                       acaAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGat                                                                                                                                                    >  2:2964023/1‑112 (MQ=9)
                                       acaAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGAtggtga                                                                                                               <  2:4189255/150‑2 (MQ=17)
                                       acaAATATCCGGGCCAGCTTTCAAGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGAtggtga                                                                                                               <  2:9109883/150‑2 (MQ=12)
                                          aGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGATCAAGGAAGTGCTCGATACCATGGTGGGcc                                                                                                            >  2:7043000/1‑150 (MQ=255)
                                              cccTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGc                                                                                                        <  2:8910739/150‑1 (MQ=255)
                                               ccTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGATCCGACCTCCGCACTGGACCCGTAAATGGTCAAGGAAGTGATCGATACCATGGTGGGCCTGGcc                                                                                                       >  2:3246200/1‑150 (MQ=255)
                                               ccTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGcc                                                                                                       <  2:4679820/150‑1 (MQ=255)
                                               ccTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGAGACCATGGTGGGCCTGGcc                                                                                                       >  2:7709294/1‑150 (MQ=255)
                                                cTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGAGGGCCTGGCCg                                                                                                      >  1:8471157/1‑150 (MQ=255)
                                                 tGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGa                                                                                                     <  1:4731486/150‑1 (MQ=255)
                                                  ggCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAg                                                                                                    <  1:6279306/150‑1 (MQ=255)
                                                   gCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAgg                                                                                                   <  1:7573138/150‑1 (MQ=255)
                                                      aGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATg                                                                                                <  1:4504999/150‑1 (MQ=255)
                                                      aGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATg                                                                                                <  1:6497027/150‑1 (MQ=255)
                                                         tttcaggcggcCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCa                                                                                             <  1:1891995/145‑1 (MQ=17)
                                                         tttCTGGCGGGCAGCAGCAGCGTGTAGCCATCGCGCGTTCACTCTGCGTTAAACCATCGGTCCTGCTATGCGACGAGCCAACCTCCGCACTCGACCCGGAAATGGTCAAGGAGGTGCTGGAAACGTTGGCAGACCTGGCGCAGGACGGCa                                                                                             >  1:5678537/3‑150 (MQ=255)
                                                             aggcggcCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGAtggtga                                                                                                               <  2:6292527/127‑2 (MQ=17)
                                                              ggcgggCAGCAGCAGCGTGTCGCCATAGCGCGTTCCTTGTGCATGCAGCCCGACATCATCCTGTTCGACGAGCCGACCTCAGCCCTCGACCCTGAACTGGTGCGCGAAGTTCTCGATACGATGGTGTCGCTCGCCGAGGATGGTATGAcg                                                                                        >  2:1604665/1‑149 (MQ=255)
                                                                  gcCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATga                                                                                    <  1:6135247/150‑2 (MQ=17)
                                                                   ccAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATgat                                                                                   <  2:8778056/150‑3 (MQ=17)
                                                                   ccAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGGTGGCATGACCATGct                                                                                   >  1:3305996/1‑150 (MQ=255)
                                                                    cagcagcagCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGatc                                                                                  >  1:2848954/1‑150 (MQ=17)
                                                                    cagcagcagCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGctc                                                                                  >  2:4329836/1‑150 (MQ=255)
                                                                       cagcagcGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTgt                                                                               <  2:561609/150‑2 (MQ=17)
                                                                       cagcagcGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGc                                                                               <  2:5179440/150‑1 (MQ=255)
                                                                       cagcagcGTGTCGCCATAGCGCGTTCCTTGTGCATGCAGCCCGACATCATCCTGTTCGACGAGCCGACCTCAGCCCTCGACCCTGAACTGGTGCGCGAAGTTCTCGATACGATGGTGTCGCTCGCCGAGGATGGTATGACGATGATCTgt                                                                               <  2:4568935/150‑2 (MQ=255)
                                                                        agcagcGTGTCGCCATAGCGCGTTCCTTGTGCATGCAGCCCGACATCATCCTGTTCGACGAGCCGACCTCAGCCCTCGACCCTGAACTGGTGCGCGAAGTTCTCGATACGATGGTGTCGCTCGCCGAGGATGGTATGACGATGATCTgtg                                                                              <  2:2886123/150‑3 (MQ=255)
                                                                         gcagcGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTgtgt                                                                             <  1:4600628/150‑1 (MQ=17)
                                                                           agcGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGa                                                                           <  1:1182067/150‑1 (MQ=255)
                                                                                gtgGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCAc                                                                      >  1:3990997/1‑150 (MQ=12)
                                                                                gtgGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCAc                                                                      <  2:8735234/150‑1 (MQ=255)
                                                                                  ggCGATTGCCCGTGCGCTTTTCATGAACCCGAAGGTCATGATGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATATGTGTGACGCACGa                                                                    >  1:528935/1‑150 (MQ=12)
                                                                                  ggCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGa                                                                    >  1:8555904/1‑150 (MQ=255)
                                                                                   gcgaTTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGaa                                                                   >  1:7316208/4‑150 (MQ=12)
                                                                                   gCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGaa                                                                   >  1:2780643/1‑150 (MQ=255)
                                                                                   gCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGaa                                                                   >  1:8283073/1‑150 (MQ=255)
                                                                                        tGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATggg                                                              <  1:3431596/150‑1 (MQ=17)
                                                                                        tGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATggg                                                              <  2:6223159/150‑1 (MQ=255)
                                                                                        tGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATggg                                                              <  1:622967/150‑1 (MQ=255)
                                                                                          cccGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCt                                                            <  1:5452098/150‑1 (MQ=255)
                                                                                          cccGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCt                                                            <  2:6552710/150‑1 (MQ=255)
                                                                                           gcgttccttGTGCATGCAGCCCGACATCATCCTGTTCGACGAGCCGACCTCAGCCCTCGACCCTGAACTGGTGCGCGAAGTTCTCGATACGATGGTGTCGCTCGCCGAGGATGGTATGACGATGATCTGTGTAACGCACGAGATGGGttt                                                           <  2:211952/143‑1 (MQ=255)
                                                                                           ccGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCtt                                                           >  2:3950614/1‑150 (MQ=17)
                                                                                           ccGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCtt                                                           >  2:3151367/1‑150 (MQ=255)
                                                                                              ttccttGTGCATGCAGCCCGACATCATCCTGTTCGACGAGCCGACCTCAGCCCTCGACCCTGAACTGGTGCGCGAAGTTCTCGATACGATGGTGTCGCTCGCCGAGGATGGTATGACGATGATCTGTGTAACGCACGAGATGGGTTTTGc                                                        >  1:1001552/5‑150 (MQ=255)
                                                                                              cgcgTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCAGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGc                                                        <  2:5882555/150‑1 (MQ=255)
                                                                                               gcgTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGcc                                                       <  1:5063199/150‑1 (MQ=255)
                                                                                                cgTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGccc                                                      <  1:9120536/150‑1 (MQ=255)
                                                                                                  aacagCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGt                                                    >  2:5382258/5‑150 (MQ=17)
                                                                                                   tttgCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGtc                                                   >  1:5939349/3‑149 (MQ=17)
                                                                                                    ttgCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGtca                                                  <  1:7373054/149‑3 (MQ=17)
                                                                                                    gtgCATGCAGCCCGACATCATCCTGTTCGACGAGCCGACCTCAGCCCTCGACCCTGAACTGGTGCGCGAAGTTCTCGATACGATGGTGTCGCTCGCCGAGGATGGTATGACGATGATCTGTGTAACGCACGAGATGGGTTTTGCCCGtca                                                  <  1:319475/150‑3 (MQ=12)
                                                                                                       cATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGtcaggt                                               <  2:1792560/150‑6 (MQ=17)
                                                                                                       cATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGtcaggt                                               >  2:5032360/1‑145 (MQ=17)
                                                                                                        aTGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTg                                              >  1:2517796/1‑150 (MQ=255)
                                                                                                        aTGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTg                                              >  1:1543660/1‑150 (MQ=255)
                                                                                                        aTGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTg                                              <  2:3811539/150‑1 (MQ=255)
                                                                                                               ccAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAAcc                                       >  1:2209820/1‑150 (MQ=255)
                                                                                                                 aaGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCgg                                     <  2:3329238/150‑1 (MQ=255)
                                                                                                                  aGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCggg                                    <  1:4329836/150‑1 (MQ=255)
                                                                                                                        tGCTATGCGACGAGCCAACCTCCGCACTCGACCCGGAAATGGTCAAGGAGGTGCTGGAAACGTTGGCAGACCTGGCGCAGGACGGCAtca                                                                                          >  1:1975046/1‑88 (MQ=255)
                                                                                                                            gTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCAt                          <  2:1322233/150‑1 (MQ=12)
                                                                                                                            gTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCAt                          >  1:6649297/1‑150 (MQ=12)
                                                                                                                            gTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCAt                          <  1:1259590/150‑1 (MQ=12)
                                                                                                                            gTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCAt                          <  1:5612999/150‑1 (MQ=255)
                                                                                                                             ttCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCAt                          >  2:496942/1‑149 (MQ=255)
                                                                                                                             ttCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATg                         <  1:6700837/150‑1 (MQ=255)
                                                                                                                              tCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATgg                        >  1:7185233/1‑150 (MQ=12)
                                                                                                                              tCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATgg                        <  2:6755448/150‑1 (MQ=255)
                                                                                                                                   gaGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAg                   >  2:858599/1‑150 (MQ=255)
                                                                                                                                    aGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGatcagg                  <  1:3839739/150‑6 (MQ=255)
                                                                                                                                       cGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGcc               >  2:6097367/1‑149 (MQ=255)
                                                                                                                                        gACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGcca              <  1:5764292/150‑3 (MQ=255)
                                                                                                                                        gACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCga              <  2:2520975/150‑1 (MQ=255)
                                                                                                                                           cTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGAt           >  2:8331019/1‑150 (MQ=255)
                                                                                                                                            tCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGAtt          >  1:7711202/1‑150 (MQ=255)
                                                                                                                                             cGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGCCAGAtcg         <  1:6950756/150‑3 (MQ=255)
                                                                                                                                               gCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATtgtg       <  2:7556680/150‑1 (MQ=255)
                                                                                                                                                    ggACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCa  <  2:3756823/150‑1 (MQ=255)
                                                                                                                                                    ggACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCa  <  2:1003370/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCA  >  NC_002947/1486618‑1486915

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