Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,580,393 T→A 6.8% intergenic (+105/‑219) PP_3158 → / → PP_5589 hypothetical protein/rhs repeat‑associated core domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,580,3930TA 6.8% 683.7 / 11.5 253intergenic (+105/‑219)PP_3158/PP_5589hypothetical protein/rhs repeat‑associated core domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (108/119);  new base A (21/2);  total (132/121)
Fisher's exact test for biased strand distribution p-value = 3.77e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 266 total aligned reads displayed.

GCGCAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGT  >  NC_002947/3580247‑3580542
                                                                                                                                                  |                                                                                                                                                     
gcgcAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAtttt                                                                                                                                                    >  1:3000654/1‑150 (MQ=255)
 cgcAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAttttt                                                                                                                                                   <  1:1444148/150‑1 (MQ=255)
 cgcAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTCTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAttttt                                                                                                                                                   <  2:1935144/150‑1 (MQ=255)
   cAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGc                                                                                                                                                 <  2:4106612/150‑1 (MQ=255)
     gagCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCa                                                                                                                                               <  1:7509275/150‑1 (MQ=255)
        ccTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGtt                                                                                                                                            <  1:3136153/150‑1 (MQ=255)
         cTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTg                                                                                                                                           <  2:449973/150‑1 (MQ=255)
          tGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTGCCACTTGc                                                                                                                                          >  2:554302/1‑150 (MQ=255)
          tGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAAAGGCATTTTTGCCAGTTGc                                                                                                                                          >  1:1338164/1‑150 (MQ=255)
           gTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCa                                                                                                                                         <  1:2125611/150‑1 (MQ=255)
            ttGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCGGCAACTGGCAGTTTTGGCAGTTGCAg                                                                                                                                        >  2:4637740/1‑150 (MQ=255)
               cTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGa                                                                                                                                     <  1:7486457/150‑1 (MQ=255)
                tgatgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAATGGCATTTTTTCCAGTTTCAGCGAg                                                                                                                                    >  2:5805032/1‑150 (MQ=255)
                 gatgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAATGGGATTTTTTCCAGTTGCATCGAgt                                                                                                                                   >  1:3576784/1‑150 (MQ=255)
                  atgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAgtg                                                                                                                                  >  2:4193370/1‑150 (MQ=255)
                   tgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAgtgt                                                                                                                                 >  2:320731/1‑150 (MQ=255)
                   tgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAAGGCACATACTTGAATTAGGATTGAGTTTATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGAAACTGGCATTTTTGCCAGTTGCAGCGAgtgt                                                                                                                                 <  2:7262907/150‑1 (MQ=255)
                    taaCTTCCTGCGGGTGCATTATGAGGCCTCTGCCGCACGCCCTTTGCTAAATCCCCATACTTGATTGAGGTTTGTTTTGATGTGTGTTCATGAGCTAATCAAGTTTTCGTCGGGCTGAAACTGGCTTTTTTGCCATTTGCAGCGAGTGTa                                                                                                                                <  1:2565737/149‑1 (MQ=255)
                       cTTCCTGCCGGTGCAGTAAGAGGCCTCGGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGa                                                                                                                             <  2:5914446/150‑1 (MQ=255)
                        ttCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGGGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTGCCAGTTTGAGCGAgtgtgtttc                                                                                                                            >  1:3955917/1‑145 (MQ=255)
                        ttCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTGCCAGTTGCAGCGAGTGTATTAc                                                                                                                            >  1:5932186/1‑150 (MQ=255)
                        ttCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCAATTTTTCCAGTTTCAGCGAGTGTATGAc                                                                                                                            >  1:2764402/1‑150 (MQ=255)
                            tGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGcc                                                                                                                        >  2:4967794/1‑150 (MQ=255)
                                     aGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCa                                                                                                               <  1:9122685/150‑1 (MQ=255)
                                     aGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTTCCAGGTGCAGCCAGTGGATTACAACCAAGaaaaca                                                                                                               >  1:1607789/1‑146 (MQ=255)
                                      ttAAGTGGCCTCTGCCGCACGCCCTTTGCGAAATGCCAATACATTAATGAGTTTTGAGTTTATGTGTTTTCATGCGCTAATCAAGGTTTCGTGGGGCTGAAACTGGCATTTTTGCCTTTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                              <  1:1171380/149‑1 (MQ=255)
                                      gTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACTTGAATGAGGATTGAGTTGATGAGTTTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                              <  2:3415570/150‑1 (MQ=255)
                                      gTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                              <  2:2400284/150‑1 (MQ=255)
                                      gTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTTCCAGTTTCAGCGAGGGGATGACAGCCAAGaaataaa                                                                                                              >  1:7759315/1‑145 (MQ=255)
                                      gTAAGAGGCCTCTGCCCCACGCCATTTGCTAAATGCCCATACATTAATGAGTATTGAGTTGATGAGTGTTCATGAGCCAAACAACGTTTCGTCGGGCTGAAACTGCAATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                              <  1:1032936/150‑1 (MQ=255)
                                        aaGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACTTGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGAATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAgg                                                                                                            <  2:1205552/150‑1 (MQ=255)
                                        aaGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCTTGAGCGAATCAAGGTTTCGTCGGGCTGCAACTGGAATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAgg                                                                                                            <  1:4545476/150‑1 (MQ=255)
                                        aaGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATAAATGAATGAGTATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTGGGGCTGCAACTGGAATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAgg                                                                                                            <  1:6242670/150‑1 (MQ=255)
                                          gagGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAATGGGAATTTTTGCACTTGCAGCGAGTTTTTTAAAACCAAGAATTAAAGGGa                                                                                                          >  1:6442551/1‑150 (MQ=255)
                                              ccTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATg                                                                                                      <  2:1754441/150‑1 (MQ=255)
                                              ccTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCCTTTTTGCCAGTTGCAGCGAGTGTATGTAAGCCGAGAAATGAAGGGAAATg                                                                                                      >  1:2014201/1‑150 (MQ=255)
                                                tctGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGCGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGgt                                                                                                    <  1:2657508/150‑1 (MQ=255)
                                                tctGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGgt                                                                                                    <  2:3704270/150‑1 (MQ=255)
                                                tctGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGgt                                                                                                    <  2:1414999/150‑1 (MQ=255)
                                                 ctGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGgtg                                                                                                   >  2:5315765/1‑150 (MQ=255)
                                                  tGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGGGGCAACAAGCAATTTTTCCACGTTCAGCCAGTGTATGACAACCAAGAAATAAAGGGAAAAGgtgt                                                                                                  >  1:2495053/1‑150 (MQ=255)
                                                   gCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGtt                                                                                                 <  2:4827190/150‑1 (MQ=255)
                                                    ccGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTAAAGGGAAAAGGTGTtc                                                                                                >  2:6083110/1‑150 (MQ=255)
                                                    ccGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTTCCACTTTCAGCGAGTGTATGAAAGCCAAGAAATAAAGGGAAAATGTGttg                                                                                                >  1:4228462/1‑149 (MQ=255)
                                                    ccGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTGCCAGTTGCAGCGAGTGTATGCCAGCCTAGAATTCAAGGGAAATGGGGTtc                                                                                                >  2:8305216/1‑150 (MQ=255)
                                                      gCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTTCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTtctc                                                                                              >  1:5339743/1‑150 (MQ=255)
                                                         cGCCCTTGGCGAAATGCCAATACATGAATGAGGATTGATTTGATGATTGTTCATGAGCTAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTaa                                                                                           <  1:8947793/150‑1 (MQ=255)
                                                          gCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTTCAAAAGGCAATTTTTCCAGTTGCAGCGAGTGTATTAAACCCAAGAAATAAAAGGAAAATGTGTTATATaat                                                                                          >  1:7312930/1‑149 (MQ=255)
                                                          gCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAATTTTGCCCGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGGGTTCGCGAAg                                                                                          >  2:7107621/1‑150 (MQ=255)
                                                            ccTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCGGCAACAGGCATTTTTTGCCGTTGCAGCGCGTGGATGGCCCCCCACAATTAAAGGGCAAAGGGGGTGTCTGAggg                                                                                        >  1:32976/1‑149 (MQ=255)
                                                              tttGCTAATTGCCCATACATGATTGAGTTTTTTGTTTATGAGTGTTCTTGTGCGAATAAACGTTTGGTGGGGCTGAACCTGGAATTTTTGCCTGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCa                                                                                      <  1:2516441/150‑1 (MQ=255)
                                                              tttGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCTAATCAACGTTGCGTCGGGCTGCAACTGGCTTTTTTGCCTGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCa                                                                                      <  1:4833422/150‑1 (MQ=255)
                                                                 gCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGGTGCAACAGGCCATTTTTCCACGTGCAGCGGGTTTTTGACAGCCAAGAAATGAAGGGAAAATGTGTTGTCGAAGGTCATTa                                                                                   >  1:2311774/1‑150 (MQ=255)
                                                                   gAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCAATTTTTCCAGTTTCAGCGAGTGTATTACAGCCTAGAAATCAAGGGAAATGGTGTTCTCTAAGGTCATTAAt                                                                                 >  2:2468450/1‑150 (MQ=255)
                                                                    aaaTGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCAGCAACACGCAGATTTGCCACTTGCAGCGCGTGTAAGACCGCCCACAATACAAGGGAAAAGGGGGTCACAAAGGTCAATaaac                                                                                >  2:5673081/1‑148 (MQ=255)
                                                                      aTGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCAAATTTGCCAGTTTAACCCAGTTTATTAAACACAATAATTAAATGTAAAAGGTGTTATATAAGGTAAATaaaagt                                                                              >  2:2742809/1‑146 (MQ=255)
                                                                        gCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGa                                                                            <  2:5635077/150‑1 (MQ=255)
                                                                        gCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGAAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGa                                                                            <  1:8104987/150‑1 (MQ=255)
                                                                         cccATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGaa                                                                           <  2:4767921/150‑1 (MQ=255)
                                                                               cATGAATGAGGTTTGTGTGGATGAGTGTTCATGAGCTAATCAACGTTTCGTCGGGCTGCAACTGGAATTTTTGCCTGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTc                                                                     <  1:148481/150‑1 (MQ=255)
                                                                               cATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTTCAACAGGCATTTTTGCCAGTTGCAGCGAGTGTATGAAAGCCTAGAATTCAAGGGAAAATGTGTTATCTAAAGTAATTAATCCTGAAACGTTc                                                                     >  2:5226651/1‑150 (MQ=255)
                                                                               cATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTc                                                                     <  1:672280/150‑1 (MQ=255)
                                                                                   aatgaGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCAATTTTTCCACTTGCAGCGAGTGTATGACACCCAAGAAATAAAGGGAAAATGTGTTATCAAAGGTCAATAATCCTGAAAACATCGCaa                                                                 >  1:3516440/1‑150 (MQ=255)
                                                                                     tgaGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTTCCAGTTGCAGCGAGTGTAAGACAGCCTAGAATTCAAGGGACATGGTGTTGTCTAAGGTCAATAATCGTGAAAACTTAGCAAGc                                                               >  2:2807378/1‑150 (MQ=255)
                                                                                      gaGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTGCCAGTTGCAGCGAGTGTATTACAGCCTAGAAATCAACGGAAATGGTGTTCTCTAAGGTCATTAATAGTGAAACGTTAGCAAGCg                                                              >  2:1608810/1‑150 (MQ=255)
                                                                                        ggATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCAAATTTTCCACTTGCAGCGAGAGTAAGACACCCAAGAAAAAAACGGAAAAGGTGTTCACAAAGGTCAATAAACGAGAAACCGTCCCAAGCCTa                                                            >  2:7858579/1‑150 (MQ=255)
                                                                                              aGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGcc                                                      <  2:3000654/150‑1 (MQ=255)
                                                                                                   atgaGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACg                                                 <  1:8382254/150‑1 (MQ=255)
                                                                                                   atgaGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTGCCAGTTGCAGCGAGTGTATTACAGCCAAGAATTCAAGGGAAAAGGTGGTGTCTAAGGTCATTAATAGTGAAACGTTAGCAAGCGTACAAGcccaaag                                                 >  1:712860/1‑146 (MQ=255)
                                                                                                     gaGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTGCCAGTTGCAGCGCGTGTATGACAGCCTAGAATTCAACGGAAATGGTGTTCTCTCAGGTCATTAATCGTGAAACGTTCGCAACCGTACACCCCCTACGgc                                               >  2:6716758/1‑150 (MQ=255)
                                                                                                        tgtTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCAGTTTTGCCCGTTGCAGCGAGTGTATGACAGCCGAGCATTCAAGGGAAATGGGGGTCTCGAAGGTCAATAATCGTGAAAACGTAGCAAGCGTACACGCCCAACCgcgcg                                            >  2:6371456/1‑149 (MQ=255)
                                                                                                             aTGAGCGAATCAACGTTTCGTCGGGCTGAAACTTGCATTTTTGCAAGTTGCAGCGAGTTTATGACAGCCTAGAAATCAAGGGAAATGGTGTTCTAAAAGGTCATTAATACTGAAAAGTTCGCAATCGTAAACCCCCAGCGGCGCGAGTTg                                       >  2:7562756/1‑150 (MQ=255)
                                                                                                              tGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTgg                                      <  2:888077/150‑1 (MQ=255)
                                                                                                                aGCGAATCAACGTTTCGTCGGGCTGCAAAAGGCATTTTTTCCAGTTGGAGCGAGTGTATGACAGCCTAGAATTCAAGGGAACTGGGGTTATCGAAGGTCATTAAACGTGAAAAGTTAGCAAGCGTACACGCCCTAAGGCGCGATTTgggg                                    >  2:3086764/1‑150 (MQ=255)
                                                                                                                aGCGAATCAACGTTTCGTCGGGCTGCAAAAGGCATTTTTGGCAGTTGCAGCGAGTGTATTACAGCCTAGAATTCAAGGGAAAATGTGTTGTGTAAGGTCAATAATAGTGAAAAGTTAGCAAGCGTACAAGCCCAACCGCGCGAGTTtggg                                    >  1:5183207/1‑150 (MQ=255)
                                                                                                                 gCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTTCAGCGAGTGTATGACCGCCGAGAAATCAAGGGAAATGGTGTTCTCTAAGGTCATTAATAGGGAAACGTTCGCAAGCGTACACGCCCAACGGCGCTATTTGGGGt                                   >  1:9154926/1‑150 (MQ=255)
                                                                                                                 gCGAATCAACGTTTCGTCGGGCTGCAACATGCATTTTTTCCAGTTTCAGCGAGTGTATTGCACCCTAGAAATAAAGGGAAATGGTGTTGTCTAAGGTCATTAATAGTGAAACGTTCGCAAGCGTAAACGCCCAACCGCGCTGTTTGGGGt                                   >  2:2565901/1‑150 (MQ=255)
                                                                                                                      tCAACGTTTCGTCGGGCTACAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATgg                              >  1:1592075/1‑150 (MQ=255)
                                                                                                                          cGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCAt                          <  2:4564613/150‑1 (MQ=255)
                                                                                                                              tcgtcgGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCt                      <  2:4639401/150‑1 (MQ=255)
                                                                                                                               cgtcgGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTc                     <  1:1839671/150‑1 (MQ=255)
                                                                                                                                    ggCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCgtgt                <  2:3576784/150‑1 (MQ=255)
                                                                                                                                     gCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGtt               <  1:3181315/150‑1 (MQ=255)
                                                                                                                                      cTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTg              <  1:4637740/150‑1 (MQ=255)
                                                                                                                                       tGCAACTGGCAAAAATGCCAGTTGAAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGTTCATTAATCGTGAAACTTTAGCAAGAGTACACGCCCTACGGCGCTATTTTGGTTGATGGCCATTGCTCCGTGTTGa             >  1:8481077/1‑150 (MQ=255)
                                                                                                                                        gCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGaa            <  1:554302/150‑1 (MQ=255)
                                                                                                                                         cAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAg           <  1:9153813/150‑1 (MQ=255)
                                                                                                                                          aaCTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTggg                                     >  2:9137313/1‑123 (MQ=255)
                                                                                                                                          aaCTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGt          >  1:8124158/1‑150 (MQ=255)
                                                                                                                                          aaCTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGt          <  1:7107621/150‑1 (MQ=255)
                                                                                                                                          aaCTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGt          >  1:4272006/1‑150 (MQ=255)
                                                                                                                                             tGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACa       >  2:3919467/1‑150 (MQ=255)
                                                                                                                                             tGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACa       >  2:6808212/1‑150 (MQ=255)
                                                                                                                                             tGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACa       >  2:1180536/1‑150 (MQ=255)
                                                                                                                                             tGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACa       >  2:9031893/1‑150 (MQ=255)
                                                                                                                                              ggCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTcc                    <  2:1064848/136‑1 (MQ=255)
                                                                                                                                                cATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATg    <  2:7641771/150‑1 (MQ=255)
                                                                                                                                                  tttttGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGt  <  2:2014201/150‑1 (MQ=255)
                                                                                                                                                  tttttGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGt  <  2:1338094/150‑1 (MQ=255)
                                                                                                                                                  tttttGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGt  >  2:8552709/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCGCAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGT  >  NC_002947/3580247‑3580542

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