Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 724,813 A→T 82.7% intergenic (+17/‑226) ybfA → / → rhsC DUF2517 family protein/Rhs family putative polymorphic toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273724,8130AT82.7% 574.5 / 118.9 237intergenic (+17/‑226)ybfA/rhsCDUF2517 family protein/Rhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base A (20/21);  new base T (76/120);  total (96/141)
Fisher's exact test for biased strand distribution p-value = 2.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-01

Only 100 of 241 total aligned reads displayed.

CGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTC  >  CP009273/724666‑724962
                                                                                                                                                   |                                                                                                                                                     
cgcTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCa                                                                                                                                                     >  1:1116001/1‑150 (MQ=255)
cgcTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCa                                                                                                                                                     >  1:1047607/1‑150 (MQ=255)
  cTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCAtt                                                                                                                                                   <  2:1863465/150‑1 (MQ=255)
   tGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTGTTGAATTCATTg                                                                                                                                                  >  1:424420/1‑150 (MQ=255)
    gCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGa                                                                                                                                                 <  1:954110/150‑1 (MQ=255)
           aTGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCCGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTATTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGc                                                                                                                                          >  1:354568/1‑150 (MQ=255)
              cTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAg                                                                                                                                       <  1:4496969/150‑1 (MQ=255)
                     ggAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAAc                                                                                                                                >  1:2710939/1‑150 (MQ=255)
                     ggAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAAc                                                                                                                                <  2:3365722/150‑1 (MQ=255)
                      gAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCCGCGGAAACg                                                                                                                               >  2:779912/1‑150 (MQ=255)
                      gAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACg                                                                                                                               >  2:780209/1‑150 (MQ=255)
                      gAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACg                                                                                                                               <  1:1636596/150‑1 (MQ=255)
                         aGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACGGGTGATTTTTATTGAATGCTGAGTTTTGACTTCATAGAAAACCGCCCGCGGAAACGCTg                                                                                                                            >  2:2189624/1‑150 (MQ=255)
                           aCCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGc                                                                                                                          <  2:15224/150‑1 (MQ=255)
                            ccGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCTAAAACCAGCGATAAACCGGTGTGGATGGATCAGGACGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCg                                                                                                                         <  2:4419367/150‑1 (MQ=255)
                            ccGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGGGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCg                                                                                                                         >  2:887638/1‑150 (MQ=255)
                               cgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATACACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTGTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGtt                                                                                                                      >  2:2911969/1‑150 (MQ=255)
                               cgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGGGTGGATGGATCAGGCCGAAATGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGtt                                                                                                                      >  2:1782100/1‑150 (MQ=255)
                                   cGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGtttttt                                                                                                                  >  2:1112227/1‑150 (MQ=255)
                                    gCTTTGACAGCTACCTGCATCTCGTCTGGTCGAAAACCAGCGATAAACCGGTTTTGATGGATCAGGCCGAAAATGCAACTGTTTATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGTAAACGCTGGCGGTTTTTTg                                                                                                                 <  1:2198950/150‑1 (MQ=255)
                                      ttttACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTa                                                                                                               <  2:4485978/150‑1 (MQ=255)
                                       tttACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTAt                                                                                                              <  2:4152124/150‑1 (MQ=255)
                                        ttACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGttta                                                                                                             >  1:568868/1‑147 (MQ=255)
                                         tACAGCTACCTGAATCGCGTCTTGTCTAAAACCAGCGATAAACCGGTGTGTATGTTTCAGGCCGAAATGGCAACTGGTTTTTTTTATTGATTGCTGATTTTTTATTTCATTGAAACCCGCCTGCGGAAACGCTGGCGGTTTTTTGTATCa                                                                                                            <  1:459776/150‑1 (MQ=255)
                                               tACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCCGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGa                                                                                                      >  1:2859703/1‑150 (MQ=255)
                                               tACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTTTATCAGGAAGa                                                                                                      >  1:267317/1‑150 (MQ=255)
                                               tACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTGTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGa                                                                                                      >  1:2248854/1‑150 (MQ=255)
                                                aCCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGAt                                                                                                     <  2:3640698/150‑1 (MQ=255)
                                                 ccTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATg                                                                                                    <  2:333790/150‑1 (MQ=255)
                                                    gCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTc                                                                                                 <  2:2464491/150‑1 (MQ=255)
                                                     cATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCa                                                                                                <  2:4259820/150‑1 (MQ=255)
                                                           gTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGaa                                                                                          >  1:2601582/1‑150 (MQ=255)
                                                               ggTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATAtt                                                                                      >  2:4130433/1‑150 (MQ=255)
                                                                  cGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAg                                                                                   >  2:2000374/1‑150 (MQ=255)
                                                                   gAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGa                                                                                  <  2:3485105/150‑1 (MQ=255)
                                                                   gAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGa                                                                                  >  2:2539816/1‑150 (MQ=255)
                                                                   gAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGa                                                                                  >  2:272350/1‑150 (MQ=255)
                                                                   gAAAACCAGCGATAAACCGGAGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGa                                                                                  <  1:757900/150‑1 (MQ=255)
                                                                     aaaCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATAAGGAATATGTTCAGAGGAATATTTATACa                                                                                >  2:3147837/1‑150 (MQ=255)
                                                                       aCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATc                                                                              >  2:1403473/1‑150 (MQ=255)
                                                                         cAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGt                                                                            >  2:3212385/1‑150 (MQ=255)
                                                                           gCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGTTGATTTTTATTGATGGCTGATTTTTGAATTCATTGAAAACCGCCAGCGTAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTta                                                                          <  1:2758992/150‑1 (MQ=255)
                                                                           gCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTta                                                                          <  1:1380349/150‑1 (MQ=255)
                                                                            cGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTTAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTtat                                                                         <  1:620240/150‑1 (MQ=255)
                                                                             gATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATAGGtata                                                                        >  1:4106948/1‑150 (MQ=255)
                                                                              aTAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATAc                                                                       <  2:829544/150‑1 (MQ=255)
                                                                               tAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTTATTTTTATTGATTGCTGAGTTTTGAATTAATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATAcc                                                                      <  1:1256468/150‑1 (MQ=255)
                                                                                aaaCCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATAAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATAACa                                                                     >  2:3032105/1‑150 (MQ=255)
                                                                                   ccGGTGTGGATTGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATc                                                                  <  2:2695414/150‑1 (MQ=255)
                                                                                    cGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCg                                                                 <  2:166847/150‑1 (MQ=255)
                                                                                    cGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTAGTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCg                                                                 <  2:168475/150‑1 (MQ=255)
                                                                                     ggTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGa                                                                <  1:3602647/150‑1 (MQ=255)
                                                                                        gtgGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGAGTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAAc                                                             <  2:4448521/150‑1 (MQ=255)
                                                                                         tggatggaTCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACa                                                            <  2:3202985/150‑1 (MQ=255)
                                                                                         tggatggaTCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTTTATAAAGAAAATTTTAAGATGAATAATTATAAATAGTTATAACAAAAGaaaaa                                                            >  1:1983946/1‑148 (MQ=255)
                                                                                         tggatggaACAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACa                                                            <  2:3284083/150‑1 (MQ=255)
                                                                                           gatggatCAGGCCTAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACata                                                          <  1:3939856/150‑1 (MQ=255)
                                                                                              ggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatatat                                                       >  2:2122740/1‑150 (MQ=255)
                                                                                              ggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACAGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATCTTTAGACATCGTTATACCCATCGAAACatatat                                                       >  2:2122676/1‑150 (MQ=255)
                                                                                               gatCAGGCCTAAACGGCAACTGGTTATTTTTATTGATTGCTTTTTTTTGAATTCATTTAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATAtt                                                      <  1:789766/150‑1 (MQ=255)
                                                                                               gatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCATTCGAAACATATAtt                                                      <  2:3753003/150‑1 (MQ=255)
                                                                                               gatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTTTATCAGGAAGATGTTAAGAGGAATATTTAGACAAAGTTATAACAATCGAAACATATAtt                                                      >  1:1617785/1‑150 (MQ=255)
                                                                                               gatCAGGCAGAAAAGGCAACTGGTGATTTTTATGTATTGCTGAGTTTTGAATTAATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATAtt                                                      <  1:180479/150‑1 (MQ=255)
                                                                                                atCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGACTATTTAGACATCGTTATACCAATCGAAACATATATTc                                                     >  2:334541/1‑150 (MQ=255)
                                                                                                  cAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAATATTTTCAGAGGAATATTTAGACATCGTTATACCAATAGAAACATATAATCAt                                                   >  2:3179160/1‑150 (MQ=255)
                                                                                                  cAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCAt                                                   <  2:467815/150‑1 (MQ=255)
                                                                                                      ccGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGaaa                                               <  1:3249824/150‑1 (MQ=255)
                                                                                                      ccGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGaaa                                               <  1:4374485/150‑1 (MQ=255)
                                                                                                       cGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAat                                              <  1:1694557/150‑1 (MQ=255)
                                                                                                          aaaGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGCAGATGTTCAGAGGAATATTTAGAAATAGTTATAACAAAAGAAAAATAAATTAATTAAatata                                           >  1:2554765/1‑150 (MQ=255)
                                                                                                          aaaGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCCGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATGTTTAGACATCGTTATACCAATCGTACCATATATTAATGAAatata                                           >  2:510740/1‑150 (MQ=255)
                                                                                                          aaaGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGATGAATATTTAGACATCGTTATACCAATAGAAAAAAATATTCATGAAatata                                           >  2:2679664/1‑150 (MQ=255)
                                                                                                          aaaGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATAAGGAAGATGTTCAGAGGAATATTTATAAATCGTTATAACAATCGAAACATATATTCATTAAatata                                           >  2:415642/1‑150 (MQ=255)
                                                                                                            aGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAatatata                                         <  2:691347/150‑1 (MQ=255)
                                                                                                            aGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAatatata                                         <  2:4008503/150‑1 (MQ=255)
                                                                                                             ggCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAAAATATATTAATGAAatttatat                                        >  2:3950037/1‑150 (MQ=255)
                                                                                                               cAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATaa                                      <  1:3920158/150‑1 (MQ=255)
                                                                                                                  cTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAAAATTTAGACATCGTTAAACCAATCGAAACAAAAAATCATTAAATATATATAAata                                   >  1:284746/1‑150 (MQ=255)
                                                                                                                    ggTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGTAATATATATAAATAtt                                 >  2:2963329/1‑150 (MQ=255)
                                                                                                                     gTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATAttt                                <  2:2246983/150‑1 (MQ=255)
                                                                                                                      tGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAAAGCTGGCGGTTTTTTGTATAAGGTAGATTTTAAGACGAATAATTATACATAGTTATAACAAACCAAAAATATATTAATGAAATATATATAAAAAtttt                               >  1:1560739/1‑150 (MQ=255)
                                                                                                                       gATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTc                              <  2:1047607/150‑1 (MQ=255)
                                                                                                                        aTTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTcc                             <  2:1975502/150‑1 (MQ=255)
                                                                                                                         tttttATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCt                            >  2:595185/1‑150 (MQ=255)
                                                                                                                           tttATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTaa                          <  1:3802591/150‑1 (MQ=255)
                                                                                                                            ttATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAAt                         <  1:928253/150‑1 (MQ=255)
                                                                                                                             tATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATTTTCAGAGGAAAATTTAGACCTCGTTATACCAATCGAAACATATATTCATGAAATAAATATAAATATTTTCCTAAtt                        >  2:4322002/1‑150 (MQ=255)
                                                                                                                              attgattgCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTTTATCAGGAATTTGTTAAGAGGAATATTTATAAAAAGTTATACAAAAAGAAACAAATATTCATGAAAAATATATAAATATATTACAAAttt                       >  1:4476383/1‑149 (MQ=255)
                                                                                                                               ttgattgCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGt                      <  1:3311184/150‑1 (MQ=255)
                                                                                                                                tgattgCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATtttt                     >  2:3890813/1‑150 (MQ=255)
                                                                                                                                  attgCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCt                   <  2:1810379/150‑1 (MQ=255)
                                                                                                                                   tgggtGAGTTTTTAATTACTTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCTTTATACCTATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCtt                  <  1:143324/146‑1 (MQ=255)
                                                                                                                                    tgCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATAAGGAAGATTTTAAGATGAATAATTAGACATAGTTATACCAATCGAAAAAAATATTAATGAAATATATATAAATATTTTCATAATTTttatta                 >  1:3976982/1‑150 (MQ=255)
                                                                                                                                      cTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATc               >  2:3714878/1‑150 (MQ=255)
                                                                                                                                          gTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGAc           <  2:1244898/150‑1 (MQ=255)
                                                                                                                                            tttGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAg         <  1:780209/150‑1 (MQ=255)
                                                                                                                                               gAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGtta      >  2:1460799/1‑150 (MQ=255)
                                                                                                                                               gAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGata      >  2:2204037/1‑150 (MQ=255)
                                                                                                                                                 aTACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATAtc    >  2:638720/1‑150 (MQ=255)
                                                                                                                                                   aCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATAtctc  <  2:4138704/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTC  >  CP009273/724666‑724962

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