Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,955,252 T→A 11.8% intergenic (+39/‑154) PP_1752 → / → PP_1753 hypothetical protein/DMT superfamily permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,955,2520TA11.8% 214.9 / 18.4 76intergenic (+39/‑154)PP_1752/PP_1753hypothetical protein/DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base T (36/31);  new base A (4/5);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01

CCCCTTGGACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGA  >  NC_002947/1955107‑1955398
                                                                                                                                                 |                                                                                                                                                  
ccccTTGGACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTAcc                                                                                                                                                <  2:876292/150‑1 (MQ=255)
ccccTTGGACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTAcc                                                                                                                                                <  1:418794/150‑1 (MQ=255)
 cccTTGGACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTAccc                                                                                                                                               <  2:680528/150‑1 (MQ=255)
 cccTTGGACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTAccc                                                                                                                                               >  2:1158482/1‑150 (MQ=255)
  ccTTGGACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCcg                                                                                                                                              >  2:915011/1‑150 (MQ=255)
  ccTTGGACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCcg                                                                                                                                              >  1:50633/1‑150 (MQ=255)
      ggACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTAAACCCGCGaa                                                                                                                                          <  1:1102493/150‑1 (MQ=255)
        aCGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGACAg                                                                                                                                        <  1:915011/150‑1 (MQ=255)
         cGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAgg                                                                                                                                       <  2:50633/150‑1 (MQ=255)
          gCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTCGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAggg                                                                                                                                      >  1:1083736/1‑150 (MQ=255)
           cTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTGCTGACCCTTTCGCGGGTAAACCCGCGaagcggc                                                                                                                                     <  1:848438/150‑6 (MQ=37)
            tggtggAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATATCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTc                                                                                                                                    <  1:1242311/150‑1 (MQ=255)
             ggtggAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCa                                                                                                                                   <  1:1104034/150‑1 (MQ=255)
                    gCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAgg                                                                                                                            >  1:860748/1‑150 (MQ=255)
                    gCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAgg                                                                                                                            >  2:764726/1‑150 (MQ=255)
                       aaTGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCaa                                                                                                                         <  2:860748/150‑1 (MQ=255)
                             aaGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTc                                                                                                                   >  1:1419175/1‑150 (MQ=255)
                             aaGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTc                                                                                                                   >  1:997606/1‑150 (MQ=255)
                                cgcTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCAc                                                                                                                <  1:764726/150‑1 (MQ=255)
                                    gggggggggTGAAGGTGGTGCGAAAAATAAACGAGGCCCGTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAg                                                                                                            <  2:1083736/144‑1 (MQ=255)
                                          gggTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTAAACCCGCGaa                                                                                                                                          <  2:530070/114‑1 (MQ=255)
                                          gggTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTAAACCCGCGaa                                                                                                                                          >  1:530070/1‑114 (MQ=255)
                                                 ggtggtGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCaaaa                                                                                               <  1:1158482/150‑1 (MQ=255)
                                                      tgCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGGCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATg                                                                                          >  1:758038/1‑150 (MQ=255)
                                                         gcACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTTGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGAcc                                                                                       <  1:246541/150‑1 (MQ=255)
                                                         gcACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGAcc                                                                                       >  2:246541/1‑150 (MQ=255)
                                                         gcACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGAc                                                                                        >  1:908972/1‑149 (MQ=255)
                                                         gcACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGAc                                                                                        <  2:908972/149‑1 (MQ=255)
                                                          cagaaTCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAgg                                                                                                                                       <  2:1302788/98‑1 (MQ=255)
                                                          cacaATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAgg                                                                                                                                       >  1:1302788/1‑101 (MQ=255)
                                                          cacaATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCt                                                                                      >  2:1393938/1‑150 (MQ=255)
                                                           acaATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCtg                                                                                     >  2:458144/1‑150 (MQ=255)
                                                             aaTCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCtgat                                                                                   <  1:1393938/150‑1 (MQ=255)
                                                             aaTCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCtgat                                                                                   >  1:628099/1‑150 (MQ=255)
                                                             aaTCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCtgat                                                                                   <  1:458144/150‑1 (MQ=255)
                                                                   ccAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTc                                                                             >  2:92340/1‑150 (MQ=255)
                                                                         gcgtGAGGCGGGTTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTTAGCACAGGCAACACCTCCACCCAGCCAGGCCCCAAAAGAATGACCTGAGGAATTCTTTtct                                                                       <  1:445673/147‑1 (MQ=255)
                                                                                     gTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGCATGACATGATGAATTCTTTTGTAGTCGCtctcgg                                                           >  2:442051/1‑148 (MQ=255)
                                                                                         gCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGACACTCTCCACCGt                                                       >  1:1313051/1‑150 (MQ=255)
                                                                                           cAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTTGCTCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGtct                                                     <  1:74138/150‑1 (MQ=255)
                                                                                            aaGCACAAGAAATAGCCGAAGGTTGGGTTGCCTGTGCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGtctc                                                    >  1:32656/1‑150 (MQ=255)
                                                                                             aGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTctct                                                   <  1:1270753/150‑1 (MQ=255)
                                                                                                  aaGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGc                                                                                                           >  1:1267269/1‑89 (MQ=255)
                                                                                                  aaGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGc                                                                                                           <  2:1267269/89‑1 (MQ=255)
                                                                                                   aGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGAACAGGCAAAACCTCCACCCAGCCATGCACCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCAACGTCTCTGtcctc                                             >  2:439061/1‑150 (MQ=255)
                                                                                                    gAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCtt                                                                           <  2:43870/119‑1 (MQ=255)
                                                                                                    gAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCtt                                                                           >  1:43870/1‑119 (MQ=255)
                                                                                                           ccGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCtt                                                                           <  1:759626/112‑1 (MQ=255)
                                                                                                           ccGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCtt                                                                           >  2:759626/1‑112 (MQ=255)
                                                                                                            cGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAAccccat                                    >  1:967311/1‑148 (MQ=255)
                                                                                                                ggTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGaa                                >  1:231032/1‑150 (MQ=255)
                                                                                                                      gTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTcaac                          >  2:732841/1‑150 (MQ=255)
                                                                                                                       ttGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGtc                                                <  1:364345/127‑1 (MQ=255)
                                                                                                                       ttGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGtc                                                >  2:364345/1‑127 (MQ=255)
                                                                                                                          ccTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAAcc                      <  1:732841/150‑1 (MQ=255)
                                                                                                                             gTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAAccccc                   >  1:609006/1‑150 (MQ=255)
                                                                                                                                cTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGcccc                                                                                                   <  1:1056085/67‑1 (MQ=255)
                                                                                                                                cTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGcccc                                                                                                   >  2:1056085/1‑67 (MQ=255)
                                                                                                                                  ggCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCcaccac              >  2:712171/1‑150 (MQ=255)
                                                                                                                                  gaccctttCGCGGGTAAACCCGCGAAGGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCcaccac              >  2:612728/7‑150 (MQ=255)
                                                                                                                                   gCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAAcccc                                      >  1:1280286/1‑125 (MQ=255)
                                                                                                                                   gCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAAcccc                                      <  2:1280286/125‑1 (MQ=255)
                                                                                                                                   gCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACg             >  2:283755/1‑150 (MQ=255)
                                                                                                                                    ccGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCATAGAATTCAACAACCCCCACCACgg            <  2:967311/150‑1 (MQ=255)
                                                                                                                                     cctttCGCGGGTAAACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGcc                                                                                                     >  2:1148005/4‑60 (MQ=25)
                                                                                                                                     cctttCGCGGGTAAACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGcc                                                                                                     <  1:1148005/57‑1 (MQ=25)
                                                                                                                                     cGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACggg           >  1:530177/1‑150 (MQ=255)
                                                                                                                                      gCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGACCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGc          >  1:1144548/1‑150 (MQ=255)
                                                                                                                                       tttCGCGGGTAAACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACTCCGTAGAATTCAACAAcc                      <  2:1128905/136‑1 (MQ=255)
                                                                                                                                       tttCGCGGGTAAACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAAcc                      >  1:1128905/2‑137 (MQ=255)
                                                                                                                                        ttCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAAcc                      >  1:886096/1‑136 (MQ=255)
                                                                                                                                        ttCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAAcc                      <  2:886096/136‑1 (MQ=255)
                                                                                                                                           gcgGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTa                                   <  1:874596/120‑1 (MQ=255)
                                                                                                                                           gcgGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTa                                   >  2:874596/1‑120 (MQ=255)
                                                                                                                                              gggttACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGCCCTGATCAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCCCCAACCCCGTAAAATTCACCAACCCCCGCCACGGGCTTGATCGa  <  2:628099/147‑1 (MQ=255)
                                                                                                                                              ggTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGa  <  2:231032/150‑1 (MQ=255)
                                                                                                                                              ggTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGa  <  1:210256/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CCCCTTGGACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGA  >  NC_002947/1955107‑1955398

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