Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,929,928 T→A 27.5% intergenic (‑1011/‑157) USA300HOU_1794 ← / → USA300HOU_1796 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,929,9280TA27.5% 284.7 / 130.7 156intergenic (‑1011/‑157)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (63/50);  new base A (24/19);  total (87/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

TATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCT  >  CP000730/1929790‑1930064
                                                                                                                                              |                                                                                                                                               
   ttttGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGCCATCAGTGGCTGTGCCATTGCGTTTTTTTcgt                                                                                                                                            <  2:895126/140‑3 (MQ=255)
      cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAACACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                                        >  2:518595/5‑141 (MQ=255)
        ccGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTgg                                                                                                                                       >  1:1519713/1‑141 (MQ=255)
         agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                                                                                     >  1:303314/2‑141 (MQ=255)
          gatgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCt                                                                                                                                     <  2:199270/141‑1 (MQ=255)
             gatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAg                                                                                                                                  <  1:753997/141‑1 (MQ=255)
              gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                                                                >  2:863701/2‑141 (MQ=255)
              atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                                                                                                 >  1:562697/1‑141 (MQ=255)
               tataATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGGAGAGAGGAGGTGTATGAGGTGATGCTTATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGCTGTGACATTGCGTTTTTTTCTTATTGGCTAAGTa                                                                                                                                <  2:319493/141‑1 (MQ=255)
               tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                               >  2:748652/1‑140 (MQ=255)
               tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                               <  2:260133/141‑2 (MQ=255)
                  aataatAGAGACGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                                                                            >  2:578323/1‑141 (MQ=255)
                   ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                           <  2:1369602/137‑1 (MQ=255)
                   ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACg                                                                                                                            >  2:816708/1‑141 (MQ=255)
                    tagtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGTCTAAGTAGACGc                                                                                                                           >  2:650653/1‑141 (MQ=255)
                    taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                                          >  2:827990/4‑141 (MQ=255)
                     agtagAGTCGCCTATCTCGCAGGCGTCAAATTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCa                                                                                                                          >  2:473945/1‑141 (MQ=255)
                      atagAGTCGCCTATCTCTCAGGCGTCAAGT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                                                                                        <  2:1651970/140‑2 (MQ=255)
                       tagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAt                                                                                                                        >  1:393789/1‑141 (MQ=255)
                        agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                      >  1:483866/1‑138 (MQ=255)
                         gagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑Ac                                                                                                                    <  1:1643263/141‑1 (MQ=255)
                          agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                    >  1:659176/1‑136 (MQ=255)
                          agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                    <  1:1173820/141‑6 (MQ=255)
                          agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                    <  1:333207/141‑6 (MQ=255)
                           gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACaa                                                                                                                  <  1:1601841/141‑1 (MQ=255)
                           gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                                   >  2:1063340/1‑135 (MQ=255)
                           gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                                   <  1:1116210/141‑7 (MQ=255)
                           gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                                   >  1:569305/1‑135 (MQ=255)
                            tCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACaa‑‑a                                                                                                               <  1:1464183/141‑1 (MQ=255)
                              gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGCGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AAt                                                                                                             >  2:1102600/1‑141 (MQ=255)
                              gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AAt                                                                                                             >  2:1452066/1‑141 (MQ=255)
                              gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AAt                                                                                                             >  2:1261408/1‑141 (MQ=255)
                              gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                <  1:1499478/141‑10 (MQ=255)
                              gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                <  1:107660/141‑10 (MQ=255)
                              gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                <  2:377700/141‑10 (MQ=255)
                                cTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AAAAg                                                                                                           >  1:828419/1‑141 (MQ=255)
                                cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                              <  2:1009074/141‑12 (MQ=255)
                                cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGGAAACGCaacaaataattaa                                                                                                              >  1:389714/1‑130 (MQ=255)
                                cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                              <  2:434342/141‑12 (MQ=255)
                                         aGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACatat                                                                                                 <  1:1058421/141‑1 (MQ=255)
                                           gCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAg                                                                                              >  2:398026/1‑141 (MQ=255)
                                              tCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAC‑‑AATAGGTAACATATACCCg                                                                                          >  2:593632/1‑141 (MQ=255)
                                               cAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGc                                                                                         >  1:234803/1‑141 (MQ=255)
                                                  tttAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCAcc                                                                                      <  2:82037/141‑1 (MQ=255)
                                                   ttAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCa                                                                                     <  1:616746/141‑1 (MQ=255)
                                                    tAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCaa                                                                                    >  2:307028/1‑141 (MQ=255)
                                                     aGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAAt                                                                                   <  1:1087759/141‑1 (MQ=255)
                                                      gACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATa                                                                                  >  2:398111/1‑141 (MQ=255)
                                                        cGCAGAGAGGAGGTGTATAAGTTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATaaa                                                                                <  2:148404/141‑1 (MQ=255)
                                                         gCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATaaaa                                                                               <  1:1004352/141‑1 (MQ=255)
                                                         gCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATaaaa                                                                               >  1:71364/1‑141 (MQ=255)
                                                         gCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTAATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATaaaa                                                                               >  2:1318084/1‑141 (MQ=255)
                                                                ggaggTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCt                                                                        >  2:1583716/1‑141 (MQ=255)
                                                                  aggTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCa                                                                      >  1:1449911/1‑141 (MQ=255)
                                                                        ataAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCg                                                                <  2:778769/141‑1 (MQ=255)
                                                                           aaGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCa                                                              >  1:482454/1‑141 (MQ=255)
                                                                           aaGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCATCACTACCGCa                                                              >  2:1164721/1‑141 (MQ=255)
                                                                              ggTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTatttgcgg                                                            <  1:433283/141‑8 (MQ=255)
                                                                              ggTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaaa                                                            <  1:1002215/141‑1 (MQ=255)
                                                                              ggTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaaa                                                            <  2:273451/141‑1 (MQ=255)
                                                                              ggTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaaa                                                            <  2:69780/141‑1 (MQ=255)
                                                                                    aTGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCCCTACCGCAAATAGt                                                        >  2:1552005/1‑141 (MQ=255)
                                                                                       cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACaa‑‑aa                                                                                                              <  1:1587018/87‑1 (MQ=255)
                                                                                         tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                   <  1:633162/141‑1 (MQ=255)
                                                                                         tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                   <  2:286852/141‑1 (MQ=255)
                                                                                          aTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                                  >  2:189264/1‑141 (MQ=255)
                                                                                           ttttCGTTCGCATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                 <  2:893844/141‑1 (MQ=255)
                                                                                           ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGTCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTATCGCAAATAGTGAGGGGa                                                 <  2:1382912/141‑1 (MQ=255)
                                                                                           ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                 >  1:1374596/1‑141 (MQ=255)
                                                                                           ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                 <  1:1315681/141‑1 (MQ=255)
                                                                                           ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                 <  2:1357559/141‑1 (MQ=255)
                                                                                             ttCGTTCACATCATAGCACCAGTAATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAtt                                               <  1:180279/141‑1 (MQ=255)
                                                                                              tCGTTCACATCATAGCCCCATCCATCCGTGGTTTTGCCATGTCGGTGGTTTTTTATTTTCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTg                                              <  1:1562916/141‑1 (MQ=255)
                                                                                                  tCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgtgt                                          <  1:1504494/141‑1 (MQ=255)
                                                                                                     catcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGtata                                       >  1:1042125/1‑141 (MQ=255)
                                                                                                      atcatAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGCGGGGATTGgtg                                           >  2:1641040/1‑141 (MQ=255)
                                                                                                          tAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTaa                                  >  1:1625092/1‑141 (MQ=255)
                                                                                                          tAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTaa                                  >  1:670369/1‑141 (MQ=255)
                                                                                                          tAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTaa                                  >  1:266107/1‑141 (MQ=255)
                                                                                                             cACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAt                                                                                                                        <  2:1492794/59‑1 (MQ=255)
                                                                                                             cACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAt                                                                                                                        >  1:1492794/1‑59 (MQ=255)
                                                                                                             cACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATa                               >  2:1654052/1‑141 (MQ=255)
                                                                                                              aCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATAc                              <  1:781169/141‑1 (MQ=255)
                                                                                                              aCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATAc                              <  1:359744/141‑1 (MQ=255)
                                                                                                              aCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATAc                              <  2:627953/141‑1 (MQ=255)
                                                                                                               ccAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACt                             <  1:785120/141‑1 (MQ=255)
                                                                                                                cAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTTTTTTTTTTTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACtt                            >  2:1303069/1‑141 (MQ=255)
                                                                                                                  gTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAt                          >  1:820262/1‑141 (MQ=255)
                                                                                                                     atcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtttt                       <  1:759162/141‑1 (MQ=255)
                                                                                                                      tccgTGGCTGTGCCATTGCGTTTTTTTCTTATTTGTTTAGTAGACGCGAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATTGTGAGGGAGATGGTGTAAAAGTATGTACGTAtgtgt                      >  2:745494/4‑137 (MQ=255)
                                                                                                                         gTGGCTGTGCCATTGCGTTTTTTTTTTTTTTTTTTAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCCCCAATAACAATCCCCTCACTACCGCAAATAGTTAGGGGATTGGTGTATACTTAAATTCTTTTTTTCGtt                   >  2:845157/1‑141 (MQ=255)
                                                                                                                         gTGGCTGTGCCATTGCGTTTTTTTTTTTTTTTTTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCttt                   >  2:921495/1‑141 (MQ=255)
                                                                                                                           ggCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGt                 >  1:1484354/1‑141 (MQ=255)
                                                                                                                            gCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGCGGATTGTTGTATAAGTAAATACTTATTTTCGTTGTc                <  1:793150/141‑1 (MQ=255)
                                                                                                                               gtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTa             >  2:769257/1‑141 (MQ=255)
                                                                                                                                tgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGc                                                                                             >  1:1185039/1‑61 (MQ=255)
                                                                                                                                tgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGc                                                                                             <  2:1185039/61‑1 (MQ=255)
                                                                                                                                tgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTaa            >  1:182714/1‑141 (MQ=255)
                                                                                                                                        cGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTg    <  2:1066434/141‑1 (MQ=255)
                                                                                                                                          tttttttCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCt  <  2:350340/141‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
TATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCT  >  CP000730/1929790‑1930064

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: